poj2992数论与组合数学,略水。。。

这题刚开始拿起来就果断用唯一分解定理,结果发现TLE。。。。。发现貌似数据超大,给做一个预处理。。。可是怎嘛预处理?不知道。后来看了别人的评论惊奇的发现前几天学过的一个定理用上了场,就是判断N!中有几个质因子:P(N!)=N/i+N/i^2+N/i^3+.....N/i^m(i为一个质因子,m是使N/i^m=0的最小值)

这样的化就可以预处理了,然后利用公式C(N,K)=N!/(K!*(N-K)!)和约数定理,就可以得到答案。。。。下面附代码:

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>

using namespace std;

int a[450],b[450],p[450][450];
int m;

int pj(int n,int x)
{
    int r=0,y=x;
    while (n>=y)
    {
          r+=n/y;
          y*=x;
    }
    return r;
}

int main()
{
    int n,k,i,j;
    for (i=2;i<=431;i++)
        if (!a[i])
        {
                  for (j=i*2;j<=431;j+=i) a[j]=1;
                  b[++m]=i;
        }
    for (i=1;i<=431;i++)
        for (j=1;j<=m;j++)
            p[i][j]=pj(i,b[j]);
    while (scanf("%d",&n)!=EOF)
    {
          scanf("%d",&k);
          __int64 r=1;
          for (i=1;i<=m;i++)
              if (n>=b[i])
              {
                 r*=(p[n][i]-p[n-k][i]-p[k][i]+1);
              }
          printf("%I64d\n",r);
    }
    return 0;
}

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值