poj 2992 Cipher

最新推荐文章于 2021-03-30 21:38:04 发布

原创最新推荐文章于 2021-03-30 21:38:04 发布 · 375 阅读

0 ·

CC 4.0 BY-SA版权

数论专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一种质因数分解算法的实现，通过预先计算所有小于432的质数，并利用递归方法来计算特定整数范围内每个质数的出现次数，进而计算组合数。该算法适用于竞赛编程和其他需要高效质因数分解的应用场景。

质因数分解。

#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
int isPrime[432],pChart[432];
long long cntC[432][432];
int cot;
void getPrime()
{
    memset(pChart,0,sizeof(pChart));
    int i,j;
    for(i = 2;i < 432;i ++)
    {
        for(j = 2;j * i < 432;j ++)
        pChart[i*j] = 1;
    }
    j = 0;
    for(i = 2;i < 432;i ++) if(pChart[i] == 0) isPrime[j ++] = i;
    cot = j;
}
int cal(int n,int p)
{
    if(n < p) return 0;
    else return n / p + cal(n / p, p);
}
int main()
{
    int n,k,i,j;
    getPrime();
    //for(i = 0;i < cot;i ++) cout<<isPrime[i]<<endl;
    /*for(i = 1;i < 432;i ++)
    {
        for(j = 1;j < i;j ++)
        {
            cntC[i][j] = 1;
            for(k = 0;k < cot&&isPrime[k] <= i;k ++)
            cntC[i][j] *= (cal(i,isPrime[k]) - cal(i-j,isPrime[k]) -cal(j,isPrime[k]) +1);
        }
    }*/
    while(scanf("%d%d", &n, &k) != EOF)
    {
        __int64 cnt = 1;
        //if(k > n / 2) k = n - k;
        for(i = 0; isPrime[i] <= n && i < cot; i ++) {
            int t = cal(n, isPrime[i]) - cal(n - k, isPrime[i]) - cal(k, isPrime[i]);
            cnt *= (t + 1);
        }
        printf("%I64d\n", cnt);
    }
    return 0;
}