1. 程式人生 > >【分塊】【bitset】hdu6085 Rikka with Candies

【分塊】【bitset】hdu6085 Rikka with Candies

com die 過程 names -1 memset sign with pri

給你數組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;
}

【分塊】【bitset】hdu6085 Rikka with Candies