BZOJ 3140: [Hnoi2013]消毒网络流最小点覆盖+状压压维

最新推荐文章于 2024-08-16 23:08:28 发布

dilu0653

最新推荐文章于 2024-08-16 23:08:28 发布

阅读量120

点赞数

文章标签：数据结构与算法人工智能

原文链接：http://www.cnblogs.com/G-hsm/p/11329810.html

版权

title

BZOJ 3140
LUOGU 3231
Description

最近在生物实验室工作的小T遇到了大麻烦。
由于实验室最近升级的缘故，他的分格实验皿是一个长方体,其尺寸为abc，a、b、c 均为正整数。为了实验的方便，它被划分为abc个单位立方体区域，每个单位立方体尺寸为111。用(i,j,k)标识一个单位立方体，1 ≤i≤a，1≤j≤b，1≤k≤c。这个实验皿已经很久没有人用了，现在，小T被导师要求将其中一些单位立方体区域进行消毒操作（每个区域可以被重复消毒）。而由于严格的实验要求，他被要求使用一种特定的F试剂来进行消毒。这种F试剂特别奇怪，每次对尺寸为xyz的长方体区域（它由xyz个单位立方体组成）进行消毒时，只需要使用min{x,y,z}单位的F试剂。F试剂的价格不菲，这可难倒了小 T。现在请你告诉他，最少要用多少单位的F试剂。(注：min{x,y,z}表示x、y、z中的最小者。)

Input

第一行是一个正整数D，表示数据组数。接下来是D组数据，每组数据开头是三个数a,b,c表示实验皿的尺寸。接下来会出现a个b 行c列的用空格隔开的01矩阵，0表示对应的单位立方体不要求消毒，1表示对应的单位立方体需要消毒；例如，如果第1个01矩阵的第2行第3列为1，则表示单位立方体(1,2,3)需要被消毒。输入保证满足abc≤5000,T≤3。

Output

仅包含D行，每行一个整数，表示对应实验皿最少要用多少单位的F试剂。

Sample Input

1
4 4 4
1 0 1 1
0 0 1 1
0 0 0 0
0 0 0 0
0 0 1 1
1 0 1 1
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0

Sample Output

3

HINT

对于区域(1,1,3)-(2,2,4)和(1,1,1)-(4,4,1)消毒，分别花费2个单位和1个单位的F试剂。2017.5.26新加两组数据By Leoly，未重测.

analysis

星期一晚上晚自习没去，和 \(Chdy\) 一起连续讨论了几道网络流题目。

对于这道题，我比较快地看懂了题意，由于觉得和切糕模型很像，又加上刚写过交换棋子，便想着把两个模型放在一起，稍加处理应该是可以的，然而 \(Chdy\) 指出：虽然切糕模型在题目中是三维的，但是在实际处理中其实是只考虑切面，也就是二维的，所以这道题可能是要状压压维的。蒟蒻表示同意。

其后，\(Chdy\) 想到这会不是一个最小割转最大流模型，因为我们其实可以一层层消，这样每消一个小长方体（长宽高任意一个为 1 ），那么花费就是 1 了，那么我们其实可以把这些小长方体之间连边，跑出最小割，就是能够把这些小长方体切断联系的最小费用，应该就是答案了。

然而当我们打开题解看思路是否正确时，额，清一色的二分图匹配：匈牙利算法。不禁懵逼。还好翻到最后， \(⚡cdecl⚡\) 写了种网络流写法，嗯，然而...他怎么分析的是最小路径覆盖？

额，这样的话，我们分析出来的模型岂不是错了。怎么回事呢？为了更好理解，我们跳到了那些使用匈牙利算法的大佬那里看了分析，

这道题的正确模型应该是抽象成：

二维平面有一些方块被染了黑色，你每次可以选择 \((x,y)\) 的区域染成白色，代价是 \(\min(x,y)\)，问你付出的最小代价

所以你肯定不会直接把整个面直接涂白，那样太浪费了，肯定是只把这些点所在的那一条边（这样说好理解些吧）涂白，那就是最小路径覆盖了：用最少的边覆盖所有点。

这样一来，我们的最小割模型的问题就是我们只看到了三维情况，并没有来处理其中所包含的二维情况，导致如果我们真的那样处理，还是需要看二维状态下怎么处理，额，可见，简化问题很重要，三维的一般能转成二维。

