sicily1231 The Embarrassed Cryptography

最新推荐文章于 2020-07-30 03:54:59 发布

guozhchun

最新推荐文章于 2020-07-30 03:54:59 发布

阅读量490

点赞数

分类专栏： sicily 文章标签： sicily 1231 大数求模素数

本文链接：https://blog.csdn.net/guozhengchun/article/details/41412145

版权

sicily 专栏收录该内容

11 篇文章 0 订阅

订阅专栏

题目链接：http://soj.sysu.edu.cn/1231

题目大意：给出一个（大）整数K，K是两个素数的乘积，另外给出一个整数L，要求判断K是否能被小于 L 的素数整除

解题思路：大数求模

依据：（1）模的运算：(a * b + c) % d = ((a * b) % d + c % d) % d = ((a % d) * (b % d) + c % d) % d

（2）n位数的表示方法：a1a2a3...an = (((a1 * 10) + a2) * 10 + a3) + a4...) + an

#include <iostream>
#include <math.h>
#include <string>
using namespace std;

int prime[100000];			//  存储素数的数组
int ProdecePrime(void)		//  用于产生素数的函数，返回值是数组的长度
{
	int count = 0, i, j, tem;
	prime[count++] = 2;
	prime[count++] = 3;
	//  找出所有小于1000000的素数存于prime数组中
	for (i = 5; i < 1000000; i++)
	{
		tem = (int)sqrt(i);
		for (j = 2; j <= tem; j++)
			if (i % j == 0)
				break;
		if (j == tem + 1)
			prime[count++] = i;
	}
	return count;
}

int main()
{
	int L, count,length, mod;
	string s;
	bool flag;
	count = ProdecePrime();		//  产生素数并获取数组长度
	while (1)
	{
		cin >> s >> L;
		if (s == "0" && L == 0)
			break;
		length = s.size();
		flag = true;			//  用于标记是否应该输出GOOD，值为true时输出
		for (int i = 0; i < count; i++)
		{
			if (prime[i] >= L)
				break;
			mod = 0;
			for (int j = 0; j < length; j++)	//  大数求模
				mod = (mod * 10 + (s[j] - '0')) % prime[i];
			if (mod == 0 && prime[i] < L)
			{
				cout << "BAD " << prime[i] << endl;
				flag = false;		//  改变标记，避免重复输出
				break;
			}
		}
		if (flag)
			cout << "GOOD" << endl;
	}
	return 0;
}