DP练习之——HDU

DP练习之——HDU

前言:虽然做了将近50道题,但是感觉dp这方面还是不怎么样,状压dp,树形dp,插头dp这种比较难的dp问题都跳过了。
总的来说:也还是有所收获的,至少思维有所提升,未来的路还很长,加油!

day 1:

一、HDU 1025:Constructing Roads In JGShining’s Kingdom(题目链接)

题意:
开始看了好久的题意,毫无思路…其实就是求最长上升子序列

题解:
开始直接用的普通的O(n^2)的算法,结果TLE…
于是便利用了贪心的策略,通过二分查找计算最长上升子序列长度

具体思想:
假设要寻找最长上升子序列的序列是a[n],然后寻找到的递增子序列放入到数组b中。

(1)当遍历到数组a的第一个元素的时候,就将这个元素放入到b数组中,以后遍历到的元素都和已经放入到b数组中的元素进行比较;

(2)如果比b数组中的每个元素都大,则将该元素插入到b数组的最后一个元素,并且b数组的长度要加1;

(3)如果比b数组中最后一个元素小,就要运用二分法进行查找,查找出第一个比该元素大的最小的元素,然后将其替换。

在这个过程中,只重复执行这两步就可以了,最后b数组的长度就是最长的上升子序列长度。例如:如该数列为:

5 9 4 1 3 7 6 7

那么:

5 //加入
5 9 //加入
4 9 //用4代替了5
1 9 //用1代替4
1 3 //用3代替9
1 3 7 //加入
1 3 6 //用6代替7
1 3 6 7 //加入

最后b中元素的个数就是最长递增子序列的大小,即4。

要注意的是最后数组里的元素并不就一定是所求的序列,
例如如果输入 2 5 1
那么最后得到的数组应该是 1 5
而实际上要求的序列是 2 5

(即:只能用来计算最长上升子序列的长度,不能用来构造最长上升子序列)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=500000;
int a[maxn+10];
int b[maxn+10];
int n,x,y,b_len;
//二分查找,返回位置
int binSearch(int num){
    int l=1,r=b_len;
    while(l<=r){
        int mid=(l+r)/2;
        if(num>=b[mid])
            l=mid+1;
        else
            r=mid-1;
    }
    return l;
}
int main(){
    int t=1;
    while(scanf("%d",&n)!=EOF){
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        for(int i=1;i<=n;i++){
            scanf("%d%d",&x,&y);
            a[x]=y;
        }
        b_len=1;
        b[1]=a[1];
        for(int i=2;i<=n;i++){
            if(a[i]>b[b_len]){
                b[++b_len]=a[i];
            }
            else{
                int pos=binSearch(a[i]);
                b[pos]=a[i];
            }
        }
        printf("Case %d:\n",t++);

        //还有一个坑点,注意road和roads。好吧,我眼瞎了
        if(b_len==1){
            printf("My king, at most 1 road can be built.\n");
        }
        else{
            printf("My king, at most %d roads can be built.\n",b_len);
        }
        printf("\n");
    }
    return 0;
}

二、HDU 1058:Humble Numbers(题目链接)

题意:
给定一个数n,要你求第n个丑数的值(丑数:即素因子只能为2,3,5,7)

题解:
开始还没想到什么策略,其实只要从2,3,5,7的倍数入手即可。
如果一个数是Humble Numbers,那么它的2倍、3倍、5倍、7倍也都是Humble Numbers
于是状态转移方程为:
dp[n]=min(dp[i]*2,dp[j]*3,dp[k]*5,dp[l]*7)
i,j,k,l在被选择后相互移动

(有个英语知识点,居然让我wa了好几次…)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=5843;
int a[maxn+10];
int n;
//要手写min函数
int min(int num1,int num2,int num3,int num4){
    int minnum=num1;
    if(minnum>num2){
        minnum=num2;
    }
    if(minnum>num3){
        minnum=num3;
    }
    if(minnum>num4){
        minnum=num4;
    }
    return minnum;
}
void solve(){
    int i,j,k,l;
    i=j=k=l=1;
    a[1]=1;
    for(int t=2;t<=maxn;t++){
        a[t]=min(2*a[i],3*a[j],5*a[k],7*a[l]);
        if(a[t]==2*a[i]){
            i++;
        }
        if(a[t]==3*a[j]){
            j++;
        }
        if(a[t]==5*a[k]){
            k++;
        }
        if(a[t]==7*a[l]){
            l++;
        }
    }
}
int main(){
    solve();
    while(~scanf("%d",&n)){
        if(n==0)
            break;
        //需要注意的是,并非只有11,12,13要注意。这种英语都犯错了....
        if(n%10==1&&n%100!=11){
            printf("The %dst humble number is %d.\n",n,a[n]);
        }
        else if(n%10==2&&n%100!=12){
            printf("The %dnd humble number is %d.\n",n,a[n]);
        }
        else if(n%10==3&&n%100!=13){
            printf("The %drd humble number is %d.\n",n,a[n]);
        }
        else{
            printf("The %dth humble number is %d.\n",n,a[n]);
        }
    }
    return 0;
}

三、HDU 1059:Dividing(题目链接)

题意:
每组给出六个数,表示每个值的个数,其中相应的位置表示数的值,同样也表示你相应的代价

题解:
其实就是多重背包的模板题,看网上的题解要用二进制优化,但是我没用二进制优化也过了…

开始真的没想到这是多重背包,题目还是做少了

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=120000;
int volume[maxn+10],value[maxn+10],c[maxn+10];
int dp[maxn+10];
int n,v;
//0-1背包
void ZeroOnepark(int val,int vol) {
    for(int j=v; j>=vol; j--) {
        dp[j]=max(dp[j],dp[j-vol]+val);
    }
}
//完全背包
void Completepark(int val,int vol){
    for(int j=vol;j<=v;j++){
        dp[j]=max(dp[j],dp[j-vol]+val);
    }
}
//多重背包
void Multiplepark(int val,int vol,int amount) {
    if(vol*amount>=v){
        Completepark(val,vol);
    }
    else{
        int k=1;
        while(k<amount) {
            ZeroOnepark(k*val,k*vol);
            amount-=k;
            k<<=1;
        }
        ZeroOnepark(amount*val,amount*vol);
    }
}
int main() {
    int t=1;
    int a1,a2,a3,a4,a5,a6;
    while(~scanf("%d%d%d%d%d%d",&a1,&a2,&a3,&a4,&a5,&a6)) {
        if(a1==0&&a2==0&&a3==0&&a4==0&&a5==0&&a6==0) {
            break;
        }
        int sum=a1*1+a2*2+a3*3+a4*4+a5*5+a6*6;
        if(sum%2) {
            printf("Collection #%d:\n",t++);
            printf("Can't be divided.\n");
            printf("\n");
        } else {
            n=6,v=sum/2;
            for(int i=1; i<=n; i++) {
                value[i]=i;
                volume[i]=i;
            }
            c[1]=a1,c[2]=a2,c[3]=a3,c[4]=a4,c[5]=a5,c[6]=a6;
            memset(dp,0,sizeof(dp));
            for(int i=1; i<=n; i++) {
                Multiplepark(value[i],volume[i],c[i]);
            }
            if(dp[v]!=v) {
                printf("Collection #%d:\n",t++);
                printf("Can't be divided.\n");
                printf("\n");
            } else {
                printf("Collection #%d:\n",t++);
                printf("Can be divided.\n");
                printf("\n");
            }
        }
    }
}

