bzoj1408: [Noi2002]Robot

传送门
题目大半的篇幅在介绍$\phi$和$\mu$
不过不是普通的φ和μ，有些变形- -
新定义的φ(1)=0，新定义的μ只计算奇质数，含有2为因子的数都按照μ值为零处理
我们首先求出第一问和第二问，即μ值不等于0的部分
由于μ的定义，μ值不等于0当且仅当每个质因数的次数都是1次
因此我们枚举每个奇质数 计算加上这个奇质数之后φ值之和多出来的部分
由于φ是积性函数 所以前两问可以在O(n)时间内出解
第三问可以用总答案减掉前两问的答案
由于Σ[d|n]φ(d)=n，因此总答案就是m-1(此题1不算其它数的因数，去掉1)
用m-1减掉前两问的答案就是第三问的答案

#include<cstring>  
#include<cmath>   
#include<cstdio>  
#include<iostream>  
#include<cstdlib>   
#include<algorithm>
#define mo 10000
using namespace std;
int qp(int x,int y){
    if (y==1) return x;
    int t=qp(x,y/2);
    t=(t*t)%mo;
    if (y%2) t=(t*x)%mo;
    return t;
}
int n,m,x,y,u,v,ans1,ans2,ans3;
int main(){
    scanf("%d",&n);
    m=1;
    for (int i=1;i<=n;i++){
        scanf("%d%d",&x,&y);
        m=(m*qp(x,y))%mo;
        if (x==2) continue;
        u=(ans1+ans2*(x-1))%mo;
        v=(ans2+(ans1+1)*(x-1))%mo;
        ans1=u;
        ans2=v;
    }
    ans3=(m-1-ans1-ans2+3*mo)%mo;
    printf("%d\n%d\n%d\n",ans1,ans2,ans3);
}