程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
 程式師世界 >> 編程語言 >> C語言 >> C++ >> C++入門知識 >> poj 3696 The Luckiest number 歐拉函數在解a^x=1modm的應用

poj 3696 The Luckiest number 歐拉函數在解a^x=1modm的應用

編輯:C++入門知識

poj 3696 The Luckiest number 歐拉函數在解a^x=1modm的應用


題意:

給一個L,求長度最小的全8數滿足該數是L的倍數。

分析:

轉化為求方程a^x==1modm。之後就是各種數學論證了。

代碼:

 

//poj 3696
//sep9
#include 
#include 
using namespace std;
typedef long long ll;
ll L;
ll factor[65536];

ll mul(ll x,ll y,ll p)
{
	ll ret=0;
	while(y){
		if(y&1)
			ret=(ret+x)%p;
		x=(x+x)%p;
		y>>=1;
	}
	return ret;
}

ll pow(ll x,ll n,ll m)
{
	ll ret=1;
	x%=m;
	while(n){
		if(n&1)
			ret=mul(ret,x,m);
		x=mul(x,x,m);
		n>>=1;
	}
	return ret;
}

ll euler(ll n)
{
	ll ret=n;
	for(ll i=2;i*i<=n;++i)
		if(n%i==0){
			ret=ret/i*(i-1);
			while(n%i==0)
				n/=i;
		}
	if(n>1)
		ret=ret/n*(n-1);
	return ret;	
}

ll cal()
{
	ll m=L*9;
	for(int i=0;i<3;++i)
		if(m%2==0)
			m/=2;
		else
 			break;	
	if(m%2==0||m%5==0)
		return 0;
	ll phi=euler(m);
	int idx=0;
	for(ll i=1;i*i<=phi;++i)
		if(phi%i==0)
			factor[idx++]=i,factor[idx++]=phi/i;	
	sort(factor,factor+idx);
	for(int i=0;i

 

  1. 上一頁:
  2. 下一頁:
Copyright © 程式師世界 All Rights Reserved