小测验:震惊——题解

这是一道快速幂的题……
题面大(简)致(化)如下:
给你很多个分数,求他们的和对p取模的值。
其中(a/b)%p=a*(b在模p的意义下的乘法逆元)
p=1e9+7

恩我们一看乘法逆元,我们想到了什么哇?
费马小定理!
假如p是质数,且gcd(a,p)=1,那么 a^(p-1)≡1(mod p)
首先我们知道,1e9+7是质数,且a%p小于p,所以gcd(a,p)=1
那么a^(p-1)≡1(mod p)导出a^(p-1)%p=1
再导出a^(p-2)%p*a%p=1
其次我们设a在模p的意义下的乘法逆元为x,则有a%p=1/x
那么移项可得到x=a^(p-2)%p
好的我们已经会了这道题。

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<iostream>
using namespace std;
const long long p=1e9+7;
long long qpow(long long k,long long n){
    if(n==0)return 1;
    if(n==1)return k%p;
    int q=qpow(k,n/2)%p;
    if(n%2==0)return q%p*q%p;
    else return q%p*q%p*k%p;
}
long long sum=0;
int main(){
    int n,a,b;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d %d",&a,&b);
        sum+=a*qpow(b,1e9+5)%p;
        sum%=p;
    }
    printf("%lld",sum);
    return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章