hdu 3944 DP? (Lucas 定理)
生活随笔
收集整理的這篇文章主要介紹了
hdu 3944 DP? (Lucas 定理)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
仔細觀察楊輝三角后可以發現從最高點到第n行第k個數的最短路為c(n+1,k);
根據Lucas定理可以求出,一般來說要求答案模去一個質數p且p的范圍不大于10^5則可用Lucas.
Lucas(n,m,p)=cm(n%p,m%p)* Lucas(n/p,m/p,p) Lucas(x,0,p)=1; 另外注意當k>n/2時,必須令k=n-k。 #include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> using namespace std; const int maxn=10005; int n,k,p; int frc[1300][maxn]; int id[maxn]; bool vis[maxn]; int prim[maxn]; int cnt; void init(){for(int i=2;i<maxn;i++)if(!vis[i]){prim[cnt]=i;id[i]=cnt;for(int j=i*2;j<maxn;j+=i)vis[j]=1;p=i;frc[cnt][0]=1;for(int j=1;j<p;j++)frc[cnt][j]=frc[cnt][j-1]*j%p;cnt++;} } int POW(int x,int n,int p){int res=1;while(n){if(n&1)res=res*x%p;x=x*x%p;n>>=1;}return res; } int cm(int n,int m,int p){if(m>n)return 0;int x=id[p];int ans=frc[x][n]*POW(frc[x][n-m]*frc[x][m]%p,p-2,p)%p;return ans; } int lucas(int n,int m,int p){if(m==0)return 1;int ans=cm(n%p,m%p,p)*lucas(n/p,m/p,p)%p;return ans; } int main() { // freopen("in","r",stdin); init();int cas=1;while(scanf("%d%d%d",&n,&k,&p)>0){if(k>n/2)k=n-k;int ans=lucas(n+1,k,p)+n-k;printf("Case #%d: %d\n",cas++,ans%p);}return 0; }?
轉載于:https://www.cnblogs.com/wshh/p/3987407.html
總結
以上是生活随笔為你收集整理的hdu 3944 DP? (Lucas 定理)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: sql 系统函数取表基本信息
- 下一篇: 只有对经过自己努力才得到的东西,我才会高