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

POJ2635 The Embarrassed Cryptographer 简单数论

时间:2015-04-14 21:34:28      阅读:114      评论:0      收藏:0      [点我收藏+]

标签:

题目链接

看到这题的示意图也是醉了~题意:给你一个k,他是两个素数之积,然后给了一个数字L,然后找到具体是哪两个素数相乘等于k,较小的那个素数是否小于L,若小于L就输出 "BAD"外加较小的那个素数,否则就输出“GOOD”,

刚拿到这题目,有些钻牛角尖外加题意没看清楚,一开始纠结于 K很大,若想具体找出两个素数不可能,因为总有一个很大很大,求出其中一个素数 是否在10^6内是可以的,但是那时候我觉得还需要证明 k/prime  的值必须也是素数才符合要求,其实题目都规定好了 k必定是两个素数之积 而不会由其它合数相乘构成,所以就简单了许多

我不知道这个是什么定理的,但是对于一个很长的数字 比如  123456789,若把它每位上数字存在数组num中,123456789%m也等于

int now = 0;

for(int i=0;i<9;i++) now = (now * 10 + num[i])%m;

如果知道这个的话就好做了,对于K其实可以转化为大进制的数字存在数组里,比如转化为千进制万进制,先保存好,最后 计算的时候再乘以一千或者一万,上面那个for是相对于十进制的,把那个now*10改成对应进制即可,剩下的枚举素数只需要打印素数表即可,一开始害怕K太大存不下来而没有去细算直接转化为万进制,结果一直WA,后来转化为千进制就过了,原因是 for循环时每次取余的时候,余数最大可达10^6,然后再乘10000就会爆int,然后又改成了 long long,转化为万进制,也是可以过的

string s;

int L;

int nnum[100000 + 5];

int cnt ;

bool isprime[1000000 + 55];
int prime[1000000 + 55];
int k;

void make_prime() {
	memset(isprime,false,sizeof(isprime));
	for(int i=2;i<1000055 ;i++)
		if(!isprime[i])
			for(int j=i*2;j<1000055;j+=i)
				isprime[j]=true;
	for(int i=2;i<1000055;i++)
		if(!isprime[i])
			prime[k++]=i;
}

void init() {
	memset(nnum,0,sizeof(nnum));
	cnt = 0;
}

bool input() {
	while(cin>>s>>L) {
		if(s == "0" && L == 0)break;
		return false;
	}
	return true;
}

void slove() {
	int len = s.length();
	for(int i=len - 1;i >= 0;i-=3) {
		for(int j= i - 2;j<=i;j++) {
			if(j < 0)j = 0;
			nnum[cnt] = nnum[cnt] * 10 + s[j] - '0';
		}
		cnt++;
	}
}

bool gao(int x) {
	int ret = 0;
	for(int i= cnt - 1;i>=0;i--) {
		ret = (ret * 1000 + nnum[i])%x;
	}
	if(ret)return false;
	return true;
}

void cal() {
	slove();
	for(int i=0;i<k && prime[i] < L;i++) {
		if(gao(prime[i])) {
			printf("BAD %d\n",prime[i]);
			return ;
		}
	}
	puts("GOOD");
}

void output() {

}

int main() {
	make_prime();
	while(true) {
		init();
		if(input())return 0;
		cal();
		output();
	}
	return 0;
}



POJ2635 The Embarrassed Cryptographer 简单数论

标签:

原文地址:http://blog.csdn.net/yitiaodacaidog/article/details/45047373

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