关于最小路径覆盖的建图，基本都会（行列模型，当然要先维护一下长宽高，将最小的维度当成多出的一维），但是多了一个维度，怎么搞？

那就把这一维度状压掉，暴力枚举所有状态中的某一层是直接削掉还是一会儿再处理（只有这两种情况，别的都不好，想一想为什么）。

大概状态数为 \(z \leq 17\)，因为 \(x*y*z\leq 5000\)，这样 \(z\leq \sqrt[3]{5000}\)，而 \(\sqrt[3]{5000} \approx 17\)。

对于没有被直接削掉的层，我们把它们剥离出来，然后拍扁成二维平面上的问题求解。

吼，跑一遍最大流就好了。

在此期间， \(Chdy\)大佬发现这和之前的一道错题非常像：2711 小行星。这道题便是简单的行列模型连边，求出最小割即可（怪不得你能想到最小割，原来如此，(￣^￣)。）不过已经被 \(Chdy\) 证明是错题了，我们这种连边是无法识别这条流是从哪里流出的，因此会出现多流的情况，具体的解释看帖子吧：此题有误或者该题数据有误。

只能说还好这不是某谷原创题，不然丢脸丢大了。

code

~~先上个简单的，小行星（好像不太一样）~~
很多人说需要拆点，不过应该是不需要的，标记一下就行。

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+10,maxm=1e6+10,inf=0x3f3f3f3f;

char buf[1<<15],*fs,*ft;
inline char getc() { return (ft==fs&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),ft==fs))?0:*fs++; }
template<typename T>inline void read(T &x)
{
    x=0;
    T f=1, ch=getchar();
    while (!isdigit(ch) && ch^'-') ch=getchar();
    if (ch=='-') f=-1, ch=getchar();
    while (isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48), ch=getchar();
    x*=f;
}

template<typename T>inline void write(T x)
{
    if (!x) { putchar('0'); return ; }
    if (x<0) putchar('-'), x=-x;
    T num=0, ch[20];
    while (x) ch[++num]=x%10+48, x/=10;
    while (num) putchar(ch[num--]);
}

int ver[maxm],edge[maxm],Next[maxm],head[maxn],len=1;
inline void add(int x,int y,int z)
{
    ver[++len]=y,edge[len]=z,Next[len]=head[x],head[x]=len;
    ver[++len]=x,edge[len]=0,Next[len]=head[y],head[y]=len;
}

int s,t;
int dist[maxn];
inline bool bfs()
{
    queue<int>q;
    memset(dist,0,sizeof(dist));
    q.push(s),dist[s]=1;
    while (!q.empty())
    {
        int x=q.front();
        q.pop();
        for (int i=head[x]; i; i=Next[i])
        {
            int y=ver[i];
            if (edge[i] && !dist[y])
            {
                dist[y]=dist[x]+1;
                if (y==t) return 1;
                q.push(y);
            }
        }
    }
    return 0;
}

inline int get(int x,int low)
{
    if (x==t) return low;
    int tmp=low;
    for (int i=head[x]; i; i=Next[i])
    {
        int y=ver[i];
        if (edge[i] && dist[y]==dist[x]+1)
        {
            int a=get(y,min(tmp,edge[i]));
            if (!a) dist[y]=0;
            edge[i]-=a;
            edge[i^1]+=a;
            if (!(tmp-=a)) break;
        }
    }
    return low-tmp;
}

bool v1[510],v2[510];
int main()
{
    int n;read(n);s=0,t=1501;
    for (int i=1,x,y,z; i<=n; ++i)
    {
        read(x),read(y),read(z);
        if (!v1[x]) add(s,x,1),v1[x]=1;
        if (!v2[z]) add(z+1000,t,1),v2[z]=1;
        add(x,y+500,1),add(y+500,z+1000,1);
    }
    int ans=0;
    while (bfs()) ans+=get(s,inf);
    write(ans);
    return 0;
}

~~难的来了~~
\(T\) 了快一版（加上 \(Chdy\) 交的就一版了），所以切记：开 \(O^2\) 。

#include<bits/stdc++.h>
using namespace std;
const int maxn=5e3+10,maxm=5e4+10,inf=0x3f3f3f3f;

