题意: 给出K1,求一个12位数(不含前导0)K2,使得K1^K2 mod (10^12) = K2.
解法: 求不动点问题。
有一个性质: 如果12位数K2满足如上式子的话,那么K2%1,K2%10,K2%100,…,K2%10^12都会满足如上式子。那么我们可以dfs从后往前一个一个找出这个数的每一位。
转载至:http://www.bubuko.com/infodetail-587732.html
这个不知名的大佬不知道是谁,弱弱把这个代码解读了一遍。
有这么一个问题 在这个dfs里面 跑了13层。也就是最后得到k2是一个13位的,为什么不是在12位结束呢,明明只需要求一个12位的数字啊,为什么要这样做呢?
我觉得是因为在代码里到12位的时候有可能是0开头的情况 这样会得到一个只有11位的数,然而很糟糕的是,如果只填完12位就结束,这样有可能会让k1^k2==k2 这就不符合规定的12位了。所以我们把k2遍历到第13位。那么无论如何结束的也得在12位结束。因为最后判断的是 tmp>=12位的就行就行,就是11个0。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#define SMod 1000000000000
#define ll long long
using namespace std;
#define N 10007
ll K1,K2;
inline long long mul(long long x,long long y,long long mod)
{
long long tmp=(x*y-(long long)((long double)x/mod*y+1.0e-8)*mod);
return tmp<0 ? tmp+mod : tmp;
}
ll fastm(ll a,ll b,ll m) {
ll res = 1LL;
while(b) {
if(b&1LL) res = mul(res,a,m);
a = mul(a,a,m);
b >>= 1;
}
return res;
}
ll wei[16],ans;
bool dfs(int c,ll now) {
if(c == 13) {
if(now >= wei[11]) { ans = now; return true; }
return false;
}
ll W = wei[c];
for(ll i=0;i<=9;i++) {
ll tmp = W*i+now;
if(fastm(K1,tmp,W) != tmp%W) continue;
// if(c==12)
// printf("%lld\n",tmp );
if(dfs(c+1,tmp)) return true;
}
return false;
}
int main()
{
int t,cs = 1;
wei[0] = 1LL;
for(int i=1;i<=13;i++) wei[i] = wei[i-1]*10LL;
//预处理1代表1e 12代表12e 也就是13位 12个0
while(scanf("%lld",&K1)!=EOF && K1) {
dfs(0,0);
printf("Case %d: Public Key = %lld Private Key = %lld\n",cs++,K1,ans%SMod);
}
return 0;
}