【POJ 3666】Making the Grade(简单DP)

2024-09-07 23:32
文章标签 简单 dp poj grade making 3666

本文主要是介绍【POJ 3666】Making the Grade(简单DP),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

首先可以明确一个方面,那就是如果将X改成Y,那么Y肯定是这N个数中的某一个(为什么仔细想想)

之后可以得到一个状态转移那就是dp[i][j]代表已经考虑了i位的情况下,结尾为j的最小更改数。

状态转移为dp[i][j] = min(dp[i-1][k] + abs(a[i] - b[j])) 这样的话可以写出一个初步的代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int  INF = 1 << 30;
const int maxn = 2005;
int arr[maxn];
int Hash[maxn],cnt;
int dp[maxn][maxn];
int n;
void Hash_init(){sort(Hash,Hash + cnt);cnt = unique(Hash,Hash + cnt) - Hash;//printf("%d\n",cnt);
}
void solve(){dp[0][0] = 0;for(int i = 1; i < n; i++){for(int j = 0; j < cnt; j++){dp[i][j] = INF;for(int k = 0; k <= j; k++){dp[i][j] = min(dp[i][j],dp[i - 1][k] + abs(Hash[j] - arr[i]));}}}int ans = dp[n - 1][0];for(int i = 1; i < cnt; i++)ans = min(ans,dp[n - 1][i]);printf("%d\n",ans);
}
int main(){while(scanf("%d",&n) != EOF){cnt = 0;for(int i = 0; i < n; i++){scanf("%d",&arr[i]);Hash[cnt++] = arr[i];}Hash_init();solve();}return 0;
}


继续看,这里3个for循环对于上千的数据肯定是不行的,既然这样我们优化,首先我们可以知道每次取min(dp[i - 1][t])就可以,这里t <= j,这里我们用一个数组pre[]来代表这个结果,这样直接去了一层循环,pre[t]每次循环需要更新


#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int  INF = 1 << 30;
const int maxn = 2005;
int arr[maxn],pre[maxn];
int Hash[maxn],cnt;
int dp[maxn][maxn];
int n;
void Hash_init(){sort(Hash,Hash + cnt);cnt = unique(Hash,Hash + cnt) - Hash;
}
void solve(){dp[0][0] = 1;for(int i = 1; i < n; i++){for(int j = 0; j < cnt; j++){dp[i][j] = pre[j] + abs(Hash[j] - arr[i]);if(!j)pre[j] = dp[i][j];elsepre[j] = min(pre[j - 1],dp[i][j]);}}int ans = dp[n - 1][0];for(int i = 1; i < cnt; i++)ans = min(ans,dp[n - 1][i]);printf("%d\n",ans);
}
int main(){while(scanf("%d",&n) != EOF){cnt = 0;for(int i = 0; i < n; i++){scanf("%d",&arr[i]);Hash[cnt++] = arr[i];}Hash_init();solve();}return 0;
}


这样已经可以AC了,但是我们还可以利用滚动数组去优化空间

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 2005;
int arr[maxn],pre[maxn];
int Hash[maxn],cnt;
int dp[2][maxn];
int n;
void Hash_init(){sort(Hash,Hash + cnt);cnt = unique(Hash,Hash + cnt) - Hash;
}
void solve(){dp[0][0] = 1;int k = 0;for(int i = 1; i < n; i++){for(int j = 0; j < cnt; j++){dp[k][j] = pre[j] + abs(Hash[j] - arr[i]);if(!j)pre[j] = dp[k][j];elsepre[j] = min(pre[j - 1],dp[k][j]);}k ^= 1;}k ^= 1;int ans = dp[k][0];for(int i = 1; i < cnt; i++)ans = min(ans,dp[k][i]);printf("%d\n",ans);
}
int main(){while(scanf("%d",&n) != EOF){cnt = 0;for(int i = 0; i < n; i++){scanf("%d",&arr[i]);Hash[cnt++] = arr[i];}Hash_init();solve();}return 0;
}
恩,就是这样~

这篇关于【POJ 3666】Making the Grade(简单DP)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1146478

相关文章

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu2289(简单二分)

虽说是简单二分,但是我还是wa死了  题意:已知圆台的体积,求高度 首先要知道圆台体积怎么求:设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ

hdu4826(三维DP)

这是一个百度之星的资格赛第四题 题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1004&cid=500 题意:从左上角的点到右上角的点,每个点只能走一遍,走的方向有三个:向上,向下,向右,求最大值。 咋一看像搜索题,先暴搜,TLE,然后剪枝,还是TLE.然后我就改方法,用DP来做,这题和普通dp相比,多个个向上

hdu1011(背包树形DP)

没有完全理解这题, m个人,攻打一个map,map的入口是1,在攻打某个结点之前要先攻打其他一个结点 dp[i][j]表示m个人攻打以第i个结点为根节点的子树得到的最优解 状态转移dp[i][ j ] = max(dp[i][j], dp[i][k]+dp[t][j-k]),其中t是i结点的子节点 代码如下: #include<iostream>#include<algorithm

hdu4865(概率DP)

题意:已知前一天和今天的天气概率,某天的天气概率和叶子的潮湿程度的概率,n天叶子的湿度,求n天最有可能的天气情况。 思路:概率DP,dp[i][j]表示第i天天气为j的概率,状态转移如下:dp[i][j] = max(dp[i][j, dp[i-1][k]*table2[k][j]*table1[j][col] )  代码如下: #include <stdio.h>#include

usaco 1.3 Prime Cryptarithm(简单哈希表暴搜剪枝)

思路: 1. 用一个 hash[ ] 数组存放输入的数字,令 hash[ tmp ]=1 。 2. 一个自定义函数 check( ) ,检查各位是否为输入的数字。 3. 暴搜。第一行数从 100到999,第二行数从 10到99。 4. 剪枝。 代码: /*ID: who jayLANG: C++TASK: crypt1*/#include<stdio.h>bool h

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

poj 1511 Invitation Cards(spfa最短路)

题意是给你点与点之间的距离,求来回到点1的最短路中的边权和。 因为边很大,不能用原来的dijkstra什么的,所以用spfa来做。并且注意要用long long int 来存储。 稍微改了一下学长的模板。 stack stl 实现代码: #include<stdio.h>#include<stack>using namespace std;const int M