HDU 5666 快速乘

Segment

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 360    Accepted Submission(s): 134

Problem Description

Silen August does not like to talk with others.She like to find some interesting problems.

Today she finds an interesting problem.She finds a segment x+y=q

.The segment intersect the axis and produce a delta.She links some line between (0,0)

and the node on the segment whose coordinate are integers.

Please calculate how many nodes are in the delta and not on the segments,output answer mod P.

Input

First line has a number,T,means testcase number.

Then,each line has two integers q,P.

q

is a prime number,and 2≤q≤1018,1≤P≤1018,1≤T≤10.

Output

Output 1 number to each testcase,answer mod P.

Sample Input

1

2 107

Sample Output

Source

BestCoder Round #80

题意：判断直线x+y=q与坐标轴围成的三角形“内”的整数点的个数 很容易推出公式  （q-1）*（q-2）/2

题解： 直接乘会爆long long 

         快速乘算法 处理

         类比快速幂模拟

http://www.2cto.com/kf/201505/396902.html

ll kuai(ll q,ll num,ll mod){
    ll ans=0;
    ll base=q;
    while(num){
        if(num%2) ans=(ans+base)%mod;
        num/=2;
        base=(base+base)%mod;
    }
    return ans%mod;
}

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #define ll __int64
 using namespace std;
 int t;
 ll q,p;
 ll ans1,ans2;
 ll re;
 ll kuai(ll q,ll num,ll mod){
     ll ans=;
     ll base=q;
     while(num){
         if(num%) ans=(ans+base)%mod;
         num/=;
         base=(base+base)%mod;
     }
     return ans%mod;
 }
 int main()
 {
     while(scanf("%d",&t)!=EOF)
     {
         for(int i=;i<=t;i++)
         {
             scanf("%I64d %I64d",&q,&p);
             if(q%==)
             {
                 ans1=(q-)/;
                 ans1=ans1%p;
                 ans2=q-;
             }
             else
             {
                 ans1=(q-)/;
                 ans1=ans1%p;
                 ans2=q-;
             }
           printf("%I64d\n",kuai(ans1,ans2,p));
         }
     }
     return ;
 }