问题 C: NOIP2005普及组第4题 循环
时间限制: 1 Sec 内存限制: 128 MB提交: 34 解决: 7
[ 提交][ 状态][ 讨论版][命题人: 外部导入]
题目描述
乐乐是一个聪明而又勤奋好学的孩子。他总喜欢探求事物的规律。一天,他突然对数的正整数次幂产生了兴趣。
众所周知,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位都相同。
输入
只有一行,包含两个整数n(1 <= n < 10100)和k(1 <= k <= 100),n和k之间用一个空格隔开,表示要求n的正整数次幂的最后k位的循环长度。
输出
包括一行,这一行只包含一个整数,表示循环长度。如果循环不存在,输出-1。
【数据规模】
对于30%的数据,k <= 4;
对于全部的数据,k <= 100。
样例输入
32 2
样例输出
4
思路:我们从最后一位数字开始考虑,我们设最后一位数字的循环节位L1,如果要求两位数字的循环节的话,肯定要保证最后一位数字不变,设两位数字的循环节L2,那么L2=mL1,m可能是1到10,如果超出11还不是,说明不存在循环节,以此类推下去,我们就可以求出k位的循环节了,因为结果非常大,需要用大数来表示。
代码:
#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;
const int L=110000;
int l[11]={1,1,4,4,2,1,1,4,4,2};
int k;
void multiplyh(int x[], int y[], int z[]) {//高精度高乘
int up = 0;
for (int ii = 1; ii <= k; ii++) {
for (int j = 1; j <= k; j++)
{
z[ii + j - 1] += (x[j] * y[ii] + up) % 10;
up = (x[j] * y[ii] + up) / 10;
}
up = 0;
}
for (int ii = 1; ii <= k; ii++) {//进位
z[ii + 1] += z[ii] / 10;
z[ii] %= 10;
}
}
void multiplyl(int x[], int yy, int z[]) {//高精度低乘
int up = 0;
for (int ii = 1; ii <= k; ii++) {
z[ii] = (x[ii] * yy + up) % 10;
up = (x[ii] * yy + up) / 10;
}
}
int main()
{
string a;
int b[205],c[205],d[205],now[205],last[205],ss[205];
memset(last, 0, sizeof last);
memset(b, 0, sizeof b); memset(now, 0, sizeof now);
memset(d, 0, sizeof d); memset(c, 0, sizeof c); memset(ss, 0, sizeof ss);
cin>>a;
cin>>k;
int n=a.size();
int z=0;
for(int i=n-1;i>=n-k;i--)
c[++z]=a[i]-'0';
for(int i=1;i<=k;i++)
b[i]=c[i];
for(int i=1;i<l[c[1]];i++)
{
memset(d,0,sizeof(d));
multiplyh(b,c,d);
for(int j=1;j<=k;j++)b[j]=d[j];
}
ss[1]=l[c[1]];
for(int i=1;i<=k;i++)now[i]=b[i];
for(int i=1;i<=k;i++)last[i]=now[i];
int pos=2;
int t=0;
while(pos<=k)
{
for(int i=1;i<=k;i++)
{
b[i]=c[i];
}
t=0;
while(t<11)
{
t++;
memset(d,0,sizeof(d));
multiplyh(b,now,d);
for(int j=1;j<=k;j++)
{
b[j]=d[j];
}
if(b[pos]==c[pos])
{
break;
}
memset(d,0,sizeof(d));
multiplyh(last,now,d);
for(int j=1;j<=k;j++)
{
last[j]=d[j];
}
}
if(t>=11){cout<<-1<<endl;return 0;}
for(int j=1;j<=k;j++)now[j]=last[j];
memset(d,0,sizeof(d));
multiplyl(ss,t,d);
for (int i = 1; i <= 100; i++) ss[i] = d[i];
pos++;
}
int flag = 0;
for (int i = 100; i >= 1; i--) {
if (ss[i]) flag = 1;
if (flag) cout << ss[i];
}
cout<<endl;
return 0;
}