观察可知∑nl=1∑nr=lf(l,r)∑l=1n∑r=lnf(l,r)的意义为长度为r−l+1r−l+1的所有区间的和与wr−l+1wr−l+1的乘积。
枚举ll、势必会超时,于是我们找一下规律。
以n=7n=7为例。
区间长度为1时,如下图所示:
每一条粗线表示每个元素被加的次数,也就是说,每个元素都被加了1次。
区间长度为2时,如下图所示:
第1、7个元素被加了1次,其余元素被加了2次。
区间长度为3时,如下图所示:
第1、7个元素被加了1次,第2、6个元素被加了2次,第3、4、5个元素被加了3次。
区间长度为4时,如下图所示:
区间长度为5时,如下图所示:
区间长度为6时,如下图所示:
区间长度为7时,如下图所示:
可以发现,各数组元素加的次数先增后减,且是对称的。各区间长度的和也是对称的。用前缀和处理即可做到O(n)O(n)求解。
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
typedef long long ll;
const int MOD=1e9+7;
int n,pnt;ll a[310101],w[310101],suma[310101],sumb[310101];
ll sumc[310101];
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%lld",&a[i]);
suma[i]=suma[i-1]+a[i];
suma[i]%=MOD;
}
for(int i=1;i<=n;i++)scanf("%lld",&w[i]);
sumb[1]=suma[n];
for(int i=2;i<=n;i++)
{
if(n-i+1>i-1)
{
sumb[i]=sumb[i-1]+(suma[n-i+1]-suma[i-1]);
sumb[i]%=MOD;
}
else
{
pnt=i-1;
break;
}
}
sumb[n]=suma[n];
for(int i=n-1;i>=pnt+1;i--)
{
sumb[i]=sumb[n-i+1];
}
ll ans=0;
for(int i=1;i<=n;i++)
{
ans+=(sumb[i]*w[i])%MOD;
ans%=MOD;
}
ans%=MOD;
printf("%lld\n",ans);
return 0;
}