本问题可以分成三部分求解:
#include
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
ll n = 1001733993063167141, d = 212353, C = 20190324;
ll extgcd(ll a, ll b, ll &x, ll &y)
{
if (b == 0)
{
x = 1, y = 0;
return a;
}
ll ret = extgcd(b, a % b, y, x);
y -= (a / b) * x;
return ret;
}
ll quickmi(ll x, ll n, ll m)
{
ll ans = 1;
while (n > 0)
{
if (n & 1)
ans = (__int128)ans * x % m;
x = (__int128)x * x % m;
n >>= 1;
}
return ans;
}
int main()
{
// 求p, q
ll p, q;
for (p = 2; p * p <= n; p++)
if (n % p == 0)
{
q = n / p;
break;
}
cout << "p = " << p << " "
<< "q = " << q << endl;
// 求e
ll MOD = (p - 1) * (q - 1), x, y;
extgcd(d, MOD, x, y);
ll e = (x % MOD + MOD) % MOD;
cout << "e = " << e << endl;
// 求原文X
cout << quickmi(C, e, n) << endl;
return 0;
}
注:
long long
会溢出,所以要采用更大的类型,如 __int128
或 __int256
。完