来源:互联网 发布:sql delete一条数据 编辑:程序博客网 时间:2024/04/28 14:58

不是给清北打广告!!!!!!
这里写图片描述
其实就是个等比数列求和,只是除不满足同余的性质,需要求一下逆元
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; } 
原创粉丝点击