0
点赞
收藏
分享

微信扫一扫

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


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


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;
}



举报

相关推荐

0 条评论