POJ 2992 质因子分解

简介:

题意要求出C(n,k)的因子个数 所以就需要把C(n,k)质因子分解 这题数据非常变态 我的思路没问题但是无数次TLE之后我才发现应该先打出一个 n!的一个表 不然就超时 用C(n,k)=n!/(k!*(n-k)!) 这个公式就可以求出

#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define max 435
int prime[max],nprime;
bool isprime[max];
void getprime()
{
    long long i,j;
    memset(isprime,1,sizeof(isprime));
    isprime[1]=0;
    nprime=0;
    for(i=2; i<max; i++)
        if(isprime[i])
        {
            prime[++nprime]=i;
            for(j=i*i; j<max; j+=i)
                isprime[j]=0;
        }
}
long long a[435][435];
void getnum()
{
    memset(a,0,sizeof(a));
    for(int s=1; s<=431; s++)
    {
        int x=s;
        memcpy(a[s],a[s-1],sizeof(a[s]));
        if(isprime[x])
        {
            a[s][x]++;
            continue;
        }
        for(int i=1; prime[i]<=x,x>1; i++)
            while(x%prime[i]==0&&prime[i]>0)
            {
                a[s][prime[i]]++;
                x/=prime[i];
            }
    }
}
int main()
{
    getprime();
    getnum();
    int n,k;
    long long ans;
    while(~scanf("%d%d",&n,&k))
    {
        ans=1;
        for(int i=2; i<=n; i++)
            if(a[n][i]-a[k][i]-a[n-k][i])
                ans*=(a[n][i]-a[k][i]-a[n-k][i]+1);
        printf("%lld\n",ans);
    }
    return 0;
}


目录
相关文章
|
6月前
POJ-2245-Lotto
POJ-2245-Lotto
28 0
|
容器
POJ 3640 Conformity
POJ 3640 Conformity
58 0
|
算法框架/工具
POJ 2262 Goldbach's Conjecture
POJ 2262 Goldbach's Conjecture
137 0
|
算法 数据建模 机器学习/深度学习
poj 3620
题意:给出一个矩阵,其中有些格子干燥、有些潮湿。       如果一个潮湿的格子的相邻的四个方向有格子也是潮湿的,那么它们就可以构成更大       的湖泊,求最大的湖泊。       也就是求出最大的连在一块儿的潮湿的格子的数目。
573 0
|
人工智能 BI
|
人工智能
POJ 1936 All in All
Description You have devised a new encryption technique which encodes a message by inserting between its characters randomly generated strings in a clever way.
791 0