刚好10点,hdu关了…一天才做3道题,还是太懒散了

------------------------------------------------------我是分割线-----------------------------------------------------------------

day 2:

又是给社会做贡献的一天(hhh) 希望疫情早日得以解决,武汉加油,中国加油!!,今天边看题解,边写写画画,终于做了10道题,但是10点hdu关了,只好明天写小结了…

一、HDU 1069:Monkey and Banana(题目链接)

题意:
多组数据。每组数据给定一个n,表示有多少个长方体,其中每个长方体给出对应的长宽高,求最多能叠多高。
要求:上面的长方体的底面长宽均要大于下面的长方体的底面长宽

题解:
其实就是最长递增子序列LIS的长方体模型
由于给定长宽高不唯一,即:可以根据三个数值构造6个长方体,预先进行排序即可

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1000;
//lis的长方体堆叠模型
struct node{
    int l,w,h;     //长、宽、高
};
node a[maxn+10];
int dp[maxn+10];
//底面积升序
bool cmp(node x,node y){
    if(x.l!=y.l){
        return x.l<y.l;
    }
    else{
        return x.w<y.w;
    }
}
int main(){
    int n;
    int t=1;
    while(~scanf("%d",&n)){
        if(n==0){
            break;
        }
        int x,y,z;
        int k=1;
        for(int i=1;i<=n;i++){
            scanf("%d%d%d",&x,&y,&z);
            a[k].l=x,a[k].w=y,a[k].h=z,k++;
            a[k].l=x,a[k].w=z,a[k].h=y,k++;
            a[k].l=y,a[k].w=x,a[k].h=z,k++;
            a[k].l=y,a[k].w=z,a[k].h=x,k++;
            a[k].l=z,a[k].w=x,a[k].h=y,k++;
            a[k].l=z,a[k].w=y,a[k].h=x,k++;
        }
        k-=1;       //由于最后加了一个1
        sort(a+1,a+k+1,cmp);
        for(int i=1;i<=k;i++){
            dp[i]=a[i].h;
        }
        for(int i=1;i<=k;i++){
            for(int j=1;j<=i;j++){
                if(a[j].l<a[i].l&&a[j].w<a[i].w){
                    dp[i]=max(dp[i],dp[j]+a[i].h);
                }
            }
        }
        int maxnum=0;
        for(int i=1;i<=k;i++){
            maxnum=max(maxnum,dp[i]);
        }
        printf("Case %d: maximum height = %d\n",t++,maxnum);
    }
    return 0;
}

二、HDU 1074:Doing Homework(题目链接)

只知道是状压DP,但是状压DP太难了。。。。
不会做啊
先放在这吧,等做完基础的了,再来看

三、HDU 1078:FatMouse and Cheese(题目链接)

题意:
类似于走迷宫。
胖老鼠走迷宫,其中每次最多走k步,迷宫中每个点都有对应的价值,求最多可以获得多少价值,要求:胖老鼠每次停下的地方的价值要大于原先停下来的价值
多组数据,每组数据给出迷宫的大小n(其中长宽相同),k。以及每个位置的价值

题解:
这种类似于迷宫问题,直接用记忆化搜索即可,需要注意的是,对于k这个条件要处理一下

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1000;
int n,k;
int mp[maxn+10][maxn+10];
int dp[maxn+10][maxn+10];
int vis[maxn+10][maxn+10];
int dir[4][2]={0,1,1,0,0,-1,-1,0};
int check(int x,int y){
    if(x<1||x>n||y<1||y>n||vis[x][y]==1)
        return 0;
    else
        return 1;
}
int dfs(int x,int y){
    if(dp[x][y]!=0){
        return dp[x][y];
    }
    int cur=0;
    for(int i=0;i<4;i++){
        //四个方向分别遍历k步
        for(int j=1;j<=k;j++){
            int fx=x+dir[i][0]*j;
            int fy=y+dir[i][1]*j;
            if(check(fx,fy)&&mp[fx][fy]>mp[x][y]){
                vis[fx][fy]=1;
                cur=max(cur,dfs(fx,fy));
                vis[fx][fy]=0;
         	}
        }
    }
    dp[x][y]=cur+mp[x][y];
    return dp[x][y];
}
int main(){
    while(~scanf("%d%d",&n,&k)){
        if(n==-1&&k==-1)
            break;
        memset(mp,0,sizeof(mp));
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                scanf("%d",&mp[i][j]);
            }
        }
        dfs(1,1);
        printf("%d\n",dp[1][1]);
    }
    return 0;
}

四、HDU 1080:Human Gene Functions(题目链接)

题意:
可以理解为DNA匹配。
提前给出碱基相互对应的价值,然后每组数据给出两个字符串,判断最大的匹配值。

题解:
就是洛谷的P1140的改编题
刚好前段时间有做到,

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1000;
int t;
const int table[5][5]={
    {5,-1,-2,-1,-3},
    {-1,5,-3,-2,-4},
    {-2,-3,5,-2,-2},
    {-1,-2,-2,5,-1},
    {-3,-4,-2,-1,0}
};
string sa,sb;
int la,lb;
int a[maxn+10],b[maxn+10];
int dp[maxn+10][maxn+10];
int main(){
    scanf("%d",&t);
    while(t--){
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        cin>>la>>sa;
        cin>>lb>>sb;
        for(int i=1;i<=la;i++){
            for(int j=1;j<=lb;j++){
                dp[i][j]=-inf;
            }
        }
        for(int i=1;i<=la;i++){
            if(sa[i-1]=='A'){
                a[i]=0;
            }
            else if(sa[i-1]=='C'){
                a[i]=1;
            }
            else if(sa[i-1]=='G'){
                a[i]=2;
            }
            else if(sa[i-1]=='T'){
                a[i]=3;
            }
        }
        for(int i=1;i<=lb;i++){
            if(sb[i-1]=='A'){
                b[i]=0;
            }
            else if(sb[i-1]=='C'){
                b[i]=1;
            }
            else if(sb[i-1]=='G'){
                b[i]=2;
            }
            else if(sb[i-1]=='T'){
                b[i]=3;
            }
        }
        for(int i=1;i<=la;i++){
            dp[i][0]=dp[i-1][0]+table[a[i]][4];
        }
        for(int i=1;i<=lb;i++){
            dp[0][i]=dp[0][i-1]+table[b[i]][4];
        }
        for(int i=1;i<=la;i++){
            for(int j=1;j<=lb;j++){
                dp[i][j]=max(dp[i][j],dp[i][j-1]+table[b[j]][4]);
                dp[i][j]=max(dp[i][j],dp[i-1][j]+table[a[i]][4]);
                dp[i][j]=max(dp[i][j],dp[i-1][j-1]+table[a[i]][b[j]]);
            }
        }
        printf("%d\n",dp[la][lb]);
    }
    return 0;
}

五、HDU 1081:To The Max(题目链接)

题意:
n*n矩阵中,求最大子矩阵

