Codeforces1097 D. Makoto and a Blackboard（数论+dp）

题意：

解法：

$$\begin{gathered} 每次操作一定是某些质因子的个数减少或者不变. \\ 可以将每种质因子分开考虑,最后将每种质因子的贡献乘起来. \\ 首先O(\sqrt{n})质因子分解,然后分别考虑每种质因子: \\ 设当前质因子为p^{cnt}, \\ 令d[i][j]表示当前质因子进行i次操作后,剩下j个的概率, \\ 转移方程: \\ d[i][j]=\sum _{t=j}^{cnt}\frac{d[i-1][t]}{t+1} \\ 为什么分母是t+1个,因此剩下[0,t]的概率想等,区间长度为t+1, \\ dp边界为d[0][cnt]=1. \\ 当前质因子的贡献为\sum _{i=0}^{cnt}{p}^i\ast d[k][i]. \\ 最后将所有质因子的贡献乘起来就是答案. \end{gathered}$$

code：

#include<bits/stdc++.h>
#define int long long
#define PI pair<int,int>
using namespace std;
const int mod=1e9+7;
int d[11111][55];
int inv[55];
vector<PI>p;
int n,k;
int ppow(int a,int b,int mod){
    int ans=1%mod;a%=mod;
    while(b){
        if(b&1)ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
int cal(int p,int cnt){//质因子为p,个数为cnt
    for(int i=0;i<=k;i++){
        for(int j=0;j<=cnt;j++){
            d[i][j]=0;
        }
    }
    d[0][cnt]=1;
    for(int i=1;i<=k;i++){
        for(int j=0;j<=cnt;j++){
            for(int t=j;t<=cnt;t++){
                d[i][j]=(d[i][j]+d[i-1][t]*inv[t+1]%mod)%mod;
            }
        }
    }
    int ans=0;
    int base=1;
    for(int i=0;i<=cnt;i++){
        ans=(ans+base*d[k][i]%mod)%mod;
        base=base*p%mod;
    }
    return ans;
}
void solve(){
    for(int i=0;i<55;i++)inv[i]=ppow(i,mod-2,mod);
    cin>>n>>k;
    for(int i=2;i*i<=n;i++){
        if(n%i==0){
            int cnt=0;
            while(n%i==0)n/=i,cnt++;
            p.push_back({i,cnt});
        }
    }
    if(n!=1)p.push_back({n,1});
    int ans=1;
    for(auto i:p){
        ans=ans*cal(i.first,i.second)%mod;
    }
    cout<<ans<<endl;
}
signed main(){
    ios::sync_with_stdio(0);
    solve();
    return 0;
}

---