90. 64位整数乘法

原题链接

根据快速幂的思想

b=(b1*1+b2*2+b3*4...+bn*2^n-1)

res=a*(b1*1+b2*2+...)

#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;

int main()
{
    ll a,b,p;
    cin>>a>>b>>p;
    ll res=0;

    while(b)
    {
        if(b&1)
        {
            res=(res+a)%p;
        }
        b>>=1;
        a=a*2%p;
    }
    cout<<res%p;
    return 0;
}

应用模的定义出发

由于 a%b=a-(a/p)*p   （a/b是下取整)

那么 a*b%p=a*b-(a*b/p)*p

a*b和(a*b/p)*p 的差一定是在 0~p-1 之间的

如果res=a*b%p <0 res+=p;        

        res>=p res-=p;
 

#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
int main()
{
    ll a,b,p;
    cin>>a>>b>>p;
    ll c=(long double)a*b/p;
    ll res=a*b-c*p;
    if(res<0) res+=p;
    else if (res>=p) res-=p;
    cout<<res;
    return 0;
    
}