組合數(shù)取模的話,之前多少會(huì)一些,能應(yīng)付一般的題目,而這次遇到了模數(shù)為合數(shù)的題目,于是就又來(lái)學(xué)習(xí)了一發(fā).
這次看到了一個(gè)比較不錯(cuò)的blog:https://blog.csdn.net/skywalkert/article/details/52553048
在這個(gè)blog里,其1.3里的內(nèi)容,有許多不理解的地方,并且3.2及以后的內(nèi)容,并沒有去研究.
這次主要是get到了用crt解決模數(shù)為合數(shù)的問題,并且還有與其配套使用的模數(shù)為質(zhì)數(shù)的冪的問題.
復(fù)習(xí)了一下crt,crt就是去按照一個(gè)正確且比較方便的方法去構(gòu)造一個(gè)解,并且利用了數(shù)在模里模外意義不同.
下面給出解決此類問題的代碼以及代碼注釋:
(此代碼對(duì)應(yīng)于具體題目,請(qǐng)讀者抓住重點(diǎn))

#include <cstdio> #include <cstring> #include <algorithm> typedef long long LL; const int N=1000010; int n,x,y,mod,P,p,phi,fac[N],num; inline int Pow(int x,int y){int ret=1;while(y){if(y&1)ret=(LL)ret*x%P;x=(LL)x*x%P,y>>=1;}return ret; } inline int cnt(int n){return n?n/p+cnt(n/p):0;//遞歸處理n的階乘里p的個(gè)數(shù) } inline int sum(int n){return n?((LL)(n/P?Pow(fac[P],n/P):1)*fac[n%P]%P*sum(n/p)%P):1;//遞歸處理n的階乘里刨去p之后,在模P的意義下的結(jié)果//為什么要遞歸呢?//我們發(fā)現(xiàn),我們預(yù)處理的時(shí)候,如果預(yù)處理的是刨去p的,那么他根本不循環(huán).//因?yàn)閷?duì)于一個(gè)含有p的數(shù)來(lái)說(shuō),他刨去p,和他加上P之后再刨去p,是不一樣的.//所以說(shuō),我們要預(yù)處理的是不含p的數(shù)的階乘,也就是說(shuō),如果一個(gè)數(shù)含有p,那就不乘他.//那么我們首先處理的是不含p的,然后去遞歸含有一個(gè)p的,然后是兩個(gè)的……//這樣答案就對(duì)了.//這樣好騷啊,一開始由于不知道這個(gè)操作而錯(cuò)*N……//不過(guò)這玩意log^2的吧…… } //以上兩個(gè)函數(shù)每次都尼瑪可以在一開始處理階乘的時(shí)候處理出來(lái)……然后就會(huì)跑得飛快…… //不過(guò)預(yù)處理的話,必須要求n較小,但是遞歸的話,只要n或P中的一個(gè)很小就可以了. #define deal(n,a,b) int a=sum(n),b=cnt(n); inline int C(int n,int m){if(m>n)return 0;deal(n,a0,b0)deal(m,a1,b1)deal(n-m,a2,b2)b0-=b1+b2;if(b0>=num)return 0;return (LL)a0*Pow(a1,phi-1)%P*Pow(a2,phi-1)%P*Pow(p,b0)%P; } inline int calc(){int ret=0,i,a,b,c,d;fac[0]=1;for(i=1;i<=P&&i<=n;++i)fac[i]=(i%p)?(LL)fac[i-1]*i%P:fac[i-1];a=0,b=x,c=(n-y-x)>>1,d=(n+y-x)>>1;while(c>=0){ret=(ret+(LL)C(n,a+b)*C(a+b,a)%P*C(c+d,c)%P)%P;++a,++b,--c,--d;}return ret; } int main(){//freopen("rio.in","r",stdin);scanf("%d%d%d%d",&n,&mod,&x,&y);x=std::abs(x),y=std::abs(y);if(n-y-x<0||((x&1)!=((n+y)&1))){puts("0");return 0;}int i,s=mod,ans=0;for(i=2;s>1;++i){if(i*i>s)i=s;if(s%i==0){p=i,P=1,num=0;while(s%p==0)++num,P*=p,s/=p;phi=P/p*(p-1);ans=(ans+(LL)calc()*(mod/P)%mod*Pow(mod/P,phi-1))%mod;}}printf("%d\n",ans);return 0; }

?

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

總結(jié)

以上是生活随笔為你收集整理的组合数取模学习笔记的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

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