[NOIP模拟]Array

题目描述:
题目背景: SOURCE：NOIP2015-SHY-9
给定 2 个正整数序列 A1, A2，序列长度分别为 L1, L2。
你可以进行以下的一次操作：
1. 选择两个数 K1，K2（1≤K1≤L1, 1≤K2≤L2）；
2. 移去 A1 中最后 K1 个数，得到这 K1 个数的和 S1，L1 对应减少 K1；
3. 移去 A2 中最后 K2 个数，得到这 K2 个数的和 S2，L2 对应减少 K2；
此次操作的费用为：(S1-K1) * (S2-K2)。
进行以上操作直至两个序列都为空，求最小的费用总和。
注意：序列为空当且仅当两个序列同时为空。
输入格式:
第一行是两个正整数 L1和 L2，表示 A1 与 A2 的长度。
第二行 L1 个整数，表示序列 A1[1..L1]。
第三行 L2 个整数，表示序列 A2[1..L2]。
输出格式:
输出一个整数，表示最小费用。
样例输入：
3 2
1 2 3
1 2
样例输出：
2
样例说明：
第一次选取 K1=1，K2=1。费用为 (3-1)*(2-1) = 2。
第二次选取 K1=2，K2=1。费用为 (1+2-2)*(1-1) = 0。
所以，总费用为 2。
数据范围：
对 20% 的输入数据：1≤L1,L2≤20
对 40% 的输入数据：1≤L1≤400；1≤L2≤150
对 100% 的输入数据：1≤L1,L2,A1[1..L1],A2[1..L2]≤5,000
题目分析:
两个序列记做a[l1],b[l2]。
1、操作费用的计算式为$（sum1-k1）*(sum2-k2)$,就是求和再减去个数，所以我们给原先的每个数都减一，就使得计算式变为$sum1*sum2$；
2、由乘法分配律可知：
$sum1*sum2=(a[i]+···+a[i+k1])*(b[j]+···+b[j+k2])$
$=a[i]*b[j]+···a[i]*b[j+k2]+···a[i+k1]*b[j]+···a[i+k1]*b[j+k2]$
所以你两个序列同时删除几个的贡献，就是你删的那几个对应相乘；
3、因为$(a+b)*(c+d)>a*c+b*d$，所以你在删除时，必须至少有一个队列只删一个才是最优的，否则同时删除几个，用一个一个等效替代删才是最优的。
4、综上所述，其实这是一道动态规划题，设f[i][j]，i表示l1序列剩余i个,j表示l2序列剩j个,存的是此时所需要的最小花费。有动态转移方程：
$f[i][j]=min(f[i+1][j+1],f[i+1][j],f[i][j+1])+a[i+1]*b[j+1]$
最后答案就是f[0][0]。
附代码:

#include<iostream>
#include<cstring>
#include<string>
#include<cstdlib>
#include<cstdio>
#include<ctime>
#include<cmath>
#include<cctype>
#include<iomanip>
#include<algorithm>
using namespace std;

const int N=5e3+10;
int l1,l2,a[N],b[N],f[N][N];

int readint()
{
    char ch;int i=0,f=1;
    for(ch=getchar();(ch<'0'||ch>'9')&&ch!='-';ch=getchar());
    if(ch=='-') {ch=getchar();f=-1;}
    for(;ch>='0'&&ch<='9';ch=getchar()) i=(i<<3)+(i<<1)+ch-'0';
    return i*f;
}

int main()
{
    //freopen("array.in","r",stdin);
    //freopen("array.out","w",stdout);

    l1=readint();l2=readint();
    for(int i=1;i<=l1;i++) a[i]=readint()-1;
    for(int i=1;i<=l2;i++) b[i]=readint()-1;

    memset(f,0x3f3f3f3f,sizeof(f));
    f[l1][l2]=0;
    for(int i=l1-1;i>=0;i--)
        for(int j=l2-1;j>=0;j--)
            f[i][j]=min(f[i+1][j+1],min(f[i+1][j],f[i][j+1]))+a[i+1]*b[j+1];
    printf("%d",f[0][0]);

    return 0;
}