BZOJ1408 NOI2002 Robot 快速幂+欧拉函数

Posted WDZRMPCBIT

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ1408 NOI2002 Robot 快速幂+欧拉函数相关的知识,希望对你有一定的参考价值。

题意:分别求所有质因数都不同且质因数个数为奇数个、偶数个的数的欧拉函数和,和质因数存在重复的数的欧拉函数和

题解:

说书题……前面一大串就是用一种比较有趣的语言定义欧拉函数。
显然我们只需要求出所有军人ans1和政客ans2,然后用求得的M减去独立数和就是学者的独立数和了。
由于善良的出题人已经帮我们把M给质因分解了,因此我们假定当前已经得到了ans1和ans2,那么根据欧拉函数的定义,ans1\'=ans1+ans2*(p-1),ans2\'=ans2+ans1*(p-1),因为新加入的质因子会让奇偶互换。
由于e可能很大,所以用快速幂来算。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;

const int P=10000;
const int MAXN=1000+2;
int K,M=1,ans1=1,ans2,ans3;

int Quick_Pow(int x,int y){
    int t=(y&1?x:1);
    while(y>>=1){
        x=(x*x)%P;
        if(y&1) t=(t*x)%P;
    }
    return t;
}

int main(){
    cin >> K;
    for(int i=1,t1,t2,p,e;i<=K;i++){
        scanf("%d %d",&p,&e);
        M=(M*Quick_Pow(p,e))%P;
        
        if(p==2) continue;

        t1=(ans1+ans2*(p-1))%P,t2=(ans2+ans1*(p-1))%P;
        ans1=t1,ans2=t2;
    }
    ans3=(M-ans1-ans2+2*P)%P,ans1--;
    
    cout << ans1 << endl << ans2 << endl << ans3 << endl;

    return 0;
}
View Code

 

以上是关于BZOJ1408 NOI2002 Robot 快速幂+欧拉函数的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ1408[Noi2002]Robot DP+数学

bzoj 1408 [Noi2002]Robot(欧拉函数)

BZOJ_1408_[Noi2002]Robot_数学

1408: [Noi2002]Robot|快速幂|欧拉函数

BZOJ-1407: [Noi2002]Savage (扩展欧几里得)

bzoj1407 [Noi2002]Savage