传送门
文章目录
- 题意:
- 思路:
题意:
给你两个长度为nnn的数组a,ba,ba,b,你可以至多反转一段连续区间,求∑i=1nai∗bi\sum _{i=1}^n a_i*b_i∑i=1nai∗bi最大是多少。
n<=5e3n<=5e3n<=5e3
思路:
首先我们可以通过n2n^2n2来枚举所有的区间,但是要计算翻转之后的贡献的话还需要多加一个nnn,这样复杂度是n3n^3n3,显然不可接受。
考虑能否通过小区间来递推出大区间。
定义f[i][j]f[i][j]f[i][j]为将[i,j][i,j][i,j]段区间翻转之后变换的值,我们发现如果我们知道一个小区间[l,r][l,r][l,r],那么[l−1,r+1][l-1,r+1][l−1,r+1]的区间贡献就是在[l,r][l,r][l,r]的贡献基础上再将l−1,r+1l-1,r+1l−1,r+1两个位置的值互换了一下,互换两个值l,rl,rl,r对答案的贡献为a[l]∗b[l]+a[r]∗b[r]+a[l]∗b[r]+a[r]∗b[l]=a[l]∗(b[r]−b[l])+a[r]∗(b[l]−b[r])a[l]*b[l]+a[r]*b[r] +a[l]*b[r]+a[r]*b[l]=a[l]*(b[r]-b[l])+a[r]*(b[l]-b[r])a[l]∗b[l]+a[r]∗b[r]+a[l]∗b[r]+a[r]∗b[l]=a[l]∗(b[r]−b[l])+a[r]∗(b[l]−b[r]),所以我们可以通过[l,r][l,r][l,r]扩展到[l−1,r+1][l-1,r+1][l−1,r+1],所以转移方程为f[l][r]=f[l−1][r+1]+(a[l]−a[r])∗(b[r]−b[l])f[l][r]=f[l-1][r+1]+(a[l]-a[r])*(b[r]-b[l])f[l][r]=f[l−1][r+1]+(a[l]−a[r])∗(b[r]−b[l])
复杂度O(N2)O(N^2)O(N2)
#include<bits/stdc++.h>
using namespace std;const int N=5010;typedef long long LL;int n;
LL a[N],b[N];
LL f[N][N],pre[N];int main() {cin>>n;LL ans=0;for(int i=1;i<=n;i++) cin>>a[i];for(int i=1;i<=n;i++) cin>>b[i];for(int i=1;i<=n;i++) ans+=a[i]*b[i],pre[i]=ans;for(int len=2;len<=n;len++) {for(int l=1;l+len-1<=n;l++) {int r=l+len-1;f[l][r]=f[l+1][r-1]+(a[l]-a[r])*(b[r]-b[l]);ans=max(ans,f[l][r]+pre[n]);}}printf("%lld\n",ans);return 0;
}//a[l]*b[l]+a[r]*b[r] a[l]*b[r]+a[r]*b[l]
//a[l]*(b[r]-b[l])+a[r]*(b[l]-b[r])