BBS.zip_bbs_密码学

// BBS.cpp : 此文件包含 "main" 函数。程序执行将在此处开始并结束。 // #include <iostream> #include <cstdio> #include<math.h> using namespace std; typedef long long LL; //****************************素数检验算法begin********************************* LL mulmod(LL a, LL b, LL p) { LL d = 1; a = a % p; while (b > 0) { if (b & 1) d = (d * a) % p; a = (a * a) % p; b >>= 1; } return d; } bool witness(LL a, LL n) { LL d = n - 1; if (n == 2) return true; if (!(n & 1)) return false; while (!(d & 1)) d = d / 2; LL t = mulmod(a, d, n); while ((d != n - 1) && (t != 1) && (t != n - 1)) { t = mulmod(t, 2, n); d = d << 1; } return (t == n - 1) || (d & 1); } bool isprime(LL n)//米勒拉宾算法——素性测试 { int a[3] = { 2,7,61 }; for (int i = 0; i < 3; i++) if (!witness(a[i], n)) return false; return true; } //*******************************素数检验算法end******************************* int main() { LL p = 50000, q ; for (; p <= 65536; p++) { if (isprime(p) && p % 4 == 3) { cout << "p=" << p; break; } } cout << endl; for (q=p+1; q <= 65536; q++) { if (isprime(q) && q % 4 == 3) { cout << "q=" << q; break; } } cout << endl; LL n = p * q; cout << "n=" << n << endl; LL s = 5; cout << "s=" << s << endl; LL Xi = (s * s) % n; cout << "X0=" << Xi << endl; for (int i = 1; i <= 100; i++)//这里只输出100个数作为测试 { Xi = (Xi * Xi) % n; if (Xi >= 0) cout << Xi % 2; else cout << "1"; } } // 运行程序: Ctrl + F5 或调试 >“开始执行(不调试)”菜单 // 调试程序: F5 或调试 >“开始调试”菜单 // 入门提示: // 1. 使用解决方案资源管理器窗口添加/管理文件 // 2. 使用团队资源管理器窗口连接到源代码管理 // 3. 使用输出窗口查看生成输出和其他消息 // 4. 使用错误列表窗口查看错误 // 5. 转到“项目”>“添加新项”以创建新的代码文件，或转到“项目”>“添加现有项”以将现有代码文件添加到项目 // 6. 将来，若要再次打开此项目，请转到“文件”>“打开”>“项目”并选择 .sln 文件