Hdu 5489 合肥网络赛 1009 Removed Interval

跳跃式LIS(nlogn)，在普通的转移基础上增加一种可以跨越一段距离的转移，用一颗新的树状数组维护，同时，我们还要维护跨越完一次后面的转移，所以我用了3颗树状数组。。

比赛的时候一句话位置写错了，然后就。。。雪崩

呆马：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;
int a[100][100];
int sum,ans,num;
int t,n,m,x,y;
int dp[50][50][3002];
int main()
{
    scanf("%d",&t);
    for (int k=1;k<=t;k++)
    {
        memset(a,0,sizeof(a));
        memset(dp,0,sizeof(dp));

        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
            for (int j=1;j<=m;j++)
                scanf("%d",&a[i][j]);

        dp[1][1][a[1][1]]=a[1][1]*a[1][1];

        for (int i=1;i<=n;i++)
        {
            for (int j=1;j<=m;j++)
            {
                x=i-1; y=j-1;
                if (i==1 && j==1) continue;
                for (int k=1;k<=2000;k++)
                {
                    if (dp[x][j][k]==0) continue;
                    sum=k+a[i][j];
                    num=dp[x][j][k]+a[i][j]*a[i][j];
                    if (dp[i][j][sum]==0) dp[i][j][sum]=num;
                    dp[i][j][sum]=min(dp[i][j][sum],num);
                }            
                for (int k=1;k<=2000;k++)
                {
                    if (dp[i][y][k]==0) continue;
                    sum=k+a[i][j];
                    num=dp[i][y][k]+a[i][j]*a[i][j];
                    if (dp[i][j][sum]==0) dp[i][j][sum]=num;
                    dp[i][j][sum]=min(dp[i][j][sum],num);
                }
                
            }
        }    
        
        ans=999999999;
        for (int i=1;i<=2000;i++)
        {
            if (dp[n][m][i]==0) continue;
            num=(n+m-1)*dp[n][m][i]-i*i;
            //cout<<i<<' '<<dp[n][m][i]<<' '<<num<<endl;
            ans=min(ans,num);
        }
        printf("Case #%d: %d
",k,ans);
    }
}