有一个序列是这样定义的:f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7.
给出A,B和N,求f(n)的值。
前者:Each test case contains 3 integers A, B and n on a single line (1 <= A, B <= 1000, 1 <= n <= 100,000,000).
后者:输入3个数:A,B,N。数字之间用空格分割。(-10000 <= A, B <= 10000, 1 <= N <= 10^9)
被后者坑了将近一个小时,mod 7跟计算机的 % 7 不一样,mod只能是非负数,而%则会有负数,这个坑,够深的啊……
要点:
1、找循环部分(即重复子序列),不一定是1,1开头。
2、(….%7+7)%7
51NOD-AC代码:
#include <iostream>
#include <cstring>
using namespace std;
#define LL long long
int f[1234];
int main()
{
f[1]=1;f[2]=1;
LL A,B,n;
int Flag[8][8];
int i;
while(cin>>A>>B>>n)
{
int flag2=0;
memset(Flag,0,sizeof(Flag));
if(n<3)
{
cout<<1<<endl;
continue;
}
for(i=3;i<=400;++i)
{
f[i]=((A*f[i-1]+B*f[i-2])%7+7)%7;
if(f[i]==1&&f[i-1]==1)break;
if(Flag[f[i]][f[i-1]])
{
flag2=Flag[f[i]][f[i-1]];
break;
}
Flag[f[i]][f[i-1]]=i;
}
if(i>n)
{
cout<<f[n]<<endl;
continue;
}
if(flag2)
{
cout<<f[(n-flag2)%(i-flag2)+flag2]<<endl;
continue;
}
i-=2;
n%=i;
if(!n)n=i;
cout<<f[n]<<endl;
}
return 0;
}