2018-2019 ACM-ICPC, Asia Dhaka Regional ContestC[数论]

题意很简单 求n! 的因子的因子的个数

首先对于任意一个数字 他的因子的个数等于将这个数字分解质因数后 每个质因数的次数加1 的乘积

我们将n!进行质因数分解之后 对于某个质因数p 假设他的最高次数为a 那么我们可以选择0次 1次2次3...... a次 对于每次选择 他的因子的个数就可以为次数加一 那么总的贡献就是 1+2+3+4+......+a+1= (a+1)*(a+2)/2;

那么所有的质因数的总贡献乘积就是答案 也就是

代码实现如下:

cpp 复制代码
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int mod=1e7+7;
const int N=1e6+5;
int qpow(int a,int b){
    int ans=1;
    for(;b;b>>=1){
        if(b&1)ans=ans*a%mod;
        a=a*a%mod;
    }
    return ans;
}
vector<int>minp(N,0),primes;
void sieve(int n){
    minp.assign(n+1,0);
    primes.clear();
    for(int i=2;i<=n;i++){
        if(!minp[i]){
            minp[i]=i;
            primes.push_back(i);
        }
        for(auto p:primes){
            if(i*p>n)break;
            minp[i*p]=p;
            if(p==minp[i])break;
        }
    }
}
signed main() {
    ios::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);
    sieve(N-4);
    int n;
    while(cin>>n&&n!=0){
        int ans=1;
        for(auto p:primes){
            if(p>n)break;
            int cnt=0;
            for(int j=n;j;j/=p)cnt+=j/p;
            ans=(ans*(cnt+1)%mod*(cnt+2)%mod)%mod *5000004%mod;
        }
        cout<<ans<<'\n';
    }
    return 0;
}
相关推荐
陈王卜13 分钟前
Reduce 算子优化笔记
算法
wabs66620 分钟前
关于图论【最短路径之Bellman_ford 算法(判断负权回路)|卡码网95.城市间货物运输II的思考】
数据结构·算法·图论·bellman_ford算法·卡码网·判断负权回路
YouonlyliveonceC34 分钟前
《大话数据结构》第4章实战:中缀表达式转后缀 + 后缀表达式求值(完整可运行实现)
算法
phoenix@Capricornus43 分钟前
从统计决策到贝叶斯估计
人工智能·算法·机器学习
AICDragon1 小时前
1.8%就够了:K3的896个MoE专家,为什么激活率这么低反而是好事?
人工智能·算法
VALENIAN瓦伦尼安教学设备1 小时前
ASHOOTER激光对中仪如何通过颜色确定调整是否合适
数据库·嵌入式硬件·算法
货拉拉技术1 小时前
重塑 Agent 度量衡:基于 LLM-as-a-Judge 的离线评估体系与实践
算法·设计模式
车压1 小时前
从“为什么”理解注意力机制
算法
wabs6662 小时前
关于哈希表【力扣15.三数之和的思考】
数据结构·算法·leetcode·散列表·哈希表·三数之和
退休倒计时2 小时前
【每日一题】LeetCode 88. 合并两个有序数组 TypeScript
算法·leetcode·职场和发展·typescript