ACM总结 - 动态规划 - 第六周

好题儿

贪心与动态规划相结合(C题)

根据题目的意思,老鼠的体重要递增,而老鼠的速度要递减。刚开始以为是对两个数组求最大不下降子序列。通过对两个数组(体重、速度)求LCS后,找出相同的下标个数,就是所求的那一组数据的下标和,再输出这些下标不就行了吗。但是后来发现了一个一个问题:这些数据并不是按照顺序输入的,也就是说,要先进行排序才可以。
这时就想到了用贪心算法:以weight为第一元素从小到大排序,以speed为第二元素从大到小。此时的状态转移方程为:

d p [ i ] = 1 dp[i] = 1 dp[i]=1;//因为只有一个数据时要有值,所以dp初始值设置为1

d p [ i ] = max ⁡ ( d p [ i ] , d p [ j ] + 1 ) dp[i] = \max(dp[i], dp[j] + 1) dp[i]=max(dp[i],dp[j]+1);

( j < i , p [ j ] . w e i g h t > p [ i ] . w e i g h t , p [ j ] . s p e e d < p [ i ] . s p e e d ) (j < i, p[j].weight > p[i].weight, p[j].speed < p[i].speed) (j<i,p[j].weight>p[i].weight,p[j].speed<p[i].speed);//对i进行降序循环,对j进行初始值为i的升序循环,找出满足体重升速度减的元素。
代码如下:

#include<iostream>
#include<cstdio>
#include<algorithm>

using namespace std;

struct node{
    int weight;
    int speed;
    int num;
}p[10010];

int dp[10010];

bool cmp(node A,node B){
    if(A.weight != B.weight) return A.weight < B.weight;
    else return A.speed > B.speed;
}

int main(){
    int k = 1, x, y;
    while(scanf("%d%d", &x, &y) != EOF){
        p[k].weight = x;
        p[k].speed = y;
        p[k].num = k;
        k++;
    }
    sort(p + 1, p + k + 1, cmp);
    int ans = 0;
    for(int i = k; i >= 1; i--){
            dp[i] = 1;
        for(int j = i; j <= k; j++){
            if(p[j].weight > p[i].weight && p[j].speed < p[i].speed){
                dp[i] = max(dp[i], dp[j]+1);
            }
        }
        ans = max(ans, dp[i]);
    }
    printf("%d\n", ans);
    for(int i = 1; i <= k; i++){
        if(dp[i] == ans){
            printf("%d\n", p[i].num);
            ans--;
        }
        if(ans == 0) break;
    }
	return 0;
}

G两种方式

这一题和常见的求吃金币的不一样,这一题是求走到终点的方案数。这一题可以反过来运用动态规划,有点类似《算法笔记》里的数塔问题,也就是记忆化搜索,就是从终点开始进行计算,设置终点值为1,搜索过程中记录每一点的值。此外,还可以通过正常的方式,从起始点开始,每走一格就枚举出能到达终点的所有情况。
注意,每次运行需要对p数组清零,也就是memset一下。注意memset最好只是对初始化为0或-1的数组运用,对初始化为其他数的数组最好用fill。
这里运用第二种方法,也就是比较常用的“正向动态规划”。
代码如下:

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

int dp[111][111], a[111][111];

int main(){
    int c, i, j, n, m, x, y, di, dj, t;
    scanf("%d", &c);
    while(c--){
        scanf("%d %d", &n, &m);
        for(i = 1; i <= n; i++) for(j = 1; j <= m; j++) scanf("%d",&a[i][j]);
        memset(dp, 0, sizeof(dp));
        dp[1][1] = 1;
        for(i = 1; i <= n; i++){
            for(j = 1; j <= m; j++){
                if(dp[i][j]){
                    t = a[i][j];
                    for(x = 0; x <= t; x++){
                        for(y = 0; y <= t; y++){
                            if(x + y == 0) continue; if(x + y > t) break;
                            di = x + i; dj = y + j;
                            if(di <= n && dj <= m) dp[di][dj] += dp[i][j];
                            dp[di][dj] %= 10000;
                        }
                    }
                }
            }
        }
        printf("%d\n", dp[n][m] % 10000);
    }    
    return 0;
}