题解:
直接套模板,需要注意的是long long

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1000;
int n,m;
int a[maxn+10][maxn+10];
int dp[maxn+10][maxn+10];
int main(){
    while(~scanf("%d",&n)){
        m=n;
        memset(dp,0,sizeof(dp));
        memset(a,0,sizeof(a));
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&a[i][j]);
                a[i][j]+=a[i-1][j];
            }
        }
        ll sum,ans=0;
        for(int i=1;i<=n;i++){
            for(int j=i;j<=n;j++){
                sum=0;
                for(int k=1;k<=m;k++){
                    sum+=(a[j][k]-a[i-1][k]);
                    if(sum<0)
                        sum=0;
                    if(sum>ans){
                        ans=sum;
                    }
                }
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}

六、HDU 1087:Super Jumping! Jumping! Jumping!(题目链接)

题意:
多组数据,每组数据给出n,表示数组的长度,就是求最长上子序列LIS

题解:
这里不是求最长上升子序列长度,而是求和,只要稍微转换下即可

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1000;
int n;
int a[maxn+10];
int dp[maxn+10];
int main(){
    while(~scanf("%d",&n)&&n){
        memset(a,0,sizeof(a));
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
            dp[i]=a[i];
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<i;j++){
                if(a[i]>a[j]){
                    dp[i]=max(dp[i],dp[j]+a[i]);
                }
            }
        }
        int ans=0;
        for(int i=1;i<=n;i++){
            ans=max(ans,dp[i]);
        }
        printf("%d\n",ans);
    }
    return 0;
}

七、HDU 1114:Piggy-Bank(题目链接)

题意:
有一个储钱罐,给出储钱罐自身的重量,和装了钱的重量
然后有n种硬币,每种硬币有对应的价值和重量
求:储钱罐最少有多少钱

题解:
发现就是完全背包的裸题,需要注意的是:一般是求的最大值,这里求最小值,于是只需要对dp[]数组预处理为无穷大即可,然后dp[0]=0;

代码是利用的滚动数组

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=10000;
int n,v;
int volume[maxn+10];
int value[maxn+10];
int dp[maxn+10];
int main(){
    int t;
    int w1,w2;
    scanf("%d",&t);
    while(t--){
        memset(dp,inf,sizeof(dp));     //由于求最小值,所有都置为无穷大
        memset(value,0,sizeof(value));
        memset(volume,0,sizeof(volume));
        scanf("%d%d",&w1,&w2);
        v=w2-w1;
        scanf("%d",&n);
        for(int i=1;i<=n;i++){
            scanf("%d%d",&value[i],&volume[i]);
        }
        
        dp[0]=0;     //这个非常重要
        for(int i=1;i<=n;i++){
            for(int j=volume[i];j<=v;j++){
                if(dp[j-volume[i]]!=inf)
                    dp[j]=min(dp[j],dp[j-volume[i]]+value[i]);
            }
        }
        if(dp[v]!=inf){
            printf("The minimum amount of money in the piggy-bank is %d.\n",dp[v]);
        }
        else{
            printf("This is impossible.\n");
        }
    }
    return 0;
}

八、HDU 1121:Complete the Sequence(题目链接)

这题只能说比较难,为什么还说难度是1。可能是我太菜了吧,看了题解发现是差分,也就是牛顿插值法,虽然学过,但是不会编程啊。。。
看的大佬的代码:

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=10000;
int s,c,a[maxn][maxn];
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d",&s,&c);
        for(int i=0;i<s;i++){
            scanf("%d",&a[0][i]);
        }
        for(int i=1;i<s;i++){
            for(int j=0;j<s-i;j++){
                a[i][j]=a[i-1][j+1]-a[i-1][j];
            }
        }
        for(int i=1;i<=c;i++){
            a[s-1][i]=a[s-1][0];
        }
        for(int i=s-2;i>=0;i--){
            int pos=s-i;
            for(int j=0;j<c;j++){
                a[i][j+pos]=a[i+1][j+pos-1]+a[i][j+pos-1];
            }
        }
        printf("%d",a[0][s]);
        for(int i=1;i<c;i++){
            printf(" %d",a[0][s+i]);
        }
        printf("\n");
    }
    return 0;
}

九、HDU 1158:Employment Planning(题目链接)

为什么这也是难度1。。。

看了网上的题解后,自己写,然后一直RE。。。。
现在都还迷迷糊糊的

#include 
#define ll long long
using namespace std;
int dp[20][500];//*第i阶段保留j个人;
int hiring , salary , firing; //*雇佣工资,月薪,解聘工资;
int num[13];//*阶段人数
int main()
{
    int month ;
    while(cin>>month, month)
    {
        memset(dp,0,sizeof(dp));
        memset(num , 0  ,sizeof(num));
        int maxi = - 1 ;
        cin>>hiring>>salary>>firing;
        for(int i = 1 ; i <=month ;i++)
        {
            cin>>num[i];//*每个阶段的人数;
            if(maxi<num[i])
            {
                maxi = num[i]; //*寻求最大雇佣人数 ;
            }
        }
 
        for(int i= num[1] ; i<=maxi;i++) //*更新初始状态;
        {
            dp[1][i] = i*hiring + i*salary; //*第一个月基础工资
        }
        int mini , cost ;
        for(int i= 2 ; i<= month ; i++)
        {
            for(int j = num[i] ; j<=maxi ; j++)
            {
                mini = 9999999 ;
                for(int k = num[i-1] ; k <=maxi ; k++)
                {
                    if(k<=j)
                        cost = dp[i-1][k]+(j-k)*hiring+j*salary;
                    else cost =  dp[i-1][k]+(k-j)*firing+j*salary;
                    if(mini >cost )
                        mini = cost ;
                }
                dp[i][j] = mini ;
            }
        }
        mini = 99999999 ;
        for(int i = num[month]; i<=maxi;i++)
        {
            if(dp[month][i] < mini )
            {
                mini = dp[month][i] ;
            }
        }
 
        cout << mini <<endl;
    }
    return 0 ;
}

十、HDU 1159:Common Subsequence(题目链接)

这题相比于上题有太简单了,就是求最长公共子序列的长度
也就没什么要多说的了,非常经典的模板题

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1000;
string s1,s2;
int dp[maxn+10][maxn+10];
int main(){
    while(cin>>s1>>s2){
        int len1=s1.length();
        int len2=s2.length();
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=len1;i++){
            for(int j=1;j<=len2;j++){
                if(s1[i-1]==s2[j-1]){
                    dp[i][j]=dp[i-1][j-1]+1;
                }
                else{
                    dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
                }
            }
        }
        printf("%d\n",dp[len1][len2]);
    }
    return 0;
}

做了这么多,发现简单的一维、二维的DP问题目前还是比较熟悉了
至于树形DP、状压DP等刷完简单的了,再好好恶补一下吧

day3:

今天有点懒散。。。

一、HDU 1160:FatMouse’s Speed(题目链接)

题意:
n对数,每对数有重量和速度组成,要求:在保证重量递增,速度递减的条件下;求最多能有多少对组成一序列

题解:
其实就是先排序后,求最长上升序列。需要注意的是:需要打印路径,于是结构体中要保存位置,同时还是用pre[]存储路径,直接用递归即可打印路径

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=10000;
string s1,s2;
struct node{
    int w,v;       //体重和速度
    int num;       //第几个数
};
node a[maxn+10];
int dp[maxn+10];
int pre[maxn+10];     //记录路径
//速度降序,体重升序
bool cmp(node x,node y){
    if(x.w!=y.w){
        return x.w<y.w;
    }
    else{
        return x.v>y.v;
    }
}
//打印路径
void dfs(int pos,int num){
    if(num==0)
        return ;
    dfs(pre[pos],num-1);
    printf("%d\n",a[pos].num);
}
int main(){
    int t=1;
    while(scanf("%d%d",&a[t].w,&a[t].v)!=EOF){
        dp[t]=1;
        pre[t]=0;
        a[t].num=t;
        t++;
    }
    t--;
    sort(a+1,a+t+1,cmp);
    for(int i=1;i<=t;i++){
        for(int j=1;j<i;j++){
            if(a[i].w>a[j].w&&a[i].v<a[j].v){
                int temp=dp[j]+1;
                if(temp>dp[i]){
                    dp[i]=temp;
                    pre[i]=j;
                }
            }
        }
    }
    int pos=0,maxnum=0;
    for(int i=1;i<=t;i++){
        if(dp[i]>maxnum){
            maxnum=dp[i];
            pos=i;
        }
    }
    printf("%d\n",maxnum);
    dfs(pos,maxnum);
    return 0;
}

