【分塊】【bitset】hdu6085 Rikka with Candies
阿新 • • 發佈:2017-08-08
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