poj2635

来源：互联网发布：caxa制造工程师软件编辑：程序博客网时间：2024/05/16 18:54

这个题主要是要学会高精度求模。

思路有二，首先转化成千进制，然后利用“同余模定理”，或称“同余定理”。

比较实用的几条规律：

例如，想求123%4，可以这么操作

(100%4+20%4+3%4)%4。利用了性质1、性质3。

为何把它化作千进制，因为利用性质三，可知其是正确的。化作千进制，求模就变成了不超过1000的数的模，虽然次数增多，但是可以求解大数的模。

#include <iostream>#include <string>#include <math.h>using namespace std;#define Max 1000000 + 100bool prime[Max];int p[100000];int IsPrime(){prime[0] = prime[1] = 0;prime[2] = 1;for (int i = 3; i < Max; i++){prime[i] = i % 2 == 0 ? 0 : 1;}int t = (int) sqrt(Max*1.0);for (int i = 3; i <= t; i++){if (prime[i]){for (int j = i * 2; j < Max; j += i){prime[j] = 0;}}}int j = 1;for (int k = 0; k < Max; k++){if (prime[k]){p[j++] = k;}}return j;}bool mod(const int* K, const int p, const int len){int sq = 0;for (int i = len - 1; i >= 0; i--)sq = (sq * 1000 + K[i])%p; if (!sq)   //K被整除  return false;return true;}int main(){int num_of_prime = IsPrime();int num[100];string s;int l;while (cin >> s >> l && l){memset(num, 0, sizeof(num));for (int i = 0; i < s.size(); i++)//局部顺序，全局倒序  {int pKt = (s.size() + 2 - i) / 3 - 1;num[pKt] = num[pKt] * 10 + (s[i] - '0');}int length_of_num = (s.size() + 2) / 3;bool flag = true;for (int i = 1; p[i] < l; ++i){if (!mod(num, p[i], length_of_num)){flag = false;cout << "BAD " << p[i] << endl;break;}}if (flag)cout << "GOOD" << endl;}}