HDU-1005 Number Sequence && 51NOD-1126 求递推序列的第N项

来源：互联网发布：史丹利快报的淘宝店编辑：程序博客网时间：2024/04/24 11:20

有一个序列是这样定义的：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 longint 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;}

1 0