二、HDU 1165:Eddy’s research II(题目链接)

题意:
在这里插入图片描述
给定n,m求A(m,n)

题解:
找规律,硬是用手算了好久,直接见代码,感觉如果比赛遇到这种题,还是有点逼心态的

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1000000;
ll n,m;
ll dfs(ll m,ll n){
    if(n==0){
        return dfs(m-1,1);
    }
    else if(m==1){
        return n+2;
    }
    else if(m==2){
        return 2*n+3;
    }
    else if(m==3){
        return 2*dfs(3,n-1)+3;
    }
}
int main(){
    while(~scanf("%lld%lld",&m,&n)){
        printf("%lld\n",dfs(m,n));
    }
    return 0;
}

三、HDU 1176:免费馅饼(题目链接)

题意:
开始站在位置为5的位置,每次能够接到他左右一个位置和他现在的位置的馅饼,题目会给出一段时间类,具体某个位置某个时间会掉馅饼,求最多可以获得多少个馅饼

题意:
其实还是比较容易想到状态转移方程的,但是由于正向dp不好处理,于是采用反向dp
dp[i][j]=max(dp[i+1][j],dp[i+1][j+1],dp[i+1][j-1])+a[i][j];
(其中i表示时间,j表示位置)

又由于位置0不好处理,于是直接每个位置+1
最后只需要输出dp[0][6]即可

具体见代码:

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=100000;
int dp[maxn+10][20];
int a[maxn+10][20];
int n;
int main(){
    while(~scanf("%d",&n)){
        if(n==0)
            break;
        int maxt=0;
        int pos,time;
        memset(a,0,sizeof(a));
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            scanf("%d%d",&pos,&time);
            //由于0不好处理,所以++
            a[time][++pos]++;
            maxt=max(maxt,time);
        }
        /*for(int i=1;i<=t;i++){
            for(int j=1;j<=10;j++){
                printf("%d ",a[i][j]);
            }
            printf("\n");
        }*/
        int temp1,temp2,temp3;
        //反向dp
        for(int i=maxt;i>=0;i--){
            for(int j=1;j<=11;j++){
                dp[i][j]=max(max(dp[i+1][j],dp[i+1][j-1]),dp[i+1][j+1])+a[i][j];
            }
        }
        //因为起点为5,又位置加一了,所以就输出dp[0][6]了
        printf("%d\n",dp[0][6]);
    }
    return 0;
}

四、HDU 1203:I NEED A OFFER!(题目链接)

题意:
n元钱,m个学校可以选择。每个学校都有对应的费用和能被录取上的概率,
至少得到一个offer的最大概率

题解:
乍一眼看发现就是一个0-1背包。但是需要注意的是:题目求至少得到一个offer的最大概率p1,如果直接dp,结果其实是错的
于是不如逆向思考,先转化为一个offer都没有得到的最小概率p2,然后
p1=1-p2

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=100000;
double value[maxn+10];
int volume[maxn+10];
double dp[maxn+10];
int n,v;
int main(){
    while(~scanf("%d%d",&v,&n)){
        if(n==0&&v==0){
            break;
        }
        for(int i=1;i<=n;i++){
            scanf("%d%lf",&volume[i],&value[i]);
        }
        //fill对数组整体附一样的值:1
        fill(dp,dp+maxn,1);

        //由于求至少收到一份offer的最大概率p1,不如转化为没有收到offer的最小概率p2
        //p1=1-p2;
        for(int i=1;i<=n;i++){
            for(int j=v;j>=volume[i];j--){
                dp[j]=min(dp[j],dp[j-volume[i]]*(1-value[i]));
            }
        }
        printf("%.1f%%\n",(1-dp[v])*100);
    }
    return 0;
}

五、HDU 1208:Pascal’s Travels(题目链接)

题意:
一个n*n的迷宫,每个位置都有对应的值,表示接下来走几步。题目要求从左上角走到右下角,只能向右和向下走

题解:
同样是迷宫类似问题,dp和记忆化搜索都可以,但是开始出现了点问题,看了看网上的博客才ac

记忆化搜索和dp都写了一下

dp数组版本:
注意dp[1][1]=1;不然没结果

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=100;
int n;
char mp[maxn][maxn];
int a[maxn][maxn];    //将字符转化为数字存储
ll dp[maxn][maxn];
int main(){
    while(~scanf("%d",&n)){
        if(n==-1)
            break;
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            scanf("%s",mp[i]+1);
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                a[i][j]=mp[i][j]-'0';
            }
        }
        //初始化非常重要
        dp[1][1]=1;
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(a[i][j]==0)
                    continue;
                if(j+a[i][j]<=n)
                    dp[i][j+a[i][j]]+=dp[i][j];
                if(i+a[i][j]<=n)
                    dp[i+a[i][j]][j]+=dp[i][j];
            }
        }
        printf("%lld\n",dp[n][n]);
    }
    return 0;
}

记忆化搜索版本:

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=100;
int n;
char mp[maxn][maxn];
int a[maxn][maxn];    //将字符转化为数字存储
ll dp[maxn][maxn];
int dir[2][2]={0,1,1,0};
int check(int x,int y){
    if(x<1||x>n||y<1||y>n)
        return 0;
    else
        return 1;
}
ll dfs(int x,int y){
    if(dp[x][y]){
        return dp[x][y];
    }
    else if(x==n&&y==n){
        return 1;
    }
    else if(a[x][y]==0){
        return 0;
    }
    for(int i=0;i<2;i++){
        int fx=x+dir[i][0]*a[x][y];
        int fy=y+dir[i][1]*a[x][y];
        if(check(fx,fy)){
            dp[x][y]+=dfs(fx,fy);
        }
    }
    return dp[x][y];
}
int main(){
    while(~scanf("%d",&n)){
        if(n==-1)
            break;
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            scanf("%s",mp[i]+1);
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                a[i][j]=mp[i][j]-'0';
            }
        }
        printf("%lld\n",dfs(1,1));
    }
    return 0;
}

六、HDU 1224:Free DIY Tour(题目链接)

题意:
给出一个有向图,n个顶点,每个顶点有对应的价值,m条边构成有向图
要求从1开始最终回到1,问最多可以获得多少价值,同时打印路径

