不是给清北打广告!!!!!!
其实就是个等比数列求和,只是除不满足同余的性质,需要求一下逆元
a/bmodp,在这个题中,b与p是互质的,可以用费马小定理去求逆元(还可以用拓展欧几里得,因为写了快速幂了,小费马可以直接搞,懒得弄了)
a/bmodp=a*b^p-2%p
这样就可以取模了
用上这个求和公式的方法大约的时间复杂度为O(n*logm),比暴力强多了
当然还需要注意一下,当公比是1时,要特判一下,因为分母不能为0
#include <cstdio>
#include <iostream>
using namespace std;
long long sum;
int n,m;
const int p=1e9+7;
const int mod=1e9+7;
inline void Calculation(long long x)
{
long long b=x;
for(int j=1;j<=m;j++)
{
sum=(sum+b)%p;
b=(b*x)%p;
}
}
inline long long power_mod(long long a,long long b)
{
long long ans=1;
while (b)
{
if (b%2) ans=ans*a%mod;
a=a*a%mod;
b/=2;
}
return ans;
}
int main()
{
scanf("%d%d",&n,&m);
sum=(m%mod);
for(int i=2;i<=n;i++)
{
sum=(sum+((i*(1-power_mod(i,m))%mod)%mod*(power_mod(1-i,p-2))%mod)%mod)%mod;
}
printf("%lld",sum);
return 0;
}