蒟蒻求助,为什么第11,12个点wa了
  • 板块P1593 因子和
  • 楼主Chaser2
  • 当前回复1
  • 已保存回复1
  • 发布时间2021/11/16 15:42
  • 上次更新2023/11/4 00:24:58
查看原帖
蒟蒻求助,为什么第11,12个点wa了
422992
Chaser2楼主2021/11/16 15:42
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

const int N =1e5+5, mod = 9901;

int p[N],s[N];

ll qpow(ll a, ll b)
{
    ll res=1;
    while(b)
    {
        if(b&1)res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}

int main()
{
    int a,b,cnt=0;
    cin>>a>>b;
    for(int i=2; i<=a/i; i++)
    {
        if(a%i==0)
        {
            p[++cnt]=i;
            while(a%i==0)
            {
                s[cnt]++;
                a/=i;
            }
        }
    }
    if(a>1)
    {
        p[++cnt]=a;
        s[cnt]++;
    }
    ll ans=1;
    for(int i=1; i<=cnt; i++)
    {
        //cout<<p[i]<<" "<<s[i]<<endl;
        ans=(ans*(qpow(p[i],s[i]*b+1)-1+mod)%mod*qpow(p[i]-1,mod-2)%mod+mod)%mod;
    }
    cout<<(ans+mod)%mod<<endl;
    system("pause");
    return 0;
}
2021/11/16 15:42
加载中...