POJ 2992 Divisors（约数定理）

最新推荐文章于 2019-04-02 20:53:09 发布

yphacker

最新推荐文章于 2019-04-02 20:53:09 发布

阅读量625

点赞数

CC 4.0 BY-SA版权

分类专栏： --------------数论-------------- 文章标签：约束定理

本文链接：https://blog.youkuaiyun.com/piaocoder/article/details/47955351

--------------数论-------------- 专栏收录该内容

39 篇文章

订阅专栏

本文介绍了解决POJ 2992问题的方法，通过预先计算1!到431!的素因子个数，并利用组合数C(n,k)的性质和求因子个数的公式，最终实现了高效求解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接：

http://poj.org/problem?id=2992

解题思路：

求C(n,k)的因子个数。。。

将1!~431!的每个素因子的个数打表，而C(n,k)=n!/(k!*(n-k)!)，所以再利用求因子个数的公式即可http://blog.youkuaiyun.com/piaocoder/article/details/47954385

AC代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;

typedef long long ll;
vector<int> v;
int prime[500];
int a[450][450];
int vis[500],nprime;
int sum[500];

void getprime(){
    nprime = 0;
    memset(vis,0,sizeof(vis));
    for(int i = 2; i <= 450; i++){
        int t = 450/i;
        for(int j = 2; j <=t; j++)
            vis[i*j] = 1;
    }
    for(int i = 2; i <= 450; i++){
        if(!vis[i])
            prime[nprime++] = i;
    }
}

void div(int x){
    v.clear();
    int t = x;
    for(int i = 0; i < nprime && prime[i]*prime[i] <= x; i++){
        while(t % prime[i] == 0){
            v.push_back(prime[i]);
            t /= prime[i];
        }
        if(t == 1)
            break;
    }
    if(t > 1)
        v.push_back(t);
}

void solve(){
    memset(a,0,sizeof(a));
    for(int i = 2; i <= 431; i++){
        memcpy(a[i],a[i-1],sizeof(a[i]));
        div(i);
        for(int j = 0; j < v.size(); j++)
            a[i][v[j]]++;
    }
}

int main(){
    getprime();
    solve();
    int n,k;
    while(~scanf("%d%d",&n,&k)){
        memset(sum,0,sizeof(sum));
        for(int i = 1; i <= 431; i++)
            sum[i] += a[n][i]-a[k][i]-a[n-k][i];
        ll ans = 1;
        for(int i = 2; i <= 431; i++)
            ans *= (sum[i]+1);
        printf("%lld\n",ans);
    }
    return 0;
}