序列求和（费马小定理除法取余）

题目链接：https://ac.nowcoder.com/acm/problem/15950

序列求和（费马小定理除法取余）

序：存在(a*b)%p=((a%p)*(b%p))%p
但是不存在 (a/b)%p=(a%p)/(b%p)%p

费马小定理

如果p是一个质数，而整数a不是p的倍数，则有a^p-1≡1（mod p）。

我们要求(a/b)%c
即(a*b^-1*1)%c
因为b^c-1%c=1;带入上式
得(a/b)%c = (a*b^c-2)

一般c都是一个比较大得质数，可以用快速幂来求b^c-2

题目代码

#include<iostream>
#include<cstdio>
#include<cmath>
#include<vector>
#include<string>
#include<algorithm>
#include<queue>
#include<map>
#define ll long long
#define int128 __int128
#define s_1(_1) scanf("%d",&_1)
#define s_2(_1,_2) scanf("%d%d",&_1,&_2)
const int N=1e5+10;
ll Mod=1000000007;
using namespace std;


ll kpow(ll a,ll p)
{
    ll ans=1;
    while(p > 0)
    {
        if(p&1)
            ans=(ans%Mod)*(a%Mod);
        p>>=1;
        a=(a%Mod)*(a%Mod)%Mod;
    }
    return ans%Mod;
}
int main()
{
    unsigned ll n;
    while(cin >> n)
    {
        ll b=kpow(6,Mod-2);
        
        cout << (((n%Mod)*((n+1)%Mod))%Mod)*((((2*n%Mod+1)%Mod*b)%Mod)%Mod)%Mod << endl;
    }
    return 0;
}