找啊找啊找GF (多维动态规划)

【问题描述】

"找啊找啊找GF,找到一个好GF,吃顿饭啊拉拉手,你是我的好GF.再见."

"诶,别再见啊..."

七夕...七夕...七夕这个日子,对于sqybi这种单身的菜鸟来说是多么的痛苦...虽然他听着这首叫做"找啊找啊找GF"的歌,他还是很痛苦.为了避免这种痛苦,sqybi决定要给自己找点事情干.他去找到了七夕模拟赛的负责人zmc MM,让她给自己一个出题的任务.经过几天的死缠烂打,zmc MM终于同意了.

但是,拿到这个任务的sqybi发现,原来出题比单身更让人感到无聊-_-....所以,他决定了,要在出题的同时去办另一件能够使自己不无聊的事情--给自己找GF.

sqybi现在看中了n个MM,我们不妨把她们编号1到n.请MM吃饭是要花钱的,我们假设请i号MM吃饭要花rmb[i]块大洋.而希望骗MM当自己GF是要费人品的,我们假设请第i号MM吃饭试图让她当自己GF的行为(不妨称作泡该MM)要耗费rp[i]的人品.而对于每一个MM来说,sqybi都有一个对应的搞定她的时间,对于第i个MM来说叫做time[i]. sqybi保证自己有足够的魅力用time[i]的时间搞定第i个MM^_^.

sqybi希望搞到尽量多的MM当自己的GF,这点是毋庸置疑的.但他不希望为此花费太多的时间(毕竟七夕赛的题目还没出),所以他希望在保证搞到MM数量最多的情况下花费的总时间最少.

sqybi现在有m块大洋,他也通过一段时间的努力攒到了r的人品(这次为模拟赛出题也攒rp哦~~).他凭借这些大洋和人品可以泡到一些MM.他想知道,自己泡到最多的MM花费的最少时间是多少.

注意sqybi在一个时刻只能去泡一个MM--如果同时泡两个或以上的MM的话,她们会打起来的...

【输入文件】

输入的第一行是n,表示sqybi看中的MM数量.

接下来有n行,依次表示编号为1, 2, 3, ..., n的一个MM的信息.每行表示一个MM的信息,有三个整数:rmb, rp和time.

最后一行有两个整数,分别为m和r.

【输出文件】

你只需要输出一行,其中有一个整数,表示sqybi在保证MM数量的情况下花费的最少总时间是多少.

【输入样例】

4

1 2 5

2 1 6

2 2 2

2 2 3

5 5

【输出样例】

13

【数据规模】

对于20%数据,1<=n<=10;

对于100%数据,1<=rmb<=100,1<=rp<=100,1<=time<=1000;

对于100%数据,1<=m<=100,1<=r<=100,1<=n<=100.

【提交链接】

http://www.rqnoj.cn/Submit.asp

【问题分析】

初看问题觉得条件太多,理不出头绪来,所以要将问题简化,看能否找出熟悉的模型来,如果我们只考虑钱够不够,或只考虑RP够不够。并且不考虑花费的时间。这样原问题可以简化成下面的问题:

在给定M元RMB(或R单位RP,RP该用什么单位呢?汗。。。)的前题下,去泡足够多的MM,很显然这个问题就是典型的0/1背包问题了。

可以把泡MM用的RMB(或RP看做重量),泡到MM的个数看做价值,给定的M(或R)就是背包的载重。求解这个问题很轻松喽。

但是,这个问题既要考虑RMB有要考虑RP怎么办呢?

解决这个问题很容易啊,要是你有足够的RMB去泡第i个MM而RP不够就泡不成了,要是RP够就可以。也就是在原来问题的基础上在状态加一维。

那要是在考虑上时间最小怎么办呢?

这个也很好说,在求解过程中如果花X元RMP,Y单位RP可以到Z个MM,那么在泡第i个MM时,发现可以用X-rmb[i]元,Y-rp[i]单位RP泡到的MM数加上这个MM(也就是+1)比原来Z多,就替换它(因为你的原则是尽量多的泡MM),如果和Z一样多,这是就要考虑原来花的时间多呢,还是现在花的时间多。要是原来的多,就把时间替换成现在用的时间(因为你既然可以泡到相同数量的MM当然要省点时间去出题)。

设计一个二维状态opt[j,k]表示正好花j元RMP,k单位RP可以泡到的最多的MM的数量。增加一个辅助的状态ct[k,j]表示正好花j元RMP,k单位RP可以泡到的最多MM的情况下花费的最少的时间。

