vjudge传送
这题与其说是状压dp,到更像是插头dp。
看到\(B\)那么大,肯定要先离散化,然后再预处理出来每一列牛的位置,即有一头还是两头,在上面还是下面。
转移其实比较显然,就是有点恶心。
\(dp[i][j]\)中的\(j\)分四种情况:
1.棚子只建在第一行。
2.棚子只建在第二行。
3.两个棚子分别建在两行。
4.一个棚子覆盖了两行。
那么我们需要根据这一列牛的位置考虑每一个状态怎么转移。比如如果只有上面有一头牛,那么\(dp[i][1]\)可以直接新建一个棚子,也可以从\(dp[i-1][1],dp[i-1][3]\)延伸过来。需要注意的是,延伸要算上和上一列的牛之间的距离,而新建增加的面积是1.
其他情况与此类似,详见代码吧。
debug挺久的。
#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<vector>
#include<queue>
#include<assert.h>
#include<ctime>
using namespace std;
#define enter puts("")
#define space putchar(' ')
#define Mem(a, x) memset(a, x, sizeof(a))
#define In inline
#define forE(i, x, y) for(int i = head[x], y; ~i && (y = e[i].to); i = e[i].nxt)
typedef long long ll;
typedef double db;
const ll INF = 1e14;
const db eps = 1e-8;
const int maxn = 1e3 + 5;
In ll read()
{
ll ans = 0;
char ch = getchar(), las = ' ';
while(!isdigit(ch)) las = ch, ch = getchar();
while(isdigit(ch)) ans = (ans << 1) + (ans << 3) + ch - '0', ch = getchar();
if(las == '-') ans = -ans;
return ans;
}
In void write(ll x)
{
if(x < 0) x = -x, putchar('-');
if(x >= 10) write(x / 10);
putchar(x % 10 + '0');
}
In void MYFILE()
{
#ifndef mrclr
freopen("random.in", "r", stdin);
freopen("ha.out", "w", stdout);
#endif
}
int n, m, B;
int l[maxn], pos[maxn];
int li[maxn], _n, buc[maxn], num[maxn]; //buc: 0:up 1:down 2:both
ll dp[maxn][maxn][5];
int main()
{
// MYFILE();
while(scanf("%d%d%d", &n, &m, &B) != EOF)
{
Mem(buc, 0);
for(int i = 1; i <= n; ++i)
{
l[i] = read(), pos[i] = read();
li[i] = pos[i];
}
sort(li + 1, li + n + 1);
_n = unique(li + 1, li + n + 1) - li - 1;
for(int i = 1; i <= n; ++i)
{
int p = lower_bound(li + 1, li + _n + 1, pos[i]) - li;
if(l[i] == 2) buc[p] = 1;
if(++num[p] > 1) buc[p] = 2;
}
Mem(dp, 0x3f);
dp[0][0][3] = 0;
li[0] = li[1] - 1;
for(int i = 1; i <= _n; ++i)
{
int dis = li[i] - li[i - 1];
for(int j = 1; j <= m; ++j)
{
ll Min1 = INF, Min2 = INF;
for(int k = 1; k <= 4; ++k)
{
Min1 = min(Min1, dp[i - 1][j - 1][k]);
Min2 = min(Min2, j - 2 >= 0 ? dp[i - 1][j - 2][k] : INF);
}
if(buc[i] == 0 || buc[i] == 1)
{
int k = buc[i] + 1;
dp[i][j][k] = min(dp[i][j][k], min(dp[i - 1][j][k], dp[i - 1][j][3]) + dis);
dp[i][j][k] = min(dp[i][j][k], Min1 + 1);
}
dp[i][j][3] = min(dp[i][j][3], min(dp[i - 1][j - 1][1], dp[i - 1][j - 1][2]) + dis + 1);
dp[i][j][3] = min(dp[i][j][3], min(dp[i - 1][j][3] + (dis << 1), Min2 + 2));
dp[i][j][4] = min(dp[i][j][4], min(Min1 + 2, dp[i - 1][j][4] + (dis << 1)));
}
}
ll ans = INF;
for(int j = 1; j <= 4; ++j)
ans = min(ans, dp[_n][m][j]);
write(ans), enter;
}
return 0;
}