HDU 5451 Best Solver

Problem Description

The so-called best problem solver can easily solve this problem, with his/her childhood sweetheart.

It is known that  y=(5+26√)1+2x .
For a given integer  x (0≤x<232)  and a given prime number  M (M≤46337) , print  [y]%M . ( [y]  means the integer part of  y )

 

Input

An integer  T (1<T≤1000) , indicating there are  T  test cases.
Following are  T  lines, each containing two integers  x  and  M , as introduced above.

 

Output

The output contains exactly  T  lines.
Each line contains an integer representing  [y]%M .

 

Sample Input

  
  

   
   7
0 46337
1 46337
3 46337
1 46337
21 46337
321 46337
4321 46337
  
  

 

Sample Output

  
  

   
   Case #1: 97
Case #2: 969
Case #3: 16537
Case #4: 969
Case #5: 40453
Case #6: 10211
Case #7: 17947
  
  
  
  


  
  
  
  

   
   矩阵乘法+循环节判断
  
  
  
  

   
   
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<iostream>
using namespace std;
typedef long long LL;
int T,t=0,m,f[50000];
unsigned int n;

struct martix
{
    unsigned int a[2][2];
    martix(){memset(a,0,sizeof(a));}
    martix operator*(const martix &b)
    {
        martix c;
        for (int i=0;i<2;i++)
            for (int j=0;j<2;j++)
                for (int k=0;k<2;k++)
                    (c.a[i][k]+=a[i][j]*b.a[j][k])%=m;
        return c;
    }
};

bool operator!=(const martix &a,const martix &b)
{
    if (a.a[0][0]!=b.a[0][0]) return true;
    
    if (a.a[0][1]!=b.a[0][1]) return true;
    
    if (a.a[1][0]!=b.a[1][0]) return true;
    return false;
}

int get(unsigned int n,int m)
{
    int ans=1;
    for (int i=2;n;n>>=1)
    {
        if (n&1) (ans*=i)%=m;
        i=(i*i)%m;
    }
    return (ans+1)%m;
}

martix mget(martix a,int x)
{
    martix c=a;
    if (x<=0) return c;
    for (x-=1;x;x>>=1)
    {
        if (x&1) c=c*a;
        a=a*a;
    }
    return c;
}

int query()
{
    if (f[m]) return f[m];
        martix b,c;
        int loop=1;
        b.a[0][0]=10%m;
        b.a[0][1]=1%m;
        b.a[1][0]=(m-1)%m;
        c=b*b;
        while (b!=c) {c=c*b; loop++; }
        f[m]=loop;
        return loop;
}

int main()
{
    memset(f,0,sizeof(f));
    scanf("%d",&T);
    while (T--)
    {
        scanf("%u%d",&n,&m);
        martix a,b,c;
        int loop=1;
        a.a[0][0]=98%m;
        a.a[0][1]=10%m;
        b.a[0][0]=10%m;
        b.a[0][1]=1%m;
        b.a[1][0]=(m-1)%m;
        n=get(n,query());
        c=mget(b,n-1);
        a=a*c; 
        printf("Case #%d: %d\n",++t,(a.a[0][1]+m-1)%m);
    }
    return 0;
}