边界条件 opt[0,0]=1    (按题意应该是0,但为了标记花费是否正好设为1,这样,opt[j,k]>0说明花费正好)

状态转移方程:

opt[j,k]=max{opt[j-rmb[i],k-rp[i]]+1}

(rmb[i]<=j<=m,rp[i]<=k<=r,0<i<=n,opt[j-rmb[i],k-rp[i]]>0)

ct[j,k]:=min{ct[j-rmb[i],k-rp[i]]}+time[i] (opt[j,k]=opt[j-rmb[i],k-rp[i]]+1)

时间复杂度:

阶段数 O(N)*状态数O(MR)*转移代价O(1)=    O(NMR)

注:数据挺小的。

问题拓展:

如果要加入别的条件,比如泡MM还要一定的SP,等也就是说一个价值要不同的条件确定,那么这个问题的状态就需要在加一维,多一个条件就多一维。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

int n, rmb, rp;
int g[101][3];
int dp[101][101][101];
int t[101][101][101];

int main()
{
    while(scanf("%d",&n) != EOF)
    {
        int i, j, k;
        for(i = 1; i <= n; i++)
            scanf("%d%d%d",&g[i][0],&g[i][1],&g[i][2]);
        scanf("%d%d",&rmb,&rp);

        memset(dp,0,sizeof(dp));
        memset(t,0,sizeof(t));
        for(i = 1; i <= n; i++)
        {
            for(j = 0; j <= rmb; j++)
            {
                for(k = 0; k <= rp; k++)
                {
                    if(j < g[i][0] || k < g[i][1])
                    {
                        dp[i][j][k] = dp[i-1][j][k];
                        t[i][j][k] = t[i-1][j][k];
                        continue;
                    }
                   if(dp[i-1][j][k] > dp[i-1][j-g[i][0]][k-g[i][1]] + 1)
                   {
                       dp[i][j][k] = dp[i-1][j][k];
                       t[i][j][k] = t[i-1][j][k];
                   }
                   else if(dp[i-1][j][k] < dp[i-1][j-g[i][0]][k-g[i][1]] + 1)
                   {
                       dp[i][j][k] = dp[i-1][j-g[i][0]][k-g[i][1]] + 1;
                       t[i][j][k] = t[i-1][j-g[i][0]][k-g[i][1]] + g[i][2];
                   }
                   else if(t[i-1][j][k] > t[i-1][j-g[i][0]][k-g[i][1]] + g[i][2])
                   {
                       dp[i][j][k] = dp[i-1][j-g[i][0]][k-g[i][1]] + 1;
                       t[i][j][k] = t[i-1][j-g[i][0]][k-g[i][1]] + g[i][2];
                   }
                   else
                   {
                       dp[i][j][k] = dp[i-1][j][k];
                       t[i][j][k] = t[i-1][j][k];
                   }
                }
            }
        }
        printf("%d\n",t[n][rmb][rp]);
    }
    return 0;
}


优化一下:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

int n, rmb, rp;
int g[101][3];
int dp[101][101];
int t[101][101];

int main()
{
    while(scanf("%d",&n) != EOF)
    {
        int i, j, k;
        for(i = 1; i <= n; i++)
            scanf("%d%d%d",&g[i][0],&g[i][1],&g[i][2]);
        scanf("%d%d",&rmb,&rp);

        memset(dp,0,sizeof(dp));
        memset(t,0,sizeof(t));
        for(i = 1; i <= n; i++)
        {
            for(j = rmb; j >= g[i][0]; j--)
            {
                for(k = rp; k >= g[i][1]; k--)
                {

                   if(dp[j][k] < dp[j-g[i][0]][k-g[i][1]] + 1)
                   {
                       dp[j][k] = dp[j-g[i][0]][k-g[i][1]] + 1;
                       t[j][k] = t[j-g[i][0]][k-g[i][1]] + g[i][2];
                   }

                   if(dp[j][k] == dp[j-g[i][0]][k-g[i][1]] + 1 && t[j][k] > t[j-g[i][0]][k-g[i][1]] + g[i][2])
                   {
                       dp[j][k] = dp[j-g[i][0]][k-g[i][1]] + 1;
                       t[j][k] = t[j-g[i][0]][k-g[i][1]] + g[i][2];
                   }
                }
            }
        }
        printf("%d\n",t[rmb][rp]);
    }
    return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值