码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ_1408_[Noi2002]Robot_数学

时间:2018-04-28 16:58:17      阅读:149      评论:0      收藏:0      [点我收藏+]

标签:alt   10个   esc   images   desc   bzoj   函数   efi   span   

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。


$\sum\limits_{d|n}\mu(d)=n$

因此总和为n。

只需要求约数中$\mu$为1的$\varphi$和,$\mu$为-1的$\varphi$和。

这样,我们每个质因子只有一次贡献。

这次贡献会把之前的$\mu$从1变到-1,从-1变到1。

又因为phi是积性函数,每次都需要乘上$\varphi(p)=p-1$。

要注意2不是奇质数。一个数的约数不考虑1。

 

代码:

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <stdlib.h>
using namespace std;
#define mod 10000
int n,m,p,k,f[1050],g[1050];
int qp(int x,int y) {
    int re=1;
    while(y) {
        if(y&1) re=re*x%mod;
        x=x*x%mod;
        y>>=1;
    }
    return re;
}
int main() {
    scanf("%d",&n);
    m=1;
    int i;
    f[0]=1;
    for(i=1;i<=n;i++) {
        scanf("%d%d",&p,&k); m=m*qp(p,k)%mod;
        if(p!=2) {
            f[i]=(f[i-1]+g[i-1]*(p-1)%mod)%mod;
            g[i]=(g[i-1]+f[i-1]*(p-1)%mod)%mod;
        }else {
            f[i]=f[i-1];
            g[i]=g[i-1];
        }
    }
    f[n]=(f[n]-1+mod)%mod;
    printf("%d\n%d\n%d\n",f[n],g[n],((m-f[n]-g[n]-1)%mod+mod)%mod);
}

 

BZOJ_1408_[Noi2002]Robot_数学

标签:alt   10个   esc   images   desc   bzoj   函数   efi   span   

原文地址:https://www.cnblogs.com/suika/p/8967140.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!