BZOJ1237：配对（dp）

最新推荐文章于 2019-06-28 19:58:26 发布

Fat_tu

最新推荐文章于 2019-06-28 19:58:26 发布

阅读量247

点赞数

分类专栏： dp

本文链接：https://blog.csdn.net/q582116859/article/details/79560656

版权

dp 专栏收录该内容

24 篇文章 1 订阅

订阅专栏

该博客介绍了BZOJ1237问题的解决方法，内容涉及将n个a和n个b进行两两配对，使得配对的代价（|a-b|）之和最小。博主推断出每个数不会与超过2位距离的数配对，并给出了基于这个结论的动态规划（dp）解决方案。

摘要由CSDN通过智能技术生成

题面

题意：n个a，n个b，两两配对
代价为abs(a-b)，相同的不能配对
问最小代价

有事没事排排序
被我合理演绎（猜）出的结论：每个数不会与距离超过2的数配对

然后就是，字面意思的dp

#include <iostream>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>

using namespace std;
#define mmst(a, b) memset(a, b, sizeof(a))
#define mmcp(a, b) memcpy(a, b, sizeof(b))
#define imin(a, b) ((a) < (b)) ? (a) : (b)

typedef long long LL;

const int N=100100;
const LL oo=1e15;

int n;
LL a[N],b[N],dp[N];
int A[6]={0,0,1,1,2,2};
int B[6]={1,2,0,2,0,1};
int C[6]={2,1,2,0,1,0};

LL val(int x,int y)
{
    return (x==y) ? oo : abs(x-y);
}

int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    scanf("%lld%lld",&a[i],&b[i]);

    sort(a+1,a+n+1);
    sort(b+1,b+n+1);

    dp[1]=val(a[1],b[1]);
    dp[2]=imin(dp[1]+val(a[2],b[2]),val(a[1],b[2])+val(a[2],b[1]));

    for(int i=3;i<=n;i++)
    {
        dp[i]=oo;
        dp[i]=imin(dp[i],dp[i-1]+val(a[i],b[i]));
        dp[i]=imin(dp[i],dp[i-2]+val(a[i],b[i-1])+val(a[i-1],b[i]));
        for(int j=0;j<6;j++)
        dp[i]=imin(dp[i],dp[i-3]+val(a[i-2],b[i-A[j]])+val(a[i-1],b[i-B[j]])+val(a[i],b[i-C[j]]));
    }
    if(dp[n]<oo)
    cout<<dp[n]<<endl;
    else
    cout<<-1<<endl;

    return 0;
}