題意
一條街上有N家燒烤店，從西到東編號爲1至N，第i家和i+1家之間的距離是A[i]。
Joisino有M張餐票，編號從1到M。每家燒烤店都提供M種燒烤套餐，用不同編號的餐票可以換取不同種類的套餐。 在燒烤店i，用編號爲j的餐票可以買到美味值爲B[i][j]的套餐。 每張餐票只能使用一次，但是在每個店可以使用任意數量的餐票。
Joisino希望通過從她選擇的一家店開始，然後反覆前往另一家燒烤店並在該店使用還未使用的餐票來享用總共M個燒烤套餐。
她最終的幸福指數由以下公式計算得出：（所喫套餐的總美味度）-（經過的總路程）。
找到她最大可能的幸福指數。
題解
首先貪心考慮，確定一個區間[L,R]後可以通過ST表O(len)得到答案。考慮對於每個R的最優的L，易得R增大時L不會減小，所以用分治優化，大概是每次暴力找當前區間中點的決策點，然後就可以把這個區間分成兩個小區間，複雜度O(m*n*log(n))。
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
#define N 5005
int n,m;
LL s1[N],ans,f[N];
int a[N][205],lg[N],mx[205][N][20];
inline int read()
{
int x=0,f=1;char c=getchar();
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+(c^48);c=getchar();}
return x*f;
}
inline int Max(int i,int l,int r)
{
int l1=lg[r-l+1];
return max(mx[i][l][l1],mx[i][r-(1<<l1)+1][l1]);
}
inline void dfs(int l,int r,int L,int R)
{
if(l>r)
return;
int mid=(l+r)/2,x1;
LL m1=0;
for(int i=L;i<=min(mid,R);i++)
{
LL ml=s1[i]-s1[mid];
for(int j=1;j<=m;j++)
ml+=Max(j,i,mid);
if(ml>m1)
m1=ml,x1=i;
}
f[mid]=m1;
dfs(l,mid-1,L,x1);
dfs(mid+1,r,x1,R);
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=2;i<=n;i++)
scanf("%lld",&s1[i]),s1[i]+=s1[i-1],lg[i]=lg[i>>1]+1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
mx[j][i][0]=read();
for(int i=1;i<=m;i++)
for(int k=1;k<20;k++)
for(int j=1;j+(1<<(k-1))<=n;j++)
mx[i][j][k]=max(mx[i][j][k-1],mx[i][j+(1<<(k-1))][k-1]);
dfs(1,n,1,n);
for(int i=1;i<=n;i++)
ans=max(ans,f[i]);
printf("%lld\n",ans);
}