char buf[1<<15],*fs,*ft;
inline char getc() { return (ft==fs&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),ft==fs))?0:*fs++; }
template<typename T>inline void read(T &x)
{
    x=0;
    T f=1, ch=getchar();
    while (!isdigit(ch) && ch^'-') ch=getchar();
    if (ch=='-') f=-1, ch=getchar();
    while (isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48), ch=getchar();
    x*=f;
}

template<typename T>inline void write(T x)
{
    if (!x) { putchar('0'); return ; }
    if (x<0) putchar('-'), x=-x;
    T num=0, ch[20];
    while (x) ch[++num]=x%10+48, x/=10;
    while (num) putchar(ch[num--]);
}

int ver[maxm<<1],edge[maxm<<1],Next[maxm<<1],head[maxn],len=1;
inline void add(int x,int y,int z)
{
    ver[++len]=y,edge[len]=z,Next[len]=head[x],head[x]=len;
    ver[++len]=x,edge[len]=0,Next[len]=head[y],head[y]=len;
}

int ss,tt;
int dist[maxn];
inline bool bfs()
{
    queue<int>q;
    memset(dist,0,sizeof(dist));
    q.push(ss),dist[ss]=1;
    while (!q.empty())
    {
        int x=q.front();
        q.pop();
        for (int i=head[x]; i; i=Next[i])
        {
            int y=ver[i];
            if (edge[i] && !dist[y])
            {
                dist[y]=dist[x]+1;
                if (y==tt) return 1;
                q.push(y);
            }
        }
    }
    return 0;
}

inline int get(int x,int low)
{
    if (x==tt) return low;
    int tmp=low;
    for (int i=head[x]; i; i=Next[i])
    {
        int y=ver[i];
        if (edge[i] && dist[y]==dist[x]+1)
        {
            int a=get(y,min(tmp,edge[i]));
            if (!a) dist[y]=0;
            edge[i]-=a;
            edge[i^1]+=a;
            if (!(tmp-=a)) break;
        }
    }
    return low-tmp;
}

inline void Clear()
{
    memset(head,0,sizeof(head));
    len=1;
}

inline int Count(int x)//统计二进制数中1的个数
{
    int num=0;
    while (x) ++num, x&=(x-1);
    return num;
}

struct Orz{int x,y,z;}p[maxn];
int cnt=0;
int main()
{
    int T;read(T);
    while (T--)
    {
        int a,b,c,ans=inf;
        cnt=0;
        read(a);read(b);read(c);
        for (int i=1; i<=a; ++i)
            for (int j=1; j<=b; ++j)
                for (int k=1,x; k<=c; ++k)
                {
                    read(x);
                    if (x) p[++cnt]=(Orz){i,j,k};
                }
        if (a>b)
        {
            swap(a,b);
            for (int i=1; i<=cnt; ++i) swap(p[i].x,p[i].y);
        }
        if (a>c)
        {
            swap(a,c);
            for (int i=1; i<=cnt; ++i) swap(p[i].x,p[i].z);
        }
        ss=0,tt=maxn-1;
        for (int s=0; s<=(1<<a)-1; ++s)
        {
            Clear();
            int lans=a-Count(s),mx=0;
            for (int i=1; i<=cnt; ++i)
            {
                mx=max(mx,max(p[i].y,p[i].z));
                if (1&(s>>p[i].x-1)) add(p[i].y<<1,p[i].z<<1|1,1);
            }
            for (int i=1; i<=mx; ++i) add(ss,i<<1,1),add(i<<1|1,tt,1);
            while (bfs() && lans<ans) lans+=get(ss,inf);
            ans=min(ans,lans);
        }
        write(ans),puts("");
    }
    return 0;
}

转载于:https://www.cnblogs.com/G-hsm/p/11329810.html

dilu0653

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
BZOJ 3140: [Hnoi2013]消毒网络流最小点覆盖+状压压维

titleBZOJ 3140LUOGU 3231Description最近在生物实验室工作的小T遇到了大麻烦。由于实验室最近升级的缘故，他的分格实验皿是一个长方体,其尺寸为abc，a、b、c 均为正整数。为了实验的方便，它被划分为abc个单位立方体区域，每个单位立方体尺寸为111。用(i,j,k)标识一个单位立方体，1 ≤i≤a，1≤j≤b，1≤k≤c。这个实验皿已经很久没有...
复制链接

扫一扫