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

HDU 4344

时间:2014-12-04 22:57:21      阅读:112      评论:0      收藏:0      [点我收藏+]

标签:blog   io   ar   os   sp   for   on   div   log   

其实不过是大整数分解。。。

注意两点:注意L 不能==N

但是,N却可以是素数。。。囧

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstdio>
#define LL __int64
#define MAX 1LL<<61
#define Times 8
#define N 501
#define C 101
using namespace std;
LL cop[N];
int ct;

bool cmp(LL a,LL b){
	if(a<b) return true;
	return false;
}

LL gcd(LL a,LL b){
	return b==0?a:gcd(b,a%b);
}

LL random(LL n){
	return (LL)((double)rand()/RAND_MAX*n+0.5);
}

LL multi(LL a,LL b,LL m){   //a*b%m
	LL ret=0;
	while(b){
		if(b&1)
		ret=(ret+a)%m;
		b>>=1;
		a=(a<<1)%m;
	}
	return ret;
}

LL quick(LL a,LL b,LL m){
	LL ans=1;
	a%=m;
	while(b){
		if(b&1){
			ans=multi(ans,a,m);
			b--;
		}
		b/=2;
		a=multi(a,a,m);
	}
	return ans;
}

bool Witness(LL a,LL n){
	LL m=n-1;
	int j=0;
	while(!(m&1)){
		j++;
		m>>=1;
	}
	LL x=quick(a,m,n);
	if(x==1||x==n-1)
	return false;
	while(j--){
		x=multi(x,x,n);
		if(x==n-1)
		return false;
	}
	return true;
}

bool Miller_Rabin(LL n){
	if(n<2) return false;
	if(n==2) return true;
	if(!(n&1)) return false;
	for(int i=1;i<=Times;i++){
		LL a=random(n-2)+1;
		if(Witness(a,n)) return false;
	}
	return true;
}

LL Pollard_Rho(LL n,int c){
	LL x,y,d,i=1,k=2;
	x=random(n-1)+1;
	y=x;
	while(true){
		i++;
		x=(multi(x,x,n)+c)%n;
		d=gcd(y-x,n);
		if(d>1&&d<n){
			return d;
		}
		if(y==x) return n;
		if(i==k){
			y=x;
			k=k<<1;
		}
	}
}

void find(LL n,int k){
	if(n==1) return ;
	if(Miller_Rabin(n)){
		cop[++ct]=n;
		return ;
	}
	LL p=n;
	while(p>=n)
	p=Pollard_Rho(p,k--);
	find(p,k);
	find(n/p,k);
}

int main(){
	int T;
	LL n;
	scanf("%d",&T);
	while(T--){
		scanf("%I64d",&n);
		ct=0;
		find(n,C);
		sort(cop+1,cop+1+ct,cmp);
		if(ct==1){
			printf("1 1\n");
			continue;
		}
		LL cnt=0;LL tmp=0; LL ans=0;
		cop[0]=1;
		cop[++ct]=1;
		for(int i=1;i<=ct;i++){
			if(cop[i]!=cop[i-1]){
				ans+=tmp;
				tmp=cop[i];
				cnt++;
			}
			else{
				tmp*=cop[i];
			}
		}
		if(ans==n)
		ans/=cop[1];
		printf("%I64d %I64d\n",cnt-1,ans);
	}
	return 0;
}

  

HDU 4344

标签:blog   io   ar   os   sp   for   on   div   log   

原文地址:http://www.cnblogs.com/jie-dcai/p/4143987.html

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