洛谷 P1050 循环

题目描述
乐乐是一个聪明而又勤奋好学的孩子。他总喜欢探求事物的规律。一天，他突然对数的正整数次幂产生了兴趣。
众所周知，2的正整数次幂最后一位数总是不断的在重复2，4，8，6，2，4，8，6……我们说2的正整数次幂最后一位的循环长度是4（实际上4的倍数都可以说是循环长度，但我们只考虑最小的循环长度）。类似的，其余的数字的正整数次幂最后一位数也有类似的循环现象：
循环 循环长度
2 2、4、8、6
4
3 3、9、7、1
4
4 4、6 2
5 5 1
6 6 1
7 7、9、3、1
4
8 8、4、2、6
4
9 9、1 2
这时乐乐的问题就出来了：是不是只有最后一位才有这样的循环呢？对于一个整数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。

输入输出样例
输入样例#1：
32 2
输出样例#1：
4
说明
对于30%的数据，k <= 4；
对于全部的数据，k <= 100。
NOIP2005普及组第四题

---

这题，出成普及组的真是屈才了。
首先，听说暴力模拟30分（反正我没试过）。
说正解吧，如果第k位的循环节长度是l，那么第k+1位循环节的长度一定是l的整数倍。
设：t等于这个数的后k位（因为只有后k位有用）。
d等于这个数乘d后k位相同。s表示上一次这个数的循环节次方的后k位。
每次s*d，t*s，就可以递推了。
判断-1的方法：
开一个bool数组记录这一位出现过那些数，第一次重复的是原数字，就找到了循环节，否则，就是-1.
注意：后7个点循环节也需要高精度。
（调试到再也不想写高精度）。
感谢Orion_Rigel神犇指导。

---

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
int ori[105],s[105],t[105],c[105],d[105],k,loop[1005],len;//t是积，s是乘数 
bool b[15];
char ch[105];
void mul1()
{
    memset(c,0,sizeof(c));
    for(int i=1;i<=k;i++)
        for(int j=1;j<=k;j++)
            if(i+j-1<=k)
                c[i+j-1]+=d[i]*s[j];
    for(int i=1;i<=k;i++)
        c[i+1]+=c[i]/10,c[i]%=10;
    memcpy(d+1,c+1,k*sizeof(int));
}
void mul2()
{
    memset(c,0,sizeof(c));
    for(int i=1;i<=k;i++)
        for(int j=1;j<=k;j++)
            if(i+j-1<=k)
                c[i+j-1]+=s[i]*t[j];
    for(int i=1;i<=k;i++)
        c[i+1]+=c[i]/10,c[i]%=10;
    memcpy(t+1,c+1,k*sizeof(int));
}
void mul(int x)
{
    for(int i=1;i<=loop[0];i++)
        loop[i]*=x;
    for(int i=1;i<=loop[0];i++)
        loop[i+1]+=loop[i]/10,loop[i]%=10;
    if(loop[loop[0]+1]>0)
        ++loop[0];
    while(loop[loop[0]]/10!=0)
    {
        loop[loop[0]+1]=loop[loop[0]]/10,loop[0]%=10;
        ++loop[0];
    }
}
bool dfs(int dep)
{
    if(dep>k)
        return 1;
    int y=ori[dep];
    memset(b,0,sizeof(b));
    b[ori[dep]]=1;
    memset(d,0,sizeof(d));
    d[0]=k,d[1]=1;
    for(int i=1;;i++)
    {
        mul1();//d*s
        mul2();//t*s
        if(b[t[dep]])
        {
            if(t[dep]==y)
            {
                mul(i);
                memcpy(s+1,d+1,k*sizeof(int));
                for(int i=0;i<=len;i++)
                    t[i]=ori[i];
                return dfs(dep+1);
            }
            else
                return 0;
        }
        else
            b[t[dep]]=1;
    }
}
int main()
{
    scanf("%s%d",ch+1,&k);
    len=strlen(ch+1);
    ori[0]=k;
    loop[0]=loop[1]=1;
    for(int i=1;i<=len;i++)
        ori[len-i+1]=ch[i]-'0';
    for(int i=0;i<=len;i++)
        s[i]=ori[i];
    for(int i=0;i<=len;i++)
        t[i]=ori[i];
    if(dfs(1))
        for(int i=loop[0];i>=1;i--)
            printf("%d",loop[i]);
    else
        printf("-1");
    printf("\n");
    return 0;
}