1. 程式人生 > >CODE VS 4939 尤拉函式 質因數啟發式分解

CODE VS 4939 尤拉函式 質因數啟發式分解

 

/**
CODE VS 4939 尤拉函式
質因數啟發式分解 
連結:http://codevs.cn/problem/4939/
*/

#include<bits/stdc++.h>
#define ll long long
using namespace std;

/**********************************************Head-----Template****************************************/
bool Finish_read;
template<class T>inline void read(T &x){Finish_read=0;x=0;int f=1;char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;if(ch==EOF)return;ch=getchar();}while(isdigit(ch))x=x*10+ch-'0',ch=getchar();x*=f;Finish_read=1;}
template<class T>inline void print(T x){if(x/10!=0)print(x/10);putchar(x%10+'0');}
template<class T>inline void writeln(T x){if(x<0)putchar('-');x=abs(x);print(x);putchar('\n');}
template<class T>inline void write(T x){if(x<0)putchar('-');x=abs(x);print(x);}
inline ll gcd(ll a,ll b){ return b==0?a:gcd(b,a%b);}
inline ll lcm(ll a,ll b){ll gg=gcd(a,b);a/=gg;if(a<=LLONG_MAX/b) return a*b;return LLONG_MAX;}
/********************************Head----Temlate**********************************************/


inline ll mul(ll a, ll b, ll mod){
	return a%=mod, b%=mod, (a*b-(long long)((long double)a*b/mod+1e-3)*mod+mod)%mod;
}

inline ll ksm(ll x, ll y, ll mod){
	ll ans=1;
	for (;y;y>>=1, x=mul(x, x, mod)) if (y&1) ans=mul(ans, x, mod);
	return ans;
}

/**********
Miller_Rabin 素數判定
利用費馬小定理,以prime為基,不以隨機函式進行測試;
利用a^(p-1)=1(mod p) 利用較大概率預設p為素數,也就是忽略了偽素數的情況
優化:提取p-1中2的倍數,p-1=2^(k)*t,先算出a^t mod p 
再進行不斷平方,暴力次數[0,k],直到中間過程存在ans = p-1 則所判數字為素數
*/
bool Miller_Rabin(ll n){
	static const ll prime[]={2, 3, 5, 7, 11, 13, 17, 19, 23};
	if (n<=23){ for (int i=0;i<9;i++) if (n==prime[i]) return true;return false;}
	if (~n&1) return false;
	for (int t=0;t<9;t++){
		ll y=n-1; int bit=0;
		for (;~y&1;y>>=1,bit++);
		ll x=ksm(prime[t],y,n);
		if (x==1) goto ed;//goto 減少標記次數
		for (;bit--;x=mul(x, x, n)) if (x==n-1) goto ed;//判處最後的ans
		return false;
		ed:;
	}
	return true;
}

vector<ll>vec;//儲存質因數分解結果(無序)
/*************
Pollard_Rho質因數分解:啟發式分解

時間複雜度:O(n^(1/4)),空間複雜度:lg(n)

首先確定一個數x為n的質因數,隨機找的話概率為1/p;
如果通過某種倍增的規律進行遞推x=(x*x+c)%n,可以得到兩個數x1,x2
利用gcd(abs(x1-x2),n)>1,使得x1 - x2 不斷向左式進行靠近
由此分解出來的差值為n的質因數的概率顯著增加
 
最後返回兩種結果:
op 1:因數
op 2 :x1 == x2 也就是說當前的倍增數c不滿足當前解,那麼就需要換個c進行嘗試;


換c進行嘗試的過程,必然很隨機,因得到的ans是存在重複且無序的

*/
void Pollard_Rho(ll n){
	if (Miller_Rabin(n)) return void (vec.push_back(n));
	for (ll c=1;;c++){
		for (ll i=1,x2=0,x1=c;;x1=mul(x1,x1,n)+c,i++){
			if (x1==x2) break;
			ll d=gcd(n, x1>=x2?x1-x2:x2-x1);
			if (1<d&&d<n) return Pollard_Rho(d),Pollard_Rho(n/d),void();
			if (i==(i&-i)) x2=x1;
		}
	}
}
int main(){
	for (;;){
		ll n;read(n); if(n==0) break;
		Pollard_Rho(n);
		sort(vec.begin(),vec.end()); int num = unique(vec.begin(),vec.end())-vec.begin();
		ll ans=n;
		for(int i=0;i<num;i++) ans=ans/vec[i]*(vec[i]-1);
		printf("%lld\n",ans);
		vec.clear();
	}
	return 0;
}