BZOJ3884上帝与集合的正确用法-欧拉函数
生活随笔
收集整理的這篇文章主要介紹了
BZOJ3884上帝与集合的正确用法-欧拉函数
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
剛開始我想的是歐拉降冪,可是覺得復雜度還是挺高的就去找了一下題解。
思路大方向沒有問題,仍然是使用歐拉函數降低指數然后遞歸處理。但是不是簡單的使用歐拉降冪而是應該對模數p稍微處理一下。因為底數已經確定為2,所以我們可以將p寫成p=2k*q,q為奇數的形式,則
22的無窮次方%p=2k(22的無窮次方-k%q)
這樣可以直接對新的指數直接模q的歐拉函數值進行處理,因為每次q為奇數,q的歐拉函數值一定為偶數(他的因子一定是奇素數,奇素數的歐拉函數值一定是偶數),所以每次至少將p變為原來的一半,所以復雜度比較低。因為用到的歐拉函數值不是很多,所以直接求單點的歐拉函數值即可。
需要注意的是指數-k模歐拉函數值時應該保證指數大于0
AC代碼
#include<cstdio> #include<cstring> #include<algorithm> #include<climits> #include<cmath> #include<cstdlib> #include<ctime> #include<map> #include<queue> #include<set> #include<vector>using namespace std;typedef long long ll; const int MAXN=1e5+5;ll p;ll quick_pow(ll a,ll b,ll p) {a%=p; ll ret=1;while(b){if(b&1) ret=ret*a%p;a=a*a%p; b>>=1;}return ret; }ll quick_pow(ll a,ll b) {ll ret=1;while(b){if(b&1) ret*=a;a*=a; b>>=1;}return ret; }ll phi(int x) {ll ret=x;for(int i=2;i*i<=x;i++){if(x%i==0){ret=ret-ret/i;while(x%i==0) x/=i;}if(x==1) break;}if(x>1) ret=ret-ret/x;return ret; }ll deal(ll p) {ll k=0;while(~p&1) k++,p>>=1;if(p==1) return 0;ll phip=phi(p);return quick_pow(2,(deal(phip)-k%phip+phip)%phip,p)<<k; }int main() {int T;scanf("%d",&T);while(T--){scanf("%lld",&p);printf("%lld\n",deal(p));}return 0; }總結
以上是生活随笔為你收集整理的BZOJ3884上帝与集合的正确用法-欧拉函数的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: SPOJ-VLATTICE Visibl
- 下一篇: HDU1573-模线性方程