好久没做DP了,虽说是基础DP1吧,第一题就把我海虐了
感觉DP真是看脑洞,最近脑子各种混乱睡眠也不足,还是得要好好理解的
HDU 1024
题意:最大m段区间和啊
题解:这题正常的方法都会MLE 啊TLE啊
dp[i][j]这个数组表示前j个元素中取i段的最大值,并且取了a[j]
自然是开不下二维的
于是要用一维滚动数组(这方面我真是不行,滚动数组很久以前也学过,一碰到难题就想不出)
用dp[j]表示取了a[j]情况下的当前一层(取多少组区间)的最大值,然后要开个num[j]数组,记录1-j情况下(上一层)取到的最大值
先是循环1-m,表示当前是取了多少个区间,然后元素从i-n遍历(因为前面取了i-1层,至少有i-1个元素,所以j从i开始考察)
第i层的dp[j]是第i层的dp[j-1]+a[j],(把a[j-1]和a[j]看在一起),或者是num[j-1]+a[j]。
要开个maxn记录这一层dp[1-j]的最大值,然后赋值给num[j],由于到达i+1层之后num[i]就不成立了所以一直开始maxn赋值是INF,然后maxn=max(maxn,dp[j]),在j+1的时候赋值给num[j]
另外我开LL 会TLE,然后比较虚,看discuss里说int也能过,果然
int a[MAX]; int dp[MAX]; int num[MAX]; int main(){ int n,m; while(~scanf("%d%d",&m,&n)){ for(int i=1;i<=n;i++){ scanf("%d",&a[i]); } mem0(dp); mem0(num); int maxn; for(int i=1;i<=m;i++){ maxn=-INF; for(int j=i;j<=n;j++){ dp[j]=max(dp[j-1],num[j-1])+a[j]; num[j-1]=maxn; maxn=max(dp[j],maxn); } } printf("%d\n",maxn); } return 0; }
poj 1015
在遥远的国家佛罗布尼亚,嫌犯是否有罪,须由陪审团决定。陪审团是
由法官从公众中挑选的。先随机挑选n 个人作为陪审团的候选人,然后
再从这n 个人中选m 人组成陪审团。选m 人的办法是:控方和辩方会
根据对候选人的喜欢程度,给所有候选人打分,分值从0 到20。为了
公平起见,法官选出陪审团的原则是:选出的m 个人,必须满足辩方
总分和控方总分的差的绝对值最小。如果有多种选择方案的辩方总分
和控方总分的之差的绝对值相同,那么选辩控双方总分之和最大的
方案即可。最终选出的方案称为陪审团方案。
为叙述问题方便,现将任一选择方案中,辩方总分和控方总分
之差简称为“辩控差”,辩方总分和控方总分之和称为“辩控和”。
第i 个候选人的辩方总分和控方总分之差记为V(i),辩方总分和控
方总分之和记为S(i)。现用f(j, k)表示,取j 个候选人,使其辩
控差为k 的所有方案中,辩控和最大的那个方案(该方案称为“方
案f(j, k)”)的辩控和。并且,我们还规定,如果没法选j 个人,
使其辩控差为k,那么f(j, k)的值就为-1,也称方案f(j, k)不可行。
本题是要求选出m 个人,那么,如果对k 的所有可能的取值,求
出了所有的f(m, k) (-20×m≤ k ≤ 20×m),那么陪审团方案
自然就很容易找到了。
问题的关键是建立递推关系。需要从哪些已知条件出发,
才能求出f(j, k)呢?显然,方案f(j, k)是由某个可行的方案f(j-1, x)
( -20×m ≤ x ≤ 20×m)演化而来的。可行方案f(j-1, x)能演化成
方案f(j, k)的必要条件是:存在某个候选人i,i 在方案f(j-1, x)中
没有被选上,且x+V(i) = k。在所有满足该必要条件的f(j-1, x)中,
选出 f(j-1, x) + S(i) 的值最大的那个,那么方案f(j-1, x)再加上候选人i,
就演变成了方案 f(j, k)。这中间需要将一个方案都选了哪些人都记录下来。
不妨将方案f(j, k)中最后选的那个候选人的编号,记在二维数组的
元素path[j][k]中。那么方案f(j, k)的倒数第二个人选的编号,
就是path[j-1][k-V[path[j][k]]。假定最后算出了解方案的辩控差是k,
那么从path[m][k]出发,就能顺藤摸瓜一步步求出所有被选中的候选人。
初始条件,只能确定f(0, 0) = 0。由此出发,一步步自底向上递推,
就能求出所有的可行方案f(m, k)( -20×m ≤ k ≤ 20×m)。实际解题
的时候,会用一个二维数组f 来存放f(j, k)的值。而且,由于题目中辩
控差的值k 可以为负数,而程序中数租下标不能为负数,所以,在程序中
不妨将辩控差的值都加上400,以免下标为负数导致出错,即题目描述中,
如果辩控差为0,则在程序中辩控差为400。
这题好难,自己瞎搞了好久都跪了,然后就看kuangbin大神的题解,确实厉害啊
#include <map> #include <set> #include <stack> #include <queue> #include <cmath> #include <string> #include <vector> #include <cstdio> #include <cctype> #include <cstring> #include <sstream> #include <cstdlib> #include <iostream> #include <algorithm> using namespace std; #define MAX 205 #define MAXN 10005 #define lson l,m,rt<<1 #define rson m+1,r,rt<<1|1 #define lrt rt<<1 #define rrt rt<<1|1 #define mid int m=(r+l)>>1 #define LL long long #define ull unsigned long long #define mem0(x) memset(x,0,sizeof(x)) #define mem1(x) memset(x,-1,sizeof(x)) #define meminf(x) memset(x,INF,sizeof(x)) #define lowbit(x) (x&-x) const int mod = 1000000007; const int prime = 999983; const int INF = 0x3f3f3f3f; const int INFF = 1e9; const double pi = 3.141592653589793; const double inf = 1e18; const double eps = 1e-10; int path[25][805];//记录j,k状态下选择的是第几个人 int d[MAX],p[MAX]; int dp[25][805];//第二维是辩控差,数组记录辩控和 int num[MAX]; int main(){ int n,m; int kase=0; while(scanf("%d%d",&n,&m)&&n){ kase++; for(int i=1;i<=n;i++) scanf("%d%d",&p[i],&d[i]); mem1(path); mem1(dp); int mm=m*20;//辩控差0转化为mm dp[0][mm]=0; for(int j=1;j<=m;j++){ for(int k=0;k<=2*mm;k++){ if(dp[j-1][k]>=0){//情况可行 for(int i=1;i<=n;i++){ if(dp[j-1][k]+p[i]+d[i]>dp[j][k+p[i]-d[i]]){//前一种情况加上第i个人,大于此时j个人的情况,类似于背包 int t1=j-1,t2=k; while(t1>0&&path[t1][t2]!=i){//检查前j-1个人中有没有i t2-=p[path[t1][t2]]-d[path[t1][t2]]; t1--; } if(t1==0){//没有出现过i dp[j][k+p[i]-d[i]]=dp[j-1][k]+p[i]+d[i]; path[j][k+p[i]-d[i]]=i; } } } } } } int j=mm; int i=0; int tmp; while(dp[m][j-i]<0&&dp[m][j+i]<0) i++;//看辩控差最小是多少 if(dp[m][j+i]>dp[m][j-i]) tmp=j+i;//取辩控和大的 else tmp=j-i; int tot=0; for(int i=path[m][tmp],t=tmp,j=m;i!=-1&&j>0;i=path[j][t]){ j--; t-=p[i]-d[i]; num[tot++]=i;//顺着找回去 } sort(num,num+tot); printf("Jury #%d\n",kase); printf("Best jury has value %d for prosecution and value %d for defence:\n",(dp[m][tmp]+tmp-mm)/2,(dp[m][tmp]-tmp+mm)/2); for(int i=0;i<tot;i++) printf(" %d",num[i]); printf("\n\n"); } return 0; }
HDU 1160
题意:给你一些老鼠,有质量和速度两个参数,让你找到最多的老鼠的排列,让他们质量递增,速度递减
题解:dp[i]表示选择第i只老鼠作为最后一只的时候最多能选多少只,先对老鼠的质量升序排序,速度降序,然后dp[i]=max(dp[k])k<i;
就这样遍历,o(n^2)的复杂度可以过,然后开个vector记录路径,当前dp[i]最优解的情况下的序号
#include <map> #include <set> #include <stack> #include <queue> #include <cmath> #include <string> #include <vector> #include <cstdio> #include <cctype> #include <cstring> #include <sstream> #include <cstdlib> #include <iostream> #include <algorithm> using namespace std; #define MAX 1005 #define MAXN 10005 #define lson l,m,rt<<1 #define rson m+1,r,rt<<1|1 #define lrt rt<<1 #define rrt rt<<1|1 #define mid int m=(r+l)>>1 #define LL long long #define ull unsigned long long #define mem0(x) memset(x,0,sizeof(x)) #define mem1(x) memset(x,-1,sizeof(x)) #define meminf(x) memset(x,INF,sizeof(x)) #define lowbit(x) (x&-x) const int mod = 1000000007; const int prime = 999983; const int INF = 0x3f3f3f3f; const int INFF = 1e9; const double pi = 3.141592653589793; const double inf = 1e18; const double eps = 1e-10; struct mice{ int w,s,id; bool operator < (const mice &a) const{ if(w==a.w) return s>a.s; return w<a.w; } }p[MAX]; vector<int> v[MAX]; int dp[MAX]; int main(){ int xcount=1; int a,b; while(~scanf("%d%d",&a,&b)){ p[xcount]=(mice){a,b,xcount}; xcount++; } sort(p+1,p+xcount); int ans=0; int cnt=0; for(int i=1;i<xcount;i++){ dp[i]=1; int maxn=0; int flag=-1; for(int j=1;j<i;j++){ if(p[i].w>p[j].w&&p[i].s<p[j].s){ if(maxn<dp[j]){ maxn=dp[j]; flag=j; } } } if(flag!=-1){ dp[i]+=maxn; for(int k=0;k<v[flag].size();k++){ v[i].push_back(v[flag][k]); } if(ans<dp[i]){ ans=dp[i]; cnt=i; } } v[i].push_back(p[i].id); } printf("%d\n",ans); for(int i=0;i<v[cnt].size();i++){ printf("%d\n",v[cnt][i]); } return 0; }