动态规划学习系列三

一：最少拦截系统http://acm.hdu.edu.cn/showproblem.php?pid=1257
    两种做法,一是贪心,从后往前贪;二是DP;
    if(v[i]>max{dp[j]}) (0<=j<len)
    dp[len++]=v[i];    

贪心法：

int i,n,count,a[1005],res,max;  
    while (scanf("%d",&n)!=EOF)  
    {  
        for (i=1; i<=n; ++i)  
        {  
            scanf("%d",&a[i]);  
        }  
        res=0;  
        count=n;  
        while (count)  
        {  
            max=50000;  
            res++;  
            for (i=1; i<=n; ++i)  
            {  
                if (a[i]>0 && a[i]<max)  
                {  
                    max=a[i];  
                    a[i]=-1;  
                    count--;  
                }  
            }  
        }  
        printf("%d/n",res);  

dp法：等价于求最长递增子序列的长度

 vector<int> data;
    vector<int> dp;
  while(scanf("%d",&n)!=EOF){
      data.resize(n+1);
      dp.resize(n+1);
    for(int i=1;i<=n;++i){
        cin>>data[i];
    }
    int len=0;
    dp[0]=INT_MAX;
    for(int i=1;i<=n;++i){
        int flag=0;
        for(int j=0;j<=len;++j){
            if(data[i]<=dp[j]) {dp[j]=data[i];flag=1;break;}
        }
        if(!flag) dp[++len]=data[i];
    }
    cout<<len+1<<endl;
  }

 

二：Common Subsequence http://acm.hdu.edu.cn/showproblem.php?pid=1159 
    经典DP,最长公共子序列
    Len[i][j]={len[i-1][j-1]+1,(a[i]==b[j]); max(len[i-1][j],len[i][j-1])}

 string s1,s2;
  while(cin>>s1>>s2){
      vector<vector<int> > lcs(s1.length()+1,vector<int>(s2.length()+1));
      for(int i=0;i<s1.length();++i){
         for(int j=0;j<s2.length();++j){
            if(s1[i]==s2[j]){
                lcs[i+1][j+1]=lcs[i][j]+1;
            }else{
                lcs[i+1][j+1]=max(lcs[i+1][j],lcs[i][j+1]);
            }
         }
      }
      cout<<lcs[s1.length()][s2.length()]<<endl;
  }

 三：搬寝室http://acm.hdu.edu.cn/showproblem.php?pid=1421 
    状态Dp[i][j]为前i件物品选j对的最优解
    当i=j*2时,只有一种选择即 Dp[i-2][j-1]+(w[i]-w[i-1])^2
    当i>j*2时,Dp[i][j] = min(Dp[i-1][j],Dp[i-2][j-1]+(w[j]-w[j-1])^2) 

const int N=2001;
int data[N];
int dp[N][1001];
int main() {
 //freopen("in.txt","r",stdin);
   int n,k;
  while(cin>>n>>k){
    for(int i=1;i<=n;++i){
        cin>>data[i];
    }
    for(int i=0;i<=n;++i)
    for(int j=1;j<=k;++j)
    dp[i][j]=INT_MAX;
    dp[0][0]=0;
    sort(data+1,data+n+1);
    for(int i=2;i<=n;++i){
        for(int j=1;j*2<=i;++j){
            dp[i][j]=min(dp[i-1][j],dp[i-2][j-1]+(data[i]-data[i-1])*(data[i]-data[i-1]));
        }
    }
    cout<<dp[n][k]<<endl;
  }
    return 0;
}

   四：Humble Numbers http://acm.hdu.edu.cn/showproblem.php?pid=1058 
    如果一个数是Humble Number,那么它的2倍,3倍,5倍,7倍仍然是Humble Number
    定义F[i]为第i个Humble Number
    F[n]=min(2*f[i],3*f[j],5*f[k],7*f[L]), i,j,k,L在被选择后相互移动
    (通过此题理解到数组有序特性)    

五： Doing Homework Again http://acm.hdu.edu.cn/showproblem.php?pid=1789 
    这题为贪心,经典题;
    切题角度,对于每个任务要么在截至日期前完成要么被扣分;所以考虑每个人物的完成情况即可;由于每天只能完成一个任务,所以优先考虑分值较大的任务,看看该任务能不能完成,只要能完成,即使提前完成,占了其他任务的完成日期也没关系,因为当前任务的分值最大嘛,而对于能完成的任务能拖多久就拖多久（最好是截止时间那天完成，如果那天已经有任务了，就从那天向前寻找一个空闲的天）,以便腾出更多时间完成其他任务;    

const int N=1001;
bool vis[N];
struct HW{
    int deadline;
    int score;
    int state;
}a[N];
int cmp(const HW& a,const HW& b){
    return a.score>b.score;
}

int main() {
 //freopen("in.txt","r",stdin);
   int t,n;
   cin>>t ;
   while(t--){
        cin>>n;
        memset(vis,false,sizeof(vis));
        for(int i=0;i<n;++i)
        cin>>a[i].deadline;
         for(int i=0;i<n;++i)
        {cin>>a[i].score;a[i].state=0;}
        sort(a,a+n,cmp);
        int res=0;
        for(int i=0;i<n;++i){
            if(vis[a[i].deadline]==false){
                vis[a[i].deadline]=true;a[i].state=1;
            }else{
                for(int d=a[i].deadline-1;d>=1;--d){
                   if(vis[d]==false){vis[d]=true;a[i].state=1;break;}
                }
            }
            if(a[i].state==0){
                res+=a[i].score;
            }
        }
        cout<<res<<endl;

   }

    return 0;
}

六：How Many Ways http://acm.hdu.edu.cn/showproblem.php?pid=1978 
    两种D法,一是对于当前的点,那些点可达;二是当前点可达那些点;
    明显第二种方法高,因为第一种方法有一些没必要的尝试;
    Dp[i][j]+=Dp[ii][jj]; (map[ii][jj]>=两点的曼哈顿距离)（这里i,j表示能到达的点，ii,和jj对应从发点）
    值得优化的地方,每两点的曼哈顿距离可能不止求一次,所以预处理一下直接读取    

 int dp[101][101];
    int data[101][101];

    int t,n,m;
    cin>>t ;
    while(t--) {
        cin>>n>>m;
        for(int i=0; i<n; ++i)
            for(int j=0; j<m; ++j) {
                cin>>data[i][j];
            }
        memset(dp,0,sizeof(dp));
        dp[0][0]=1;
        for(int i=0; i<n; ++i)
            for(int j=0; j<m; ++j) {
                if(i==n-1&&j==m-1) continue;
                dp[i][j]%=10000;
                for(int x=i; x<=i+data[i][j]&&x<n; ++x)
                    for(int y=j; y<=j+data[i][j]&&y<m; ++y) {
                        if(x==i&&y==j) continue;
                        if(data[i][j]>=x-i+y-j)
                            dp[x][y]+=dp[i][j];
                    }
            }
        dp[n-1][m-1]=dp[n-1][m-1]%10000;
        cout<<dp[n-1][m-1]<<endl;
    }

七：珍惜现在 感恩生活http://acm.hdu.edu.cn/showproblem.php?pid=2191 
    每个物品最多可取n件,多重背包;
    利用二进制思想,把每种物品转化为几件物品,然后就成为了0-1背包    ，直接利用二进制优化转换得到

 int dp[105];
    int v[105];
    int c[105];
    int w[105];

    int t,n,m;
    cin>>t ;
    while(t--) {
        cin>>n>>m;
        for(int i=1; i<=m; ++i) {
            cin>>v[i]>>w[i]>>c[i];
        }
        memset(dp,0,sizeof(dp));
        for(int i=1; i<=m; ++i) {
            int k=1,num=c[i];
            while(k<num) {
                for(int j=n; j>=k*v[i]; --j) {
                    dp[j]=max(dp[j],dp[j-k*v[i]]+k*w[i]);
                }
                num-=k;
                k<<=1;
            }
            for(int j=n; j>=num*v[i]; --j) {
                dp[j]=max(dp[j],dp[j-num*v[i]]+num*w[i]);
            }
        }
        cout<<dp[n]<<endl;
    }

 

八：Piggy-Bank http://acm.hdu.edu.cn/showproblem.php?pid=1114 
    完全背包;常规背包是求最大值,这题求最小值;
    只需要修改一下初始化,f[0]=0,其他赋值为+∞即可;
    状态转移方程:f[i][V]=max{f[i-1][V],f[i-1][V-k*v[i]]+k*w[i]},其中0<=k*v[i]<=V