题解:
开始觉得是不是对最短路加以修改就行,但是没有弄出来,最终还是看的网上的博客,才发现原来可以这么弄,害,还是要加油

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=100;
int t,n,m;     //m表示边数
int value[maxn+10];          // 每个点对应的价值
int way[maxn+10];            //记录路径
int mp[maxn+10][maxn+10];    //存有向图的邻接矩阵
int dp[maxn+10];    //dp数组
stack<int>sta;
int main(){
    scanf("%d",&t);
    for(int cnt=1;cnt<=t;cnt++){
        memset(value,0,sizeof(value));
        memset(mp,0,sizeof(mp));
        memset(dp,0,sizeof(dp));
        memset(way,0,sizeof(way));
        scanf("%d",&n);
        for(int i=1;i<=n;i++){
            scanf("%d",&value[i]);
        }
        scanf("%d",&m);
        int x,y;
        for(int i=1;i<=m;i++){
            scanf("%d%d",&x,&y);
            mp[x][y]=1;
        }
        for(int i=1;i<=n+1;i++){
            for(int j=1;j<i;j++){
                if(mp[j][i]&&dp[j]+value[i]>dp[i]){
                    dp[i]=dp[j]+value[i];
                    way[i]=j;
                }
            }
        }
        printf("CASE %d#\n",cnt);
        printf("points : %d\n",dp[n+1]);
        printf("circuit : ");

        int pos=n+1;
        while(way[pos]){
            sta.push(way[pos]);
            pos=way[pos];
        }
        while(!sta.empty()){
            printf("%d->",sta.top());
            sta.pop();
        }
        printf("1\n");
        if(cnt!=t){
            printf("\n");
        }
    }
}

day4:

下午和晚上有点事,就上午做了两道题。。。。

一、HDU 1227:Fast Food

题意:
n个餐馆,k个仓库。给定n个餐馆的位置,其中仓库可以位于任意k个餐馆的位置,求餐馆到达最近仓库的距离和

题解:
说实话,没怎么搞懂。。。。。

ps:在i到j取一点使它到区间每一点的距离之和最小,这一点为(i+j)/2用图形即可证明;

Dp[i][j]=max{Dp[i-1][k]+cost[k+1][j] 其中,(i-1)<=k

#include 
#include 
#include 
#include 
#include 
using namespace std;
int n,k;
int dp[32][210];//dp[i][j]前j个商店有i个仓库
int dis[210];
const int inf=0x3f3f3f3f;
 
int cost(int i,int j)
{
    int ans=0;
    int mid=dis[(i+j)/2];
    for(int k=i;k<=j;k++)
        ans+=abs(dis[k]-mid);
    return ans;
}
 
int main()
{
    int c=1;
    while(scanf("%d%d",&n,&k)!=EOF)
    {
        if(!n||!k)
            break;
        for(int i=1;i<=n;i++)
            scanf("%d",&dis[i]);
        memset(dp,inf,sizeof(dp));
        for(int i=1;i<=n;i++)
            dp[1][i]=cost(1,i);
        for(int i=2;i<=k;i++)//第i个仓库
            for(int j=1;j<=n;j++)//前j个商店
            {
                for(int k=i-1;k<=j-1;k++)
                    dp[i][j]=min(dp[i][j],dp[i-1][k]+cost(k+1,j));
            }
 
        printf("Chain %d\n",c++);
        printf("Total distance sum = %d\n",dp[k][n]);
        printf("\n");
    }
    return 0;
}

二、HDU 1243:反恐训练营

题意:
给定n表示子弹的类型和恐怖分子类型(其中:子弹类型必须与恐怖分子类型相同,才能获得恐怖分子的分数),然后给出恐怖分子序列,以及:每种类型的对应得分。然后给出子弹的类型,求最多能获得多少得分

题解:
居然开始没想到,其实就是LCS(最长公共子序列)模型啊
只是每个字符都有价值而已

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=2000;
int value[maxn+10];             //每个恐怖分子类型对应的得分数
int dp[maxn+10][maxn+10];
char s1[maxn+10],s2[maxn+10];
char s[maxn+10];                //恐怖分子类型
int main(){
    int n;
    while(~scanf("%d",&n)){
        scanf("%s",s);
        for(int i=0;i<n;i++){
            scanf("%d",&value[s[i]-'a']);
        }
        scanf("%s",s1);          //子弹序列
        scanf("%s",s2);          //恐怖分子序列
        int len1=strlen(s1);
        int len2=strlen(s2);
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=len1;i++){
            for(int j=1;j<=len2;j++){
                if(s1[i-1]==s2[j-1]){
                    dp[i][j]=dp[i-1][j-1]+value[s1[i-1]-'a'];
                }
                else{
                    dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
                }
            }
        }
        printf("%d\n",dp[len1][len2]);
    }
    return 0;
}

day 5:

一、HDU 1300:Pearls

题意:
n种珠宝,每种珠宝有对应的数量和价格,按照题目的规则求最终买入所有珠宝的总费用最小

题解:
Dp[i]=min{Dp[j]+V}, 0<=j

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1000;
int num[maxn];
int price[maxn];
int sum[maxn];     //数量前缀和
int dp[maxn];
//dp[i]=min{dp[j]+v}   v表示第j+1到第i个都以第i个价格买入的值
int main(){
    int t,n;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        memset(dp,0,sizeof(dp));
        memset(sum,0,sizeof(sum));
        for(int i=1;i<=n;i++){
            scanf("%d%d",&num[i],&price[i]);
            sum[i]=sum[i-1]+num[i];
            dp[i]=(sum[i]+10)*price[i];
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<i;j++){
                dp[i]=min(dp[i],dp[j]+(sum[i]-sum[j]+10)*price[i]);
            }
        }
        printf("%d\n",dp[n]);
    }
    return 0;
}

二、HDU 1331:Function Run Fun

题意:
求w(a,b,c)
然后规则是:
DP练习之——HDU_第1张图片

题解:
很明显直接递归肯定会tle
记忆化搜索或者直接用dp数组存都可,但是不知道为什么我的记忆化搜索tle了。。。。

于是只给出dp数组方式代码(直接预处理即可)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=110;
ll dp[maxn][maxn][maxn];
int main(){
    memset(dp,0,sizeof(dp));
    for(int i=0;i<=20;i++){
        for(int j=0;j<=20;j++){
            for(int k=0;k<=20;k++){
                if(i==0||j==0||k==0)
                    dp[i][j][k]=1;
                else if(i<j&&j<k){
                    dp[i][j][k]=dp[i][j][k-1]+dp[i][j-1][k-1]-dp[i][j-1][k];
                }
                else{
                    dp[i][j][k]=dp[i-1][j][k]+dp[i-1][j-1][k]+dp[i-1][j][k-1]-dp[i-1][j-1][k-1];
                }
            }
        }
    }
    ll a,b,c;
    while(~scanf("%lld%lld%lld",&a,&b,&c)){
        if(a==-1&&b==-1&&c==-1)
            break;
        printf("w(%lld, %lld, %lld) = ",a,b,c);
        if(a<=0||b<=0||c<=0){
            printf("1\n");
        }
        else if(a>20||b>20||c>20){
            printf("%lld\n",dp[20][20][20]);
        }
        else{
            printf("%lld\n",dp[a][b][c]);
        }
    }
    return 0;
}

三、HDU 1337:The Drunk Jailer

题意:
不知道为什么hdu把它放到dp题

题意:
直接求n的平方根即为答案

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=110;
int t;
int n;
int main(){
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        printf("%d\n",(int)sqrt((double)n));
    }
    return 0;
}

四、HDU 1355:The Peanuts

题意:
感觉像是做过的题目,给出n*m的矩阵,矩阵每个点有对应价值,然后矩阵上方有一道路,每次必须优先选择价值最多的经过,而且或者该点的价值需要一单位时间,走一单位距离也需要一单位时间,问最终在给定允许的时间内,能获得的最多价值

题解:
感觉还是不是dp啊,应该是贪心或者是模拟啊

