日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問(wèn) 生活随笔!

生活随笔

當(dāng)前位置: 首頁(yè) > 编程资源 > 编程问答 >内容正文

编程问答

bzoj4589: Hard Nim

發(fā)布時(shí)間:2025/3/20 编程问答 21 豆豆
生活随笔 收集整理的這篇文章主要介紹了 bzoj4589: Hard Nim 小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

傳送門(mén)

不難看出就是求\(n\)個(gè)小于\(m\)的質(zhì)數(shù)異或和為\(0\)的方案數(shù),可以用\(FWT\)+快速冪解決
(我的代碼跑了4500ms……不是很明白那幾位52ms的巨巨是怎么做到的……可能是我人傻常數(shù)大……也不至于這么大吧……)

//minamoto #include<cstdio> #include<cstring> #define R register #define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i) #define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i) #define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v) using namespace std; const int N=(1<<16)+5,P=1e9+7,inv=500000004; inline int add(R int x,R int y){return x+y>=P?x+y-P:x+y;} inline int dec(R int x,R int y){return x-y<0?x-y+P:x-y;} inline int mul(R int x,R int y){return 1ll*x*y-1ll*x*y/P*P;} int vis[N],m,n,k; void init(){vis[1]=1;fp(i,2,50000)if(!vis[i]){for(R int j=(i<<1);j<=50000;j+=i)vis[j]=1;} } int A[N],B[N],lim; void FWT(int *A,int ty){for(R int mid=1;mid<lim;mid<<=1)for(R int j=0;j<lim;j+=(mid<<1)) for(R int k=0;k<mid;++k){int x=A[j+k],y=A[j+k+mid];A[j+k]=add(x,y),A[j+k+mid]=dec(x,y);if(ty==-1)A[j+k]=mul(A[j+k],inv),A[j+k+mid]=mul(A[j+k+mid],inv);} } void ksm(int *A,int *B,int y){FWT(A,1),FWT(B,1);while(y){if(y&1)fp(i,0,lim-1)B[i]=mul(B[i],A[i]);fp(i,0,lim-1)A[i]=mul(A[i],A[i]);y>>=1;}FWT(B,-1); } int main(){ // freopen("testdata.in","r",stdin);init();while(~scanf("%d%d",&n,&m)){lim=1;while(lim<=m)lim<<=1;memset(A,0,sizeof(A)),memset(B,0,sizeof(B));fp(i,1,m)A[i]=B[i]=!vis[i];ksm(A,B,n-1);printf("%d\n",B[0]);}return 0; }

轉(zhuǎn)載于:https://www.cnblogs.com/bztMinamoto/p/10195537.html

總結(jié)

以上是生活随笔為你收集整理的bzoj4589: Hard Nim的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。