标签:
15 5
5 4 3 2 1 1 2 3 4 5 5 4 3 2 1
19
又是从标题就知道这是一道DP类的题。初学DP,这道题搞了有好几个小时Orz...
这次我变化下写作方式,先写下心得,毕竟耗了这莫长时间,没点心得不亏大了吗:-
心得体会:第一点——题目中给出的变量用来限界,状态方程中用其他的变量,防止引起混乱。
第二点——写DP的核心部分时,可以先写伪码,使思路更加清晰。我就是这么做的,伪码一出,分分钟就写完DP了,至于伪码的写法,每个人都可以有自己的一套方法,这个可以在写程序的过程中不断精进。
下面开始分析:设dp(m, i) 表示长度为m 的序列取到第i 个数做尾数的最小子序列和,其中1 ≤ m ≤ i ≤ n。则有递推方程(注意变量个数是需要思考的,为什么这样列方程,这样列方程为什么是正确的都是需要思考的,为了防止文章写的又臭又长,就把这部分省略了,再说也不好说清楚:):
m = 1 时,dp(1, i) = a[i], 1 ≤ i ≤ n。(注:数组a 保存了长度为n 个整个序列)
2 ≤ m ≤ k 时,dp(m, i) = min{dp(m-1, j)} + a[i] * m, m-1 ≤ j < i ≤ n。
严格按照上面给出的递推方程的边界+伪码写程序,得:
#include<stdio.h> // Time Limit Exceeded: 1010MS
#define MAXN 10000
#define INF 1000000000000000 // 大于(1+1000)*1000/2*1000000即可,inf: infinite
long long int a[MAXN+10];
long long int dp[MAXN+10];
long long int temp[MAXN+10];
int main(void)
{
int n, k;
long long int min, ans; // ans: answer
while(scanf("%d%d", &n, &k) != EOF)
{
for(int i = 1; i <= n; i++)
{
scanf("%d", &a[i]);
dp[i] = temp[i] = a[i];
}
for(int m = 2; m <= k; m++)
{
min = INF;
for(int i = m; i <= n; i++)
{
min = INF;
for(int j = m-1; j < i; j++)
{
if(min > temp[j])
min = temp[j];
}
dp[i] = min + a[i]*m;
}
for(int i = m; i <= n; i++)
temp[i] = dp[i];
}
ans = INF;
for(int i = k; i <= n; i++)
{
if(ans > dp[i])
ans = dp[i];
}
printf("%lld\n", ans);
}
return 0;
}
正如上述程序注释中所写的,直接dp 超时了,因此还是得优化一下,我没有掌握过什么斜率/队列/四边形/树形这类听起来很牛的优化方法,因此只能找出最明显的可优化的部分给做掉了:
#include<stdio.h> // 90MS
#define MAXN 10000
#define INF 1000000000000000 // 大于(1+1000)*1000/2*1000000即可,inf: infinite
long long int a[MAXN+10];
long long int dp[MAXN+10];
long long int temp[MAXN+10];
int main(void)
{
int n, k;
long long int min, ans; // ans: answer
while(scanf("%d%d", &n, &k) != EOF)
{
for(int i = 1; i <= n; i++)
{
scanf("%d", &a[i]);
dp[i] = temp[i] = a[i];
}
for(int m = 2; m <= k; m++)
{
min = INF;
for(int i = m; i <= n; i++)
{
if(min > temp[i-1])
min = temp[i-1];
dp[i] = min + a[i]*m;
}
for(int i = m; i <= n; i++)
temp[i] = dp[i];
}
ans = INF;
for(int i = k; i <= n; i++)
{
if(ans > dp[i])
ans = dp[i];
}
printf("%lld\n", ans);
}
return 0;
}
标签:
原文地址:http://www.cnblogs.com/xpjiang/p/4442885.html