直接按照优先选择最大的价值经过,于是代码就出来了

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=2500;
int t,n,m,k;
struct node{
    int x,y;
    int value;
};
node a[maxn+10];
int vis[maxn+10];       //记录点是否被访问
bool cmp(node xx,node yy){
    return xx.value>yy.value;
}
int main(){
    scanf("%d",&t);
    while(t--){
        int cnt=1;
        int num;
        memset(a,0,sizeof(a));
        memset(vis,0,sizeof(vis));
        scanf("%d%d%d",&n,&m,&k);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&num);
                if(num){
                    a[cnt].x=i;
                    a[cnt].y=j;
                    a[cnt].value=num;
                    cnt++;
                }
            }
        }
        cnt--;        //由于最后还自加了一次
        sort(a+1,a+cnt+1,cmp);

        /*for(int i=1;i<=cnt;i++){
            printf("%d %d %d\n",a[i].x,a[i].y,a[i].value);
        }*/

        int flag=0;
        node now;                  //目前访问的点
        now.x=0,now.y=a[1].y;
        int dis,cur=1;            //dis表示当前访问的点与下一点的距离,cur记录访问的点
        int sum=0;        //步数之和
        int valueSum=0;     //最终能获得的价值之和
        while(!flag){
            dis=abs(a[cur].x-now.x)+abs(a[cur].y-now.y);
            if(sum+dis+a[cur].x+1<=k){
                sum+=dis+1;        //+1是因为吃草也要时间
                vis[cur]=1;
                valueSum+=a[cur].value;
                now=a[cur];
                cur++;
            }
            else{
                flag=1;
            }
        }
        printf("%d\n",valueSum);
    }
    return 0;
}

五、HDU 1400:Mondriaan’s Dream

题意:
给定一大矩阵,求通过若干1*2的子矩阵能有多少种组合方式

题解:
据说是非常经典的状压dp题,但是我不知道做(呜呜呜…

放一个大佬博客在这,等到时候强攻状压dp的时候再好好理解

https://blog.csdn.net/hopeztm/article/details/7841917

day 6:

虽然做了好几道题,但是都不是自己开始就想到的,dp好难啊

一、HDU 1421:搬教室

题意:
n件物品,但是只需要搬2*k件过去
题目给出每个物品的重量,每次能搬两件,受到的疲劳度为两者的差的平方
求最少能受到的疲劳度

题解:
状态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[i]-w[i-1])^2)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=2000;
string str;
int n,k;
int dp[maxn+10][maxn+10];
int w[maxn+10];
/*
状态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[i]-w[i-1])^2)
*/
int main(){
    while(~scanf("%d%d",&n,&k)){
        for(int i=1;i<=n;i++){
            scanf("%d",&w[i]);
        }
        sort(w+1,w+n+1);
        memset(dp,0,sizeof(dp));
        for(int i=2;i<=n;i++){
            for(int j=1;j<=k;j++){
                if(i==2*j){
                    dp[i][j]=dp[i-2][j-1]+(w[i]-w[i-1])*(w[i]-w[i-1]);
                }
                else if(i>2*j){
                    dp[i][j]=min(dp[i-1][j],dp[i-2][j-1]+(w[i]-w[i-1])*(w[i]-w[i-1]));
                }
            }
        }
        printf("%d\n",dp[n][k]);
    }
    return 0;
}

二、HDU 1422:重温世界杯

题意:
n个城市,每个城市都对应有着花费和生活费,而且最初也有一定的生活费,问最多能逛多少个城市

题解:
1.如果剩下的钱(res)加上当前城市的生活费(w[i])>=当前城市的消费(l[i])
dp[i]=dp[i-1]+1;
res=res+w[i]-l[i];
2.否则剩下的钱置为0

同时加一个maxpos==n时,跳出循环

需要注意的是:这里是个环形的。。。。,所以数组开为两倍即可

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=100000;
string str;
int n;
int w[2*maxn+10];
int l[2*maxn+10];
int dp[2*maxn+10];
/*
1.如果剩下的钱(res)加上当前城市的生活费(w[i])>=当前城市的消费(l[i])
    dp[i]=dp[i-1]+1;
    res=res+w[i]-l[i];
2.否则剩下的钱置为0

同时加一个maxpos==n时,跳出循环
*/
int main(){
    while(~scanf("%d",&n)){
        //由于会形成环,所以相应的n+i位置等于i位置的数值
        for(int i=1;i<=n;i++){
            scanf("%d%d",&w[i],&l[i]);
            w[n+i]=w[i];
            l[n+i]=l[i];
        }
        memset(dp,0,sizeof(dp));
        int res=0,maxpos=0;            //分别表示剩下的钱和最多能参观的城市数
        for(int i=1;i<=2*n;i++){
            if(res+w[i]>=l[i]){
                dp[i]=dp[i-1]+1;
                res=res+w[i]-l[i];
            }
            else{
                res=0;
            }
            //直接一次循环即可得最终结果
            if(dp[i]>maxpos){
                maxpos=dp[i];
            }
            if(dp[i]==n){
                break;
            }
        }
        printf("%d\n",maxpos);
    }
    return 0;
}

三、HDU 1423:Greatest Common Increasing Subsequence

题意:
求最长公共递增子序列
这里其实降低了难度,只要求长度,没有要求构造

题解:
到网上看到了一个非常巧的代码,主要是代码短小精悍

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=500;
int n,m;
int t;
int a[maxn+10],b[maxn+10];
int dp[maxn+10];
int main(){
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        for(int i=0;i<n;i++){
            scanf("%d",&a[i]);
        }
        scanf("%d",&m);
        for(int i=0;i<m;i++){
            scanf("%d",&b[i]);
        }
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            int maxnum=0;
            for(int j=1;j<=m;j++){
                if(a[i-1]==b[j-1]){
                    dp[j]=maxnum+1;
                }
                if(a[i-1]>b[j-1])
                    maxnum=max(dp[j],maxnum);
            }
        }
        int ans=0;
        for(int i=1;i<=m;i++){
            ans=max(dp[i],ans);
        }
        printf("%d\n",ans);
        if(t)
            printf("\n");
    }
    return 0;
}

四、HDU 1428:漫步校园

题意:
n*n的地图,要求从走上角走到右下角。每个位置都有对应的代价,求有多少条不同的最短路径。

题解:
发现这题用简单的dp不能解决,但是自己还是无能为力,只好看了网上的题解,最短路+dp

然后用bfs处理最短路,记忆化搜索获得最短路径数

自己还是很菜啊!!!

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=500;
int n;
ll mp[55][55];
ll dp[55][55];
ll dis[55][55];
int vis[55][55];
int dir[4][2]={0,1,1,0,0,-1,-1,0};
struct node{
    ll x,y;
};
int check(ll x,ll y){
    if(x<1||x>n||y<1||y>n)
        return 0;
    else
        return 1;
}
void bfs(){
    queue<node>q;
    node now,next;
    now.x=n,now.y=n;
    dis[n][n]=mp[n][n];
    q.push(now);
    while(!q.empty()){
        now=q.front();
        q.pop();
        for(int i=0;i<4;i++){
            ll fx=now.x+dir[i][0];
            ll fy=now.y+dir[i][1];
            next.x=fx,next.y=fy;
            if(check(fx,fy)){
                if(dis[fx][fy]>dis[now.x][now.y]+mp[fx][fy]){
                    dis[fx][fy]=dis[now.x][now.y]+mp[fx][fy];
                    q.push(next);
                }
            }
        }
    }
}
ll DP(int x,int y){
    if(vis[x][y])
        return dp[x][y];
    vis[x][y]=1;
    for(int i=0;i<4;i++){
        ll fx=x+dir[i][0];
        ll fy=y+dir[i][1];
        if(check(fx,fy)){
            if(dis[x][y]>dis[fx][fy]){
                dp[x][y]+=DP(fx,fy);
            }
        }
    }
    return dp[x][y];
}
int main(){
    while(~scanf("%d",&n)){
        memset(mp,0,sizeof(mp));
        memset(dp,0,sizeof(dp));
        memset(vis,0,sizeof(vis));
        memset(dis,inf,sizeof(dis));    //因为求最短路,所以距离数组置为最大值
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                scanf("%lld",&mp[i][j]);
            }
        }
        bfs();
        dp[n][n]=1;
        vis[n][n]=1;
        printf("%lld\n",DP(1,1));
    }
}

