【JZOJ B组】机器人M号

Description

3030年，Macsy正在火星部署一批机器人。 第1秒，他把机器人1号运到了火星，机器人1号可以制造其他的机器人。 第2秒，机器人1号造出了第一个机器人——机器人2号。 第3秒，机器人1号造出了另一个机器人——机器人3号。 之后每一秒，机器人1号都可以造出一个新的机器人。第m秒造出的机器人编号为m。我们可以称它为机器人m号，或者m号机器人。 机器人造出来后，马上开始工作。m号机器人，每m秒会休息一次。比如3号机器人，会在第6，9，12，……秒休息，而其它时间都在工作。 机器人休息时，它的记忆将会被移植到当时出生的机器人的脑中。比如6号机器人出生时，2，3号机器人正在休息，因此，6号机器人会收到第2，3号机器人的记忆副本。我们称第2，3号机器人是6号机器人的老师。 如果两个机器人没有师徒关系，且没有共同的老师，则称这两个机器人的知识是互相独立的。注意：1号机器人与其他所有机器人的知识独立（因为只有1号才会造机器人），它也不是任何机器人的老师。 一个机器人的独立数，是指所有编号比它小且与它知识互相独立的机器人的个数。比如1号机器人的独立数为0，2号机器人的独立数为1（1号机器人与它知识互相独立），6号机器人的独立数为2（1，5号机器人与它知识互相独立，2，3号机器人都是它的老师，而4号机器人与它有共同的老师——2号机器人）。 新造出来的机器人有3种不同的职业。对于编号为m的机器人，如果能把m分解成偶数个不同奇素数的积，则它是政客，例如编号15；否则，如果m本身就是奇素数或者能把m分解成奇数个不同奇素数的积，则它是军人，例如编号 3, 编号165。其它编号的机器人都是学者，例如编号2, 编号6, 编号9。 第m秒诞生的机器人m号，想知道它和它的老师中，所有政客的独立数之和，所有军人的独立数之和，以及所有学者的独立数之和。可机器人m号忙于工作没时间计算，你能够帮助它吗？ 为了方便你的计算，Macsy已经帮你做了m的素因子分解。为了输出方便，只要求输出总和除以10000的余数。

Input

输入文件的第一行是一个正整数k(1<=k<=1000)，k是m的不同的素因子个数。 以下k行，每行两个整数，pi, ei,表示m的第i个素因子和它的指数(i = 1, 2, …, k)。p1, p2, …, pk是不同的素数。所有素因子按照从小到大排列，即p1

Output

输出文件包括三行。 第一行是机器人m号和它的老师中，所有政客的独立数之和除以10000的余数。 第二行是机器人m号和它的老师中，所有军人的独立数之和除以10000的余数。 第三行是机器人m号和它的老师中，所有学者的独立数之和除以10000的余数。

Sample Input

3 2 1 3 2 5 1

Sample Output

8 6 75

Hint

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

思路

先理一下题意：

三个询问：
1. 能把 m 分解成偶数个不同奇素数的积的欧拉函数值
2. m 本身就是奇素数或者能把 x 分解成奇数个不同奇素数的积的欧拉函数值
3. 剩下数的欧拉函数值

对于（1）和（2），注意到乘上一个奇质数能使（1）（2）互相转化，不妨用f[i][0]表示讨论到第i个数时满足（2）的答案，f[i][1]表示讨论到第i个数时满足（1）的答案。考虑两种情况的递推关系。

这里要用到欧拉函数的一个性质:欧拉函数是积性函数，既当gcd(a,b)=1gcd(a,b)=1时，有ϕ(a)×ϕ(b)=ϕ(a×b)ϕ(a)×ϕ(b)=ϕ(a×b)
根据这个容易得到递推式：

f[i][0]=f[i−1][1]×(p[i]−1)+f[i−1][0]f[i][0]=f[i−1][1]×(p[i]−1)+f[i−1][0]
f[i][1]=f[i−1][0]×(p[i]−1)+f[i−1][1]f[i][1]=f[i−1][0]×(p[i]−1)+f[i−1][1]
(欧拉函数的性质，当p为素数时，ϕ(p)=p−1)(欧拉函数的性质，当p为素数时，ϕ(p)=p−1)
这样前两个问题都被解决了，考虑第三个问题。这又要用到欧拉函数的一个性质：

n=∑d|nϕ(d)n=∑d|nϕ(d)
这样第三个问题就可以解决了，注意不含1。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=1077,mod=10000;
int f[maxn][2],p[maxn],e[maxn],n,M;
int power(int a,int b)
{
    int ass=1;
    while(b)
    {
        if(b&1) ass=(ass*a)%mod;
        b>>=1; a=a*a%mod;
    }
    return ass;
}
int main()
{
    scanf("%d",&n);
    for(int i=1; i<=n; i++) scanf("%d%d",&p[i],&e[i]);
    f[0][0]=1;
    if(p[1]==2) f[1][0]=1;else f[1][1]=p[n+1]-1;
    for(int i=1; i<=n; i++)
    {
        if(p[i]==2)continue;
        f[i][0]=(f[i-1][1]*(p[i]-1)%mod+f[i-1][0])%mod;
        f[i][1]=(f[i-1][0]*(p[i]-1)%mod+f[i-1][1])%mod;
    }
    M=1;
    for(int i=1; i<=n; i++) 
    M=M*power(p[i],e[i])%mod;
    M=((M-f[n][0]-f[n][1])%mod+mod)%mod;
    printf("%d\n%d\n%d",f[n][0]-1,f[n][1],M);
}