每周学习总结:第六周;
本周并没有学太多新的东西,更多的是对线性dp的精进,就是刷题了。所以这次总结主要是我对题目的一些看法和题思路。
这一周一共打了4场codeforces ,已经对比赛有了大致的的理解。
一: vjudge 的 S题
解题思路:
题目意思为求n个物品,拿k对使得消耗的体力最少,或者说是这k对物品,每一对中两件物品的质量差平方最小,所以要使得质量差的平方小,只能排序后取质量相邻两个物品作为一对;
现在设f[i][j]为前i件物品组成k对所消耗的体力最小;
这时分两种情况含有第i件物品和不含有第i件物品(即第i件物品是不是含在第j对里)
1.含有i件物品 则有 dp[i][j]=dp[i-2][j-1]+(a[i]-a[i-1])(a[i]-a[i-1]);
2.不含第i件物品则有 dp[i][j]=dp[i-1][j];
所以动态转移方程为:dp[i][j]=min(dp[i-2][j-1]+(a[i]-a[i-1])(a[i]-a[i-1]), dp[i-1][j]);
代码:
#include<iostream>
#include<algorithm>
#include<cmath>
#include <cstdio>
#include <cstring>
using namespace std;
#define INF 2147483646
int dp[2010][1010];
int main()
{ int a[2010];
int n,k,i,j;
memset(a,0,sizeof(a));
while(scanf("%d%d",&n,&k)!=EOF)
{
a[0]=0;
for(i=1; i<=n; i++)
scanf("%d",&a[i]);
sort(a+1,a+n+1);
for(i=0; i<=n; i++)
for(j=1; j<=k; j++)
dp[i][j]=INF;
dp[0][0]=0; dp[1][0]=0;
for(i=2; i<=n; i++){
dp[i][0]=0;
for(j=1; j*2<=i; j++)
dp[i][j]=min(dp[i-2][j-1]+(a[i]-a[i-1])*(a[i]-a[i-1]),dp[i-1][j]);
}
printf("%d\n",dp[n][k]);
}
return 0;
}
此题较为常规,得出状态转移方程就完成了大半。
我的收获:1、memset 函数 最好用来赋值 0,-1,不要赋值很大的数,会出错。
二 vjudge 的 L题。 (dp +离散化)。
思路:
预处理:因为a[i]的范围是1到1e9,如果直接dp序列严格递增,那么dp就得从1循环到1e9,这样无论是在空间还是时间上都是不允许的。
1、、、、
(1) 我们可以通过a[i]-=i这步操作,将序列转化为非严格递增序列,那么答案就可以在给定序列中找到,因为把一串数都变成相等的数时,最优解一定是变成他们的中位数,而中位数一定是存在于原有数中的,
(2)再将这n个数进行离散化处理,(离散化处理:可以用b[]存储给定序列a[]-=i的数,并进行排序),这样dp就最多只需要从1循环到n,那么时间复杂的就优化到了O(n^2)。
2、、、
状态表示:dp[i][j] //表示前i个数,这i个数进行操作后最大的数不超过b[j] 的操作次数的最小值(b[]即为离散化后的序列)
状态计算:有两种方法得到可以得到dp[i][j]:
(1)、、、当将第i个数操作后为b[j]时,dp[i][j]为dp[i-1][j]+abs(a[i]-b[j])。
(2)、、、当将第i个数操作后值小于b[j]时,f[i][j]为min(dp[i][k])(1<= k < j)。这样的话时间复杂度就变成了O(n^3)。不过在这里我们可以用一个优化操作:即让 dp[i][j]为dp[i][j-1]
因为虽然每次只是求了dp[i][j-1],但因为前面的dp[i][1-j-1]都求得了最小值,因此只需要算出前面一个的即可。
因此,状态转移方程为:dp[i][j]=min(dp[i][j-1],dp[i-1][j]+abs(a[i]-b[j])).
代码:
#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#define LL long long
using namespace std;
int a[3010],b[3010];
LL dp[3010][3010];
int main()
{
int n;
cin>>n;
for(int i=1;i<=n;i++) //预处理+离散化
{
cin>>a[i];
a[i]-=i;
b[i]=a[i];
}
sort(b+1,b+1+n);
for(int i=1;i<=n;i++)
{
dp[i][0]=1e18; //初始化 : 定边界。
for(int j=1;j<=n;j++)
{
dp[i][j]=min(dp[i][j-1],dp[i-1][j]+abs(a[i]-b[j]));
}
}
LL ans=1e18;
for(int i=1;i<=n;i++) //答案为f[n][i]中的一个,要遍历寻找最小值
{
ans=min(ans,dp[n][i]);
}
cout<<ans<<endl;
return 0;
}
看法:这道题用到了离散化,正好老师前不久讲过。
离散化,是把无限空间中有限的个体映射到有限的空间中去,以此提高算法的时空效率。通俗的说,离散化是在不改变数据相对大小的条件下,对数据进行相应的缩小。
或者第二种方法:用优先队列 来解题;(对原题进行优化)(我在网上找到一种很奇妙的方法,看了半天还是没看懂原理,先把代码贴上吧,回头在补)
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#define LL long long
using namespace std;
const int N=3e3+5,mod=1e9+7;
priority_queue<int>q;
int main()
{
int n;
LL ans=0;
cin>>n;
for(int i=1;i<=n;i++)
{
int x;
cin>>x;
x-=i;
q.push(x);
int max=q.top();
if(max>x)
{
ans+=max-x;
q.pop();
q.push(x);
}
}
cout<<ans<<endl;
return 0;
}
三:codeforces总结:
本周打了3场div2,一场 div 3; rating 也到了 900分。
div2 的A,B题 一般看清楚题意,稍微用合理的方法就可以过。
而之后的题 都比较难想一点,还会涉及到一些算法(目前我不用想/–/)
通过比赛发现自己的不足:
1、思想单一,转换思路不灵活,一般我只要一碰壁 这道题就没了。
2、知识量储备不足,用到一个类型变量,不能很清晰的知道它的相关函数。导致时间过长。、
自己要走的路还有很长,要永远往前看。
下周目标:
1、 每天1~2题 ,周末 2~3题 !!、
2、 还差300分绿名, 争取这一周或下一周 上绿名!!
🆗