求m的所有约数中,满足可以分解成(奇数个不同素数/偶数个不同素数/其他)的所有的phi之和。
因此我们枚举每个奇质数 计算加上这个奇质数之后φ值之和多出来的部分
由于φ是积性函数 所以前两问可以在O(n)时间内出解
第三问可以用总答案减掉前两问的答案
由于Σ[d|n]φ(d)=n,因此总答案就是m-1(此题1不算其它数的因数,去掉1)
用m-1减掉前两问的答案就是第三问的答案
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
using namespace std;
inline int read()
{
int ans,f=1;char ch;
while ((ch=getchar())<'0'||ch>'9') if (ch=='-') f=-1;ans=ch-'0';
while ((ch=getchar())>='0'&&ch<='9') ans=ans*10+ch-'0';
return ans*f;
}
const int mod=10000;
int f[1005][2],p[1005],e[1005];
int a[1005];
int power_mod(int a,int b)
{
int ans=1;
while (b)
{
if (b&1) ans=ans*a%mod;
a=a*a%mod;
b>>=1;
}
return ans;
}
int main()
{
int n=read();
int m=1;
for (int i=1;i<=n;i++)
{
p[i]=read();e[i]=read();
m=m*power_mod(p[i],e[i])%mod;
if (p[i]!=2) a[++a[0]]=p[i];
}
m--;
f[1][0]=1;f[1][1]=a[1]-1;
for (int i=2;i<=a[0];i++)
{
f[i][0]=(f[i-1][0]+f[i-1][1]*(a[i]-1)%mod)%mod;
f[i][1]=(f[i-1][1]+f[i-1][0]*(a[i]-1)%mod)%mod;
}
f[a[0]][0]--;
printf("%d\n",(mod+f[a[0]][0])%mod);
printf("%d\n",(mod+f[a[0]][1])%mod);
printf("%d",(mod*2+m-f[a[0]][0]-f[a[0]][1])%mod);
return 0;
}