此外我在网上找到了某博主写的第一种方法的代码:


#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
#define N 105
const int mod=10000;
int n,m;
int a[N][N];
int mark[N][N];
int dfs(int x,int y)
{
    if(mark[x][y]!=-1)
        return mark[x][y];
    int t,x1,y1,tmp;
    tmp=0;
    t=a[x][y];
    for(x1=0;x1<=t;x1++)
    {
        if(x+x1>=n)
            break;
        for(y1=0;x1+y1<=t;y1++)
        {
            if(y+y1>=m)
                break;
            if(x1+y1!=0)          //不能原地走
                tmp+=dfs(x+x1,y+y1);
        }
    }
    tmp%=mod;
    return mark[x][y]=tmp;
}
int main()
{
    int i,j,T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        for(i=0;i<n;i++)
        {
            for(j=0;j<m;j++)
            {
                scanf("%d",&a[i][j]);
                mark[i][j]=-1;
            }
        }
        mark[n-1][m-1]=1;
        printf("%d\n",dfs(0,0));
    }
    return 0;
}

S题的奇妙思路

这一题提交了好多次才通过,算是比较难得了。思路不亚于前两个题。
这题也就是说疲劳度与左右手重量差有关,也就是说让小红豆最轻松,我们就要让他所有次数左右手的重量差平方和最小。
那么我们可以按照重量对物体进行排序,然后选取重量相邻的两个物体来进行比较。
设f[i][j]为前i件物品组成k对的话,所消耗的体力最小;
这时分两种情况:含有第i件物品和不含有第i件物品(第i件物品是不是含在第j对里)
1.含有i件物品,则有 f [ i ] [ j ] = f [ i − 2 ] [ j − 1 ] + ( v a l [ i ] − v a l [ i − 1 ] ) ∗ ( v a l [ i ] − v a l [ i − 1 ] ) f[i][j]=f[i-2][j-1]+(val[i]-val[i-1])*(val[i]-val[i-1]) f[i][j]=f[i2][j1]+(val[i]val[i1])(val[i]val[i1])
2.不含第i件物品,则有 f [ i ] [ j ] = f [ i − 1 ] [ j ] f[i][j]=f[i-1][j] f[i][j]=f[i1][j]
所以动态转移方程为: f [ i ] [ j ] = min ⁡ ( f [ i − 2 ] [ j − 1 ] + ( v a l [ i ] − v a l [ i − 1 ] ) ∗ ( v a l [ i ] − v a l [ i − 1 ] ) , f [ i ] [ j ] = f [ i − 1 ] [ j ] ) f[i][j]=\min(f[i-2][j-1]+(val[i]-val[i-1])*(val[i]-val[i-1]), f[i][j]=f[i-1][j]) f[i][j]=min(f[i2][j1]+(val[i]val[i1])(val[i]val[i1]),f[i][j]=f[i1][j])
代码如下:

#include <iostream>
#include <cstdio>
#include <algorithm>

using namespace std;

#define INF 0xfffffff

int list[2001];
int dp[1001][2001];

int main(){
	int n, k, i, j;
	while(scanf("%d%d", &n, &k) != EOF){
		for(i = 1; i <= n; i++) cin >>list[i];
		sort(list + 1, list + n + 1);
		for(i = 1; i <= n; i++) dp[0][i] = 0;
		for(i = 1; i <= k; i++){
			for(j = 2 * i; j <= n; j++){
				if(j > 2 * i) dp[i][j] = dp[i][j - 1];
				else dp[i][j] = INF;
				if(dp[i][j] > dp[i - 1][j - 2] + (list[j] - list[j - 1]) * (list[j] - list[j - 1])) dp[i][j] = dp[i - 1][j - 2] + (list[j] - list[j - 1]) * (list[j] - list[j - 1]);
			}
		}
		cout << dp[k][n] << endl;
	}
	return 0;
}

总结

这一周以刷题为主,题目还是比较困难的。不过其中有规律性的东西,比如上面三个,贪心与动态规划结合之类的,只要掌握了,基本上就OK了。难题儿还是比较多的,其实题儿练多了就不难了,期待下一周的背包。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值