hdu3977 - fibonacci模p的周期

假设N(m)为模m的周期，则有

N(ab)=lcm(N(a),N(b))，其中gcd(a,b)=1; /*这里证明使用了fibo数的性质，

N(p^k)=N(p)*p^(k-1)，其中p为质数   /*论文里这个没有证明，但是验证了10^14内的都满足

对于质数p,其周期可以推算得出：

当p%5=1或4时，它是p-1的因子

当p%4=2或3时，它是2*(p-1)的因子

如果便只用验证因子是否是循环节就行。

求fibo数，可以使用矩阵乘法，

#include<stdio.h>
#include<string.h>
#include<math.h>
#define maxp 1000000
#define LL long long

LL P;
int pri[maxp+10];
int b[maxp+10];
int tot;

struct data {
    int p,t;
    LL r;
} q[500];
int qt;

LL fab[15];

void getprime()
{
    memset(b,0,sizeof(b));
    tot = 0;
    int i=2;
    while (i<=maxp) {
        while (b[i] && i<=maxp)
            i++;
        pri[++tot]=i;
        int j=i;
        while (j<=maxp) {