POJ - 2992

Spring.A

已于 2022-01-30 10:08:22 修改

阅读量400

点赞数 1

分类专栏：数论文章标签：算法

于 2022-01-24 22:40:17 首次发布

本文链接：https://blog.csdn.net/weixin_60484917/article/details/122676394

版权

数论专栏收录该内容

8 篇文章 0 订阅

订阅专栏

参考题目：POJ - 2992

思路：打表(n!的质因数指数大小的表)

首先素数筛，把素数筛出来。
然后把每个阶乘的的质因子的指数存下来。使用累加，从1一直到431，可以把1! ~ 431! 的表全部打出来。原理：num[i][primes[j]] = num[i - 1][primes[j]] + count1 ，就是将上一个数字的指数累加到当前数上。
组合数公式： $C^k_n = n! / ((n - k)! * k!)$ 可以根据上一步打出来的阶乘数的质因子指数表，循环质数个数次得到答案

运行时间：

在这里插入图片描述

代码：

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

typedef long long LL;

const int N = 450;

int primes[N];
bool st[N];
int cnt;
int num[N][N];

int main()
{
    ios::sync_with_stdio(0); cin.tie(0);
    // 线性筛
    for (int i = 2; i <= 431; i ++ )
    {
        if (!st[i]) primes[cnt ++ ] = i;
        for (int j = 0; primes[j] <= 431 / i; j ++ )
        {
            st[i * primes[j]] = true;
            if (i % primes[j] == 0) break;
        }
    }
    
    // 打质因数指数表
    for (int i = 1; i <= 431; i ++ )
    {
        int num1 = i;
        for (int j = 0; j < cnt; j ++ )
        {
            int count1 = 0;
            while (num1 % primes[j] == 0)
            {
                count1 ++ ;
                num1 /= primes[j];
            }
            num[i][primes[j]] = num[i - 1][primes[j]] + count1;
        }
    }
    
    int n, k;
    while (cin >> n >> k)
    {
        LL sum = 1;
        
        for (int i = 0; i < cnt; i ++ )
        {
            LL up, down;
            up = num[n][primes[i]] - num[n - k][primes[i]];
            down = num[k][primes[i]];
            sum *= (up - down + 1);
        }
        
        cout << sum << endl;
    }
    
    return 0;
}