bzoj1408 [Noi2002]Robot

http://www.elijahqi.win/archives/2953
Description

Input

Output

Sample Input

3
2 1
3 2
5 1

Sample Output

8
6
75
HINT

90号机器人有10个老师，加上它自己共11个。其中政客只有15号；军人有3号和5号；学者有8个，它们的编号分别是：2,6,9,10,18,30,45,90。

前两问根据欧拉函数是积性函数的性质 可以直接o(n)dp求出 对于第三问我们想办法求出总的这个欧拉函数 然后减去前面两个即可 我真的是sb 下午做题的时候刚复习莫比乌斯函数的一些公式就忘了$\varphi\times1=id$这个性质了

如何证明这个性质
因为$\mu\times id=\varphi(n)$
$\sum\limits_{d|n}\mu(d)*\frac{n}{d}$然后反演一下即可

#include<cstdio>
#include<cctype>
#include<algorithm>
using namespace std;
#define mod 10000
inline char gc(){
    static char now[1<<16],*S,*T;
    if(T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
inline int ksm(int b,int t){static int tmp;
    tmp=1;for (;t;b=b*b%mod,t>>=1) if (t&1) tmp=tmp*b%mod;return tmp;
}int k; 
int main(){
    freopen("bzoj1408.in","r",stdin);
    k=read();int dp1=0,dp2=0,ans=0,m=1;//dp1->odd  dp2->even
    for (int i=1;i<=k;++i){static int p,e;
        p=read();e=read();static int tmp1,tmp2;
        m=m*ksm(p,e)%mod;if (p==2) continue;
        tmp1=dp1;tmp2=dp2;
        (dp1+=(tmp2+1)*(p-1))%=mod;
        (dp2+=tmp1*(p-1))%=mod;
    }ans=((m-1-dp1-dp2)%mod+mod)%mod;
    printf("%d\n%d\n%d\n",dp2,dp1,ans);
    return 0;
}