题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6827
解题思路:本题多次使用逆元操作,使用快速幂时间复杂度为O(nlogn),使用线性求逆元时间复杂度为O(n)
线性求逆元
求
可以按规模k来寻找规律
当k=1时,分子为ans[1]=(a1+a2+…+an)/1=(sum[n+k-1]-sum[k-1])/1=p[1]/1;
k=2时,分子为ans[2]=[(a1+a2)+(a2+a3)+(a3+a4)+…(an-1 +an)]/2=[(p[1]+(sum[n+k-1]-sum[k-1])]/2=p[2]/2;
…
可以发现每次规模扩大的时候,都会多加上sum[n+k-1]-sum[k-1]
每一步进行取逆元操作,一共有n*(n-1)/2种情况,最后再对n*(n-1)/2取逆元即可
#include<iostream>
#include<cstdio>
#include<string.h>
#include<vector>
#include<string>
using namespace std;
#define ll long long
const int mod=1e9+7;
const int maxn=2e5+100;
int t,n;
int s[maxn];
ll inv[maxn];
ll ans;
ll sum[maxn];
ll p[maxn];
ll fastpow(ll base, ll n, ll mod) {
ll ans = 1;
while (n) {
if (n & 1) ans *= base % mod, ans %= mod;
base *= base, base %= mod;
n >>= 1;
}
return ans % mod;
}
int main(){
inv[1]=1;
for(int i=2;i<maxn;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
cin>>t;
while(t--){
ans=0;
cin>>n;
for(int i=1;i<=n;i++){
scanf("%d",&s[i]);
sum[i]=(sum[i-1]+s[i])%mod;
}
for(int i=1;i<=n;i++){
p[i]=(p[i-1]+(sum[n+1-i]-sum[i-1]+mod)%mod)%mod;
ans=(ans+p[i]*inv[i]%mod)%mod;
}
ll div=fastpow(1ll*(1+n)*n/2%mod,mod-2,mod);
printf("%lld\n",ans*div%mod);
}
return 0;
}