對 N! 進行素因子分解：

n! = 1*2*3*4*…*n,因此它一定含有小于等于n 的所有素數我們要對這些素數分別進行處理求出含有多少項。

例如我們求n!含有多少個素因子p

設f[n][p] 表示n!含有多少素因子p；

1*2*3*4*,,,,,*n = ( p*2p*3p*4p*,,,,*(n/p*p) ) * k ? 其中k不含有素因子p；

(注： n/p*p ?= ?(int)(n/p) * p)

這里把 n 以內的所有 p 的倍數拿出來，剩下的數組成 k ，很明顯 k 里面不含有素因子p

1*2*3*4*…*n = p^(n/p) * (n/p)! ?* k ;

已知 k 里面不含有 素因子 p ，那么此時只需要求 (n/p)! 包含的素因子 p 個數，這里又回到了母問題 ，所以用遞歸解決。

則有遞歸方程：f[n][p] = n/p + f[n/p][p]

51nod 1189 :

#include

#include

using namespace std;

const long long mod = 2e9+14;//最后ans要除以2 ，所以這里 mod 取 1e9+7 的2倍

bool isprime[1000005];

int prime[100000],tot = 0,temp,n;

long long num;

long long getNum(int n,int p){ //用int可能會超出范圍

if(n

return n/p + getNum(n/p,p);

}

int main() {

memset(isprime,true,sizeof(isprime));

for(int i=2;i<=1000;i++) {

if(isprime[i]) prime[tot++] = i;

for(int j=i*i;j<=1000000;j+=i) isprime[j] = false;

}

for(int i=1001;i<=1000000;i++) {

if(isprime[i]) prime[tot++] = i; // prime 存儲 1000000 以內素數

}

long long ans;

while(scanf("%d",&n)!=EOF){

ans = 1;

for(int i=0;i

num = getNum(n,prime[i]);

ans *= (2*num+1);

ans %= mod;

}

if(ans%2==1) ans++;

ans %= mod;

printf("%I64d\n",ans/2);

}

}

總結

以上是生活随笔為你收集整理的c语言阶乘分解成素因子,阶乘的素因子分解 51nod 1189的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。