五、HDU 1466:计算直线的交点数

题意:
n条直线,问有多少条相交方案

题解:
感觉这个dp状态转移方程还是比较难想。。。

dp[直线的总数][交点的个数] = 状态(本状态存在为1,不存在为0)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=200;            //当n=20时,交点数最多也为n*(n-1)/2 < 200
int dp[25][200];          //dp[直线的总数][交点的个数] = 状态(本状态存在为1,不存在为0)
int main(){
    //初始化所有直线都平行
    for(int i=0;i<=20;i++){
        dp[i][0]=1;
    }
    for(int i=2;i<=20;i++){
        for(int j=1;j<i;j++){          //枚举与i相交的边的数目
            for(int k=0;k<200;k++){    //枚举j条边的交点情况
                if(dp[j][k]){           //如果存在,则推论成立
                    dp[i][(i-j)*j+k]=1;
                }
            }
        }
    }
    int n;
    while(~scanf("%d",&n)){
        for(int i=0;i<maxn;i++){
            if(dp[n][i]){
                if(i==0){
                    printf("%d",i);
                }
                else{
                    printf(" %d",i);
                }
            }
        }
        printf("\n");
    }
    return 0;
}

六、HDU 1494:跑跑卡丁车

一道之前校内选拔赛老师选的一道题,结果还是没有做出来。。。。

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=100000;
int a[maxn+10],b[maxn+10],dp[maxn+10][20];
int main(){
    int l,n;
    while(~scanf("%d%d",&l,&n)){
        for(int i=1;i<=l;i++){
            scanf("%d",&a[i]);
            for(int j=1;j<=n;j++){
                a[i+j*l]=a[i];
            }
        }
        for(int i=1;i<=l;i++){
            scanf("%d",&b[i]);
            for(int j=1;j<=n;j++){
                b[i+j*l]=b[i];
            }
        }
        memset(dp,inf,sizeof(dp));
        dp[0][0]=0;
        for(int i=1;i<=l*n;i++){
            for(int j=0;j<15;j++){
                if(j!=0){
                    dp[i][j]=min(dp[i][j],dp[i-1][j-1]+a[i]);
                }
                if(j<10){
                    dp[i][j]=min(dp[i][j],dp[i-1][j+5]+b[i]);
                }
                if(j==10){
                    dp[i][j]=min(dp[i][j],dp[i-1][14]+a[i]);
                }
            }
        }
        int ans=dp[l*n][0];
        for(int i=0;i<15;i++){
            ans=min(ans,dp[l*n][i]);
        }
        printf("%d\n",ans);
    }
    return 0;
}

七、HDU 1501:Zipper

题意:
给定两个短字符串,和一个长字符串
判断长字符串是不是有这两个短字符串组成

题解:
开始还天真的觉得只要判断两个短字符串分别与长字符串求最长公共子序列即可,如果两者的结果之和等于长字符串的长度则正确yes,否则no

其实dp状态应该为

if(dp[i][j-1]&&s2[j]==str[i+j]){
       dp[i][j]=1;
}
if(dp[i-1][j]&&s1[i]==str[i+j]){
      dp[i][j]=1;
}

最后只需要判断dp[len1][len2]是不是等于1

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=100000;
char s1[210],s2[210],str[420];    //记得str数组长度为两倍
int dp[210][210];
//dp[i][j]   表示从s1中选择i个,s2中选择j个
int main(){
    int n;
    scanf("%d",&n);
    for(int k=1;k<=n;k++){
        scanf("%s%s%s",s1+1,s2+1,str+1);
        int len1=strlen(s1+1);
        int len2=strlen(s2+1);
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=len1;i++){
            if(s1[i]==str[i]){
                dp[i][0]=1;
            }
        }
        for(int i=0;i<=len2;i++){
            if(s2[i]==str[i]){
                dp[0][i]=1;
            }
        }
        for(int i=0;i<=len1;i++){
            for(int j=1;j<=len2;j++){
                if(dp[i][j-1]&&s2[j]==str[i+j]){
                    dp[i][j]=1;
                }
                if(dp[i-1][j]&&s1[i]==str[i+j]){
                    dp[i][j]=1;
                }
            }
        }
        printf("Data set %d: ",k);
        if(dp[len1][len2]){
            printf("yes\n");
        }
        else{
            printf("no\n");
        }
    }
    return 0;
}

八、HDU 1502:Regular Words

dp[i][j][k]=dp[i-1][[j][k]+dp[i][j-1][k]+dp[i][j][k-1]
(i>=j&&j>=k)

然后一顿操作猛如虎,结果wa

原来是会暴long long,也就是大整数,于是立马用Java的BigInteger

import java.math.BigDecimal;
import java.math.BigInteger;
import java.util.*;
public class Main {
    public static void main(String args[]) {
        BigInteger dp[][][]=new BigInteger[100][100][100];
        /**
         * 初始化dp数组
         */
        for(int i=0;i<=60;i++) {
            for(int j=0;j<=60;j++) {
                for(int k=0;k<=60;k++) {
                    dp[i][j][k]=BigInteger.ZERO;
                }
            }
        }
        dp[0][0][0]=BigInteger.ONE;
        
        for(int i=0;i<=60;i++) {
            for(int j=0;j<=60;j++) {
                for(int k=0;k<=60;k++) {
                    int x=i-1;
                    int y=j-1;
                    int z=k-1;
                    if(x>=0&&j>=0&&k>=0&&x>=j&&x>=k&&j>=k) {
                        dp[i][j][k]=dp[i][j][k].add(dp[i-1][j][k]);
                    }
                    if(i>=0&&y>=0&&k>=0&&i>=y&&x>=k&&y>=k) {
                        dp[i][j][k]=dp[i][j][k].add(dp[i][j-1][k]);
                    }
                    if(i>=0&&j>=0&&z>=0&&i>=j&&x>=z&&y>=z) {
                        dp[i][j][k]=dp[i][j][k].add(dp[i][j][k-1]);
                    }
                }
            }
        }
        
        Scanner scan=new Scanner(System.in);
        while(scan.hasNext()) {
            int n=scan.nextInt();
            System.out.println(dp[n][n][n]);
            System.out.println();
        }
    }
    
}

day 7:

一、HDU 1506:Largest Rectangle in a Histogram

题意:
n个矩阵,求高相同的情况最大矩形面积

题解:
说实话不是很懂。。。。

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=100000;
ll a[maxn+10];
ll l[maxn+10],r[maxn+10];
int main(){
    int n;
    while(~scanf("%d",&n)){
        if(n==0)
            break;
        for(int i=1;i<=n;i++){
            scanf("%lld",&a[i]);
        }
        a[0]=a[n+1]=-1;
        for(int i=1;i<=n;i++){
            l[i]=r[i]=i;
        }
        //dp,根据之前得到的l,r再拓展,以提高效率
        for(int i=1;i<=n;i++){
            //h[j]>=h[i]的最小j
            while(a[l[i]-1]>=a[i]){
                l[i]=l[l[i]-1];
            }
        }
        for(int i=n;i>=1;i--){
            //h[j]<=h[i]的最大j
            while(a[r[i]+1]>=a[i]){
                r[i]=r[r[i]+1];
            }
        }
        ll ans=0;
        for(int i=1;i<=n;i++){
            ans=max(ans,a[i]*(r[i]-l[i]+1));
        }
        cout<<ans<<endl;
    }
    return 0;
}

二、HDU 1505:City Game

题意:
1506的加强版,这次是求二维的面积

题解:

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=100000;
int h[maxn+10];
int l[maxn+10],r[maxn+10];
int main(){
    int t,n,m;
    char s[10];
    scanf("%d",&t);
    while(t--){
        memset(h,0,sizeof(h));
        scanf("%d%d",&n,&m);
        h[0]=h[m+1]=-1;
        int ans=0;
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%s",s);
                if(s[0]=='F'){
                    h[j]++;
                }
                else{
                    h[j]=0;
                }
            }
            for(int j=1;j<=m;j++){
                l[j]=r[j]=j;
            }
            for(int j=1;j<=m;j++){
                while(h[l[j]-1]>=h[j]){
                    l[j]=l[l[j]-1];
                }
            }
            for(int j=m;j>=1;j--){
                while(h[r[j]+1]>=h[j]){
                    r[j]=r[r[j]+1];
                }
            }
            for(int j=1;j<=m;j++){
                ans=max(ans,h[j]*(r[j]-l[j]+1));
            }
        }
        cout<<ans*3<<endl;
    }
    return 0;
}

