HDU 2817 A sequence of numbers

题目链接http://acm.hdu.edu.cn/showproblem.php?pid=2817

快速幂取模问题

套用模板就可以

View Code

 1 #include<stdio.h>
 2 #include<string.h>
 3 #define M 200907
 4 __int64 quickpow(__int64 m,__int64 n)
 5 {
 6     int b=1;
 7     while(n>0)
 8     {
 9         if(n&1)
10         b=(b*m)%M;
11         n=n>>1;//位运算 相当于n^2
12         m=(m*m)%M;
13     }
14     return b;
15 }
16 int main()
17 {
18     __int64 n,a,b,c,d,s,m,q;
19     scanf("%I64d",&n);
20     while(n--)
21     {
22         scanf("%I64d %I64d %I64d",&a,&b,&c);
23         scanf("%I64d",&d);
24         if((b-a)==(c-b))
25         {
26             q=(a+(b-a)*(d-1))%M;
27             printf("%I64d\n",q);
28         }
29         else
30         {
31                 if(d>1)
32                 {
33                     s=b/a;
34                     m=a*quickpow(s,d-1);
35                 printf("%I64d\n",m%M);
36                 }
37                 else
38                 printf("%I64d\n",a);
39         }
40     }
41     return 0;
42 }

模板 a^p mod m

int Montgomery(int a, int p, int m)
{
    if(p==0) return 1;
    int r=a%m;
    int k=1;
    while(p>1){
        if(p&1!=0){
            k=(k*r)%m;
        }
        r=(r*r)%m;
        p/=2;
    }
    return (r*k)%m;
}