HDU 5352 MZL's City

题意:
给你n个点,m次操作,和整数k
1操作,每次可以修复k个属于x所在连通块内的点
2操作,修复x与y之间的路,双向边
3操作,给你一个p,表示有p条路被毁坏了
思路:
其实说实话,如果没有想到关键部分,个人觉得不是特别好做。
第一想法比较接近于网络流或者贪心,但是因为要删边,所以很多人觉得网络流不能解决。
之后看了题解之后知道是二分图的最大匹配,才想通。

对于每一个1操作,我们将1操作当做一个点,并拆成k个点,然后每个点和x的连通块内的所以点连边。
对于每一个2操作,暴力加边。
对于每一个3操作,暴力删边。(实际上因为N比较小所以我用了邻接矩阵)

然后从最后一个操作跑匈牙利算法,这样可以保证后面的匹配数更大,从而得到字典序最小的。

//      whn6325689
//      Mr.Phoebe
//      http://blog.csdn.net/u013007900
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
#include <functional>
#include <numeric>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;

#define eps 1e-9
#define PI acos(-1.0)
#define INF 0x3f3f3f3f
#define LLINF 1LL<<62
#define speed std::ios::sync_with_stdio(false);

typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef pair<ll, ll> pll;
typedef complex<ld> point;
typedef pair<int, int> pii;
typedef pair<pii, int> piii;
typedef vector<int> vi;

#define CLR(x,y) memset(x,y,sizeof(x))
#define CPY(x,y) memcpy(x,y,sizeof(x))
#define clr(a,x,size) memset(a,x,sizeof(a[0])*(size))
#define cpy(a,x,size) memcpy(a,x,sizeof(a[0])*(size))

#define mp(x,y) make_pair(x,y)
#define pb(x) push_back(x)
#define lowbit(x) (x&(-x))

#define MID(x,y) (x+((y-x)>>1))
#define ls (idx<<1)
#define rs (idx<<1|1)
#define lson ls,l,mid
#define rson rs,mid+1,r
#define root 1,1,n

template<class T>
inline bool read(T &n)
{
    T x = 0, tmp = 1;
    char c = getchar();
    while((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
    if(c == EOF) return false;
    if(c == '-') c = getchar(), tmp = -1;
    while(c >= '0' && c <= '9') x *= 10, x += (c - '0'),c = getchar();
    n = x*tmp;
    return true;
}
template <class T>
inline void write(T n)
{
    if(n < 0)
    {
        putchar('-');
        n = -n;
    }
    int len = 0,data[20];
    while(n)
    {
        data[len++] = n%10;
        n /= 10;
    }
    if(!len) data[len++] = 0;
    while(len--) putchar(data[len]+48);
}
//-----------------------------------

const int MAXN=505*205+205;

int n,m,k;
vi g[MAXN];
int from[MAXN],tot;
int used[MAXN];
int ans[505];
int op[505];
int qq[505],cntt;
int tmstamp;

bool match(int x)
{
    int sz=g[x].size();
    for(int i=0;i<sz;i++)
        if(used[g[x][i]]!=tmstamp)
        {
            used[g[x][i]]=tmstamp;
            if(from[g[x][i]]==-1 || match(from[g[x][i]]))
            {
                from[g[x][i]]=x;
                return true;
            }            
        }
    return false;
}

int hungry()
{
    tot=0;
    CLR(from,-1);
    for(int i=cntt-1;i>=0;i--)
        for(int j=qq[i]*k;j<(qq[i]+1)*k;j++)
        {
            tmstamp++;
            if(match(j))
            {
                ans[qq[i]]++;tot++;    
            }
        }
    return tot;
}

int ma[222][222];
int vis[222],q[222],cnt;
void dfs(int u,int op)
{
    vis[u]=1;q[cnt++]=u;
    for(int i=1;i<=n;i++)
        if(ma[u][i] && !vis[i])
                dfs(i,op);
}

void init()
{
    CLR(ma,0);CLR(ans,0);cntt=0;
    for(int i=0;i<MAXN;i++)
        g[i].clear();
}

int main()
{
    //freopen("data.txt","r",stdin);
    int T;
    read(T);
    tmstamp=0;
    while(T--)
    {
        init();
        read(n),read(m),read(k);
        for(int i=1,x,y,tt;i<=m;i++)
        {
            read(op[i]);
            if(op[i]==1)
            {
                qq[cntt++]=i;
                read(x);
                CLR(vis,0);cnt=0;
                dfs(x,i);
                for(int t=0;t<cnt;t++)
                    for(int j=i*k;j<(i+1)*k;j++)
                    {
                        //cout<<i<<":("<<j<<") "<<x<<endl;
                        g[j].pb(q[t]);    
                    }
            }
            else if(op[i]==2)
            {
                read(x),read(y);
                ma[x][y]=ma[y][x]=1;
            }
            else
            {
                read(tt);
                for(int i=1;i<=tt;i++)
                {
                    read(x),read(y);
                    ma[x][y]=ma[y][x]=0;
                }
            }
        }
        int sum=hungry();
        write(sum),putchar('\n');
        bool flag=false;
        for(int i=0;i<cntt;i++)
        {
            printf("%s%d",((flag)?" ":""),ans[qq[i]]);
            flag=true;
        }
        puts("");
    }
    return 0;
} 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值