bzoj1408（欧拉函数）

  求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;
}