Hello 2019 D Makoto and a Blackboard （积性函数+DP）

题目链接：http://codeforces.com/contest/1097

对于,

状态转移方程：。

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll unsigned long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<ll,ll>
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=2e4+9;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意：给定n和k，
初始时候是n，一次变换后其等概率的变为其因子之一，
问k次变换的期望值。

题目分析：假如已知E(n)，当n变为n*p且p与n互质的时候,
分析计算E（n）最后的权重和概率，假设x出现的概率是q,
那么x变为x*p的概率是1/2^k，保持不变的概率是(2^k-1)/2^k，
相当于乘上E(p),不难发现E(n*p)=E(n)*E(p)。
可知E(n)是积性函数，这样我们质因分解后得到质数p和幂次k，
可以分析E(p^k)，转移方程很好列，再加上幂次方最高是64，
k的范围是四次方，所以时间复杂度是合适的，
转移方程：E(k,s)=sigma i=0~k E(i,s-1)。
*/
ll n,k;
ll dp[maxn][80];
int main(){
    cin>>n>>k;
    ll ans=1LL;
    mst(dp,0);
    for(ll i=2;i*i<=n;i++) if(n%i==0){
        int cnt=0;
        while(n%i==0) n/=i,cnt++;
        dp[0][0]=1LL;rep(j,1,cnt+1) dp[0][j]=dp[0][j-1]*i%mod;
        rep(j,1,cnt+1) dp[0][j]=(dp[0][j]+dp[0][j-1])%mod;
        rep(j,1,k+1){
            dp[j][0]=dp[j-1][0]%mod;
            rep(p,1,cnt+1) dp[j][p]=dp[j-1][p]*powmod(p+1,mod-2)%mod;
            if(j!=k) rep(p,1,cnt+1) dp[j][p]=(dp[j][p]+dp[j][p-1])%mod;
        }
        ans=ans*dp[k][cnt]%mod;
    }
   /// cout<<ans<<endl;
    if(n>1){
        ll inv=powmod(2LL,mod-2);
        inv=powmod(inv,k);
        n=(n-1)%mod;
        inv=inv*n%mod;
        inv=(inv+1)%mod;
        ans=ans*inv%mod;
    }
    cout<<ans<<endl;
    return 0;
}