三、HDU 1513:Palindrome

题意:
给定一字符串,并且还告诉你它的长度,求加几个字符能是其成为回文字符串

题解:
开始真的一脸懵逼,但是构造一个 字符串由给定字符串的翻转得来
求:最长公共子序列长度即为答案

但是由于字符串长度比较大,必须使用滚动数组,具体见代码

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=100000;
int n;
char s1[5010];
char s2[5010];
int dp[2][5010];       //由于5010*5010比较大,所以需要用到滚动数组
int main(){
    while(~scanf("%d",&n)){
        scanf("%s",s1);
        int t=0;
        for(int i=n-1;i>=0;i--){
            s2[t++]=s1[i];
        }
        //cout<
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(s1[i-1]==s2[j-1]){
                    dp[i&1][j]=dp[(i-1)&1][j-1]+1;
                }
                else{
                    dp[i&1][j]=max(dp[(i-1)&1][j],dp[i&1][j-1]);
                }
            }
        }
        printf("%d\n",n-dp[n&1][n]);
    }
    return 0;
}

day 8:

一、HDU 1559:最大子矩阵

题意:
mn的矩形,从中选择xy的子矩阵,问最大子矩阵为多少

题解:
很明显的dp问题,但是感觉处理起来还是有点技巧的,具体见代码

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=200000;
int t;
int m,n,x,y;
int dp[1010][1010];
int main(){
    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d%d",&m,&n,&x,&y);
        int ans=0;
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=m;i++){
            for(int j=1;j<=n;j++){
                scanf("%d",&dp[i][j]);
                dp[i][j]+=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1];
                if(i>=x&&j>=y){
                    ans=max(ans,dp[i][j]-dp[i][j-y]-dp[i-x][j]+dp[i-x][j-y]);
                }
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}

二、HDU 1574:RP问题

题意:
n个事件,每个事件对应有RP变化值a、RP门槛值b和获益值c
当RP变化值a为正,获益值c必定为负,只有你当前的RP值小于等于RP门槛值b的时候,此事件才有可能发生,当此事件发生时,你的RP值将增加|a|,获益值将减少|c|。反之,当RP变化值a为负,获益值c必定为正,只有你当前的RP值大于等于RP门槛值b的时候,此事件才有可能发生,当此事件发生时,你的RP值将减少|a|,获益值将增加|c|。

求依次经过n件事件后,最终的收益值

题解:
0-1背包的变形,但是还是有点难想

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=20000;
//0-1背包的变形
//dp[j+a[i]]=max(dp[j+a[i]],dp[j]+c[i]);
int dp[maxn+10];
int a[maxn+10],b[maxn+10],c[maxn+10];
int main(){
    int t,n;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        memset(dp,-inf,sizeof(dp));
        for(int i=1;i<=n;i++){
            scanf("%d%d%d",&a[i],&b[i],&c[i]);
        }
        for(int i=1;i<=n;i++){
            dp[10000]=0;
            if(a[i]>0){
                for(int j=b[i]+10000;j>=0;j--){
                    dp[j+a[i]]=max(dp[j+a[i]],dp[j]+c[i]);
                }
            }
            else{
                for(int j=b[i]+10000;j<=20000;j++){
                    dp[j+a[i]]=max(dp[j+a[i]],dp[j]+c[i]);
                }
            }
        }
        int ans=0;
        for(int i=0;i<=20000;i++){
            ans=max(ans,dp[i]);
        }
        printf("%d\n",ans);
    }
    return 0;
}

三、HDU 1561:The more, The Better

树形dp,先放在这

四、HDU 1565:方格取数(1)

状压dp

五、HDU 1693:Eat the Trees

插头dp

六、HDU 1712:ACboy needs your help

题意:
有n门课程,m个单位时间
给出一矩阵,a[i][j]表示:i门课程花j时间得到的价值

求在m时间内,能得到的最大价值

题解:

总觉得这题是背包问题,但是就是 没有什么思路.。。
于是又一次看了看背包九讲
这不就是典型的分组背包

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=20000;
int n,m;
int a[110][110];
int dp[110];
//分组背包问题
/*
for 所有的组k
    for v=V..0
        for 所有的i属于组k
            f[v]=max{f[v],f[v-c[i]]+w[i]}
*/
int main(){
    while(~scanf("%d%d",&n,&m)){
        if(n==0&&m==0)
            break;
        memset(a,0,sizeof(a));
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&a[i][j]);
            }
        }
        for(int i=1;i<=n;i++){         //第i组
            for(int j=m;j>=0;j--){     //背包容量
                for(int k=1;k<=j;k++){
                    dp[j]=max(dp[j],dp[j-k]+a[i][k]);
                }
            }
        }
        printf("%d\n",dp[m]);
    }
    return 0;
}

七、HDU 1723:Distribute Message

题意:
n个人排成一行,现在要求将一物品从第一个人传到第n个人
其中每个人最多可以跳多m个位置
问有多少种方案
(题意自己稍微修改了下文字,哈哈哈,觉得这样好理解些)

题意:
比较简单的dp问题

dp[n]=dp[n-1]+dp[n-2]+…+dp[n-m];

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define lowbit(x) x&(-x)
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
const double pi=acos(-1.0);
const int maxn=20000;
int n,m;
int cnt;
int dp[100];
//dp[n]=dp[n-1]+dp[n-2]+...+dp[n-m];
int main(){
    while(~scanf("%d%d",&n,&m)){
        if(n==0&&m==0)
            break;
        memset(dp,0,sizeof(dp));
        dp[1]=1;
        for(int i=2;i<=n;i++){
            for(int j=1;j<=m;j++){
                if(i-j>0)               //这个千万不能省...
                    dp[i]+=dp[i-j];
            }
        }
        printf("%d\n",dp[n]);
    }
    return 0;
}

你可能感兴趣的:(DP)