【分块】【bitset】hdu6085 Rikka with Candies

给你数组A和B，A B中的元素大小都不超过5w，且两两不同。

q次询问，每次给你个k，问你有多少对(i,j)，满足A(i)%B(j)==k。

如题目所言模拟bitset的过程，实质上是个分块，每块的大小定为63。

一个小技巧是对于最终的那个数组w，分块后记63个w数组，每个数组最前面一块是零散的部分，大小从1~63，这样比较好操作。

最后把63个w里面的每一位的值都异或起来，就是对应的k的答案。

#include<cstdio>

#include<cstring>

#include<algorithm>

using namespace std;

const int SZ=63;

typedef unsigned long long ull;

ull blo[805],b[67][805],lss[67];

int T,n,m,q,sum=794,l[805],r[805],num[50005],sz[805];

bool a[50005],c[50005];

int main(){

	int x;

	scanf("%d",&T);

	for(;T;--T){

		memset(a,0,sizeof(a));

		memset(blo,0,sizeof(blo));

		memset(b,0,sizeof(b));

		memset(lss,0,sizeof(lss));

		memset(l,0,sizeof(l));

		memset(r,0,sizeof(r));

		memset(num,0,sizeof(num));

		memset(sz,0,sizeof(sz));

		memset(c,0,sizeof(c));

		scanf("%d%d%d",&n,&m,&q);

		for(int i=1;i<=n;++i){

			scanf("%d",&x);

			a[x]=1;

		}

		r[0]=-1;

		for(int i=1;i<=sum;++i){

			l[i]=r[i-1]+1;

			r[i]=min(i*SZ-1,50000);

			sz[i]=r[i]-l[i]+1;

			for(int j=l[i],p=0;j<=r[i];++j,++p){

				if(a[j]){

					blo[i]|=((ull)1<<p);

				}

				num[j]=i;

			}

		}

		for(;m;--m){

			scanf("%d",&x);

			for(int ql=0;ql<=50000;ql+=x){

				int qr=min(ql+x-1,50000);

				if(num[ql]==num[qr]){

					int fls=ql-l[num[ql]];

					lss[qr-ql+1]^=((blo[num[ql]]>>fls)&(((ull)1<<(qr-ql+1))-(ull)1));

				}

				else{

					int ls=r[num[ql]]-ql+1;

					for(int i=num[ql]+1,j=1;i<num[qr];++i,++j){

						b[ls][j]^=blo[i];

					}

					int rs=qr-l[num[qr]]+1;

					lss[ls]^=(blo[num[ql]]>>(sz[num[ql]]-ls));

					b[ls][num[qr]-num[ql]]^=(blo[num[qr]]&(((ull)1<<rs)-(ull)1));

				}

			}

		}

		for(int i=1;i<=SZ;++i){

			for(int j=0;j<i;++j){

				c[j]^=((lss[i]>>j)&(ull)1);

			}

			int now=0,wei;

			for(int k=i,p=1;k<=50000;++k,++p,++wei){

				if(p%SZ==1){

					++now;

					wei=0;

				}

				c[k]^=((b[i][now]>>wei)&((ull)1));

			}

		}

		for(;q;--q){

			scanf("%d",&x);

			printf("%d\n",c[x]);

		}

	}

	return 0;

}