题意:
求 组合数c(n,k)的因子数量
由算术基本定理很容易求得,不过第一次却T了,加了好多预处理,o1查询,才过
#include <iostream> #include <stdio.h> #include<string.h> #include<algorithm> #include<string> #include<ctype.h> using namespace std; #define I64d lld int prime[100]; int isnotprime[500]; int fac[500][100]; int nfac[500][100]; long long ans[500][500]; int np; void setprime() { np=0; for(int i=2;i<=431;i++) { if(!isnotprime[i]) { prime[np++]=i; } for(int j=0;j<np&&i*prime[j]<=431;j++) { isnotprime[i*prime[j]]=1; if(i%prime[j]==0) break; } } return; } void setfac() { memset(fac,0,sizeof(fac)); for(int i=2;i<=431;i++) { int p=i; for(int j=0;j<np;j++) { while(p%prime[j]==0) { fac[i][j]++; p/=prime[j]; } } } return; } void setans() { memset(nfac,0,sizeof(nfac)); for(int i=2;i<=431;i++) { for(int j=0;j<np;j++) nfac[i][j]=nfac[i-1][j]+fac[i][j]; } for(int n=0;n<=431;n++) { for(int k=0;2*k<=n;k++) { long long res=1; for(int i=0;i<np;i++) { res*=nfac[n][i]-nfac[k][i]-nfac[n-k][i]+1; } ans[n][k]=ans[n][n-k]=res; } } return; } int main() { setprime(); setfac(); setans(); int n,k; while(scanf("%d%d",&n,&k)!=EOF) { printf("%I64d\n",ans[n][k]); } return 0; }
时间: 2024-10-09 22:58:42