【C++题解】【NOIP2005PJ】循环(circle.pas/cpp)

题目链接：https://gmoj.net/junior/#main/show/1130

题目描述

乐乐是一个聪明而又勤奋好学的孩子。他总喜欢探求事物的规律。一天，他突然对数的正整数次幂产生了兴趣。
众所周知，2的正整数次幂最后一位数总是不断的在重复2，4，8，6，2，4，8，6……我们说2的正整数次幂最后一位的循环长度是4（实际上4的倍数都可以说是循环长度，但我们只考虑最小的循环长度）。类似的，其余的数字的正整数次幂最后一位数也有类似的循环现象：
这时乐乐的问题就出来了：是不是只有最后一位才有这样的循环呢？对于一个整数n的正整数次幂来说，它的后k位是否会发生循环？如果循环的话，循环长度是多少呢？
注意：
1． 如果n的某个正整数次幂的位数不足k，那么不足的高位看做是0。
2． 如果循环长度是L，那么说明对于任意的正整数a，n的a次幂和a + L次幂的最后k位都相同。

输入

输入文件circle.in只有一行，包含两个整数n（1 <= n < 10^100）和k（1 <= k <= 100），n和k之间用一个空格隔开，表示要求n的正整数次幂的最后k位的循环长度。

输出

输出文件circle.out包括一行，这一行只包含一个整数，表示循环长度。如果循环不存在，输出-1。

样例输入

32 2

样例输出

4

数据范围限制

【数据规模】 对于30%的数据，k <= 4；
对于全部的数据，k <= 100。

做法

一道模拟题，需要用到高精度。

首先，既然我们要求的是后 $k$ 位的循环节，那前 $k$ 位（不包括第 $k$ 位）的循环节我们就不管了。

接着，如果我们已知后 $k-1$ 位的循环节，那么后 $k$ 位的循环节一定就是它的倍数。那我们可以一位一位地找，这里以样例为 $8123$ $4$ 为例：

8123 1
第 1 次：8123 * 8123 = 65983129
第 2 次：3129 * 8123 = 25416867
第 3 次：6867 * 8123 = 55780641
第 4 次：0641 * 8123 = 05206843
8123 ^ 4 = 4353773312630641

我们可以发现，在乘第 $4$ 次的时候，最后 $1$ 位出现了循环，循环节为 $4$ ，那么最后的答案一定就是 $4$ 的倍数。

那我们可以把乘数换一下，换成 $0641$ ，即 $8123^4$ 的后 $4$ 位。这样乘 $1$ 次 $0641$ 就等于乘了 $4$ 次 $8123$ ，而且，乘 $0641$ 还能保证最后 $1$ 位永远是 $3$ ，因为 $8123^4$ 的最后一位都是 $3$ 。

这样，我们就可以不用再去管后 $1$ 位了，而去处理后 $2$ 位。

8123 2
第 1 次：8123 * 0641 = 05206843
第 2 次：6843 * 0641 = 04386363
第 3 次：6363 * 0641 = 04078683
第 4 次：8683 * 0641 = 05565803
第 5 次：5803 * 0641 = 03719723
0641 ^ 5 = 108215668739201

最后 $2$ 位的循环节是 $5$，我们再把乘数换成 $9201$ 。

8123 3
第 1 次：8123 * 9201 = 74739723
第 2 次：9723 * 9201 = 89461323
第 3 次：1323 * 9201 = 12172923
第 4 次：2923 * 9201 = 26894523
第 5 次：4523 * 9201 = 41616123
9201 ^ 5 = 65943979755726446001

8123 4
第 1 次：8123 * 6001 = 48746123
第 2 次：6123 * 6001 = 36744123
第 3 次：4123 * 6001 = 24742123 
第 4 次：2123 * 6001 = 12740123
第 5 次：0123 * 6001 = 00738123

我们继续处理完后 $3$ 位和后 $4$ 位，得出了每 $1$ 位的循环节，答案只要把它们累乘起来就可以了，即最后的答案为 $4\times 5\times 5\times 5=500$ 。

至此，我们求出了 $8123$ 后 $4$ 位的循环长度。

由于我们在处理的时候，已经确保了后面的数字都不变，那么我们只要考虑当前这 $1$ 位的循环节就行了，而 $1$ 位的循环节最多为 $10$ 。所以，我们处理每 $1$ 位的时候，就都只要处理 $10$ 次，如果还没出结果，那就直接输出 $-1$ 并退出就行了。

因为我们每 $1$ 位最多处理 $10$ 次，每一次最多要耗费 $k^2$ 的时间，所以时间复杂度约 $O(k^3)$，妥妥地能过！

还要注意一点，最后的答案要用高精度，否则会因为超过 long long 的范围而答案错误。

#include<cstdio>
#include<cstdlib>
#include<cstring>
int a[205],b[205],s[205],d[205],f[205],ans[205];
int k,clen,cnt,t,ss=1;
char c[205];
void lc(int x)//累乘 
{
	memset(s,0,sizeof(s));
	//高精乘高精 
	for(int i=1;i<=x;i++)
	for(int j=1;j<=x;j++) s[i+j-1]+=a[i]*b[j];
	for(int i=1;i<=x;i++)
	{
		s[i+1]+=s[i]/10;
		s[i]%=10;
	}
	for(int i=1;i<=x;i++) a[i]=s[i];//累乘，把得到的积变成下一个的乘数 
	return;
}
bool bj(int x)//比较 
{
	for(int i=1;i<=x;i++)
	if(a[i]!=d[i]) return 0;//是否相同 
	return 1;
}
int cl(int x)//处理 
{
	cnt=0;
	for(int i=1;i<=x;i++) a[i]=d[i],b[i]=f[i];//重置乘数 
	do
	{
		lc(x);//累乘 
		cnt++;
		if(cnt>10)//处理超过10次就退出 
		{
			printf("-1\n");
			exit(0);
		}
	}while(!bj(x));
	for(int i=1;i<=k;i++) a[i]=0,b[i]=f[i];//重置乘数 
	a[1]=1;
	for(int i=1;i<=cnt;i++) lc(k);//累乘，得到新的乘数 
	for(int i=1;i<=k;i++) f[i]=a[i];//替换乘数 
	return cnt;
}
void cf(int x)//乘法 
{
	int w=0;
	for(int i=1;i<=ss;i++)//高精乘低精 
	{
		ans[i]=ans[i]*x+w;
		w=ans[i]/10;
		ans[i]%=10;
		if(i==ss&&w) ss++;
	}
	return;
}
int main()
{
	freopen("circle.in","r",stdin);
	freopen("circle.out","w",stdout);
	scanf("%s%d",c+1,&k);
	clen=strlen(c+1);
	for(int i=1;i<=k;i++) d[i]=f[i]=c[clen-i+1]-'0';//取后k位 
	ans[1]=1;
	for(int i=1;i<=k;i++)
	{
		t=cl(i);//获取每一位的循环节 
		cf(t);//累乘每一位的循环节 
	}
	for(int i=ss;i>=1;i--) printf("%d",ans[i]);
	return 0;
}

写文章不易，望大家给 JG_DF_ 点个赞！