HDU5740 Glorious Brilliance

题意:

给出一张不一定合法的染色图,每次可以交换相邻两点的颜色,问最少多少次能使染色图合法

合法的染色图相邻点的颜色不能相同

题解:

首先要确定原图中是否存在染色图,如果不存在直接输出\(-1\)

其次确定给定的两种颜色数量是否能够染出一张合法的染色图,如果数量不对输出\(-1\)

可以发现对于节点\(u,v\),交换两者的颜色最少需要操作两者\(dist_{u,v}\)次,并且可以保证路径中间点的颜色不变

指针先指向最左边点,然后右指针向右找到第一个和左指针所指颜色不同的点,把右侧点不断和左边的点换,然后把左指针指向右指针所指点,继续操作

所以我们需要将两种颜色的点两两匹配,使得最短路的和最小,用KM匹配即可,最后输出路径,就按上述方法输出即可

注意要对每一块连通图单独处理

view code
//#pragma GCC optimize("O3")
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<bits/stdc++.h>
using namespace std;
function<void(void)> ____ = [](){ios_base::sync_with_stdio(false); cin.tie(0); cout.tie(0);};
const int MAXN = 555;
const int INF = 0x3f3f3f3f;
int n,m,col[MAXN],tpcol[MAXN],ID;
int dist[MAXN][MAXN], pre[MAXN][MAXN];
char buf[MAXN];
vector<int> G[MAXN],pt[MAXN]; bool setColor(int x, int id){ // 染色
queue<int> que;
que.push(x);
tpcol[x] = 0;
while(!que.empty()){
int u = que.front();
que.pop();
pt[id].push_back(u);
for(int v : G[u]){
if(tpcol[v]!=-1){
if(tpcol[v]==tpcol[u]) return false;
continue;
}
tpcol[v] = tpcol[u] ^ 1;
que.push(v);
}
}
return true;
}
void getPath(int x){ // 计算最短路和最短路路径
memset(dist[x]+1,0x3f,n<<2);
dist[x][x] = 0;
pre[x][x] = -1;
queue<int> que;
que.push(x);
while(!que.empty()){
int u = que.front();
que.pop();
for(int v : G[u]){
if(dist[x][v]!=INF) continue;
dist[x][v] = dist[x][u] + 1;
pre[x][v] = u;
que.push(v);
}
}
}
int match[MAXN],visx[MAXN],visy[MAXN],lx[MAXN],ly[MAXN],slack[MAXN];
vector<int> lft,rht;
bool dfs(int u){
visx[u] = true;
for(int v : rht){
if(visy[v]) continue;
int gap = dist[u][v] - lx[u] - ly[v];
if(!gap){
visy[v] = true;
if(match[v]==-1 or dfs(match[v])){
match[v] = u;
return true;
}
}
else slack[v] = min(slack[v],gap);
}
return false;
}
int KM(int id, int tag){
lft.clear(); rht.clear();
for(int x : pt[id]){
if((col[x]^tag)==tpcol[x]) continue;
if(col[x]==0) lft.push_back(x);
else rht.push_back(x);
}
if(lft.size()!=rht.size()) return INF;
for(int x : lft){
lx[x] = INF;
for(int y : rht) lx[x] = min(lx[x],dist[x][y]);
}
for(int x : rht){
match[x] = -1;
ly[x] = 0;
}
for(int x : lft){
for(int y : rht) slack[y] = INF;
while(true){
for(int xx : lft) visx[xx] = false;
for(int yy : rht) visy[yy] = false;
if(dfs(x)) break;
int d = INF;
for(int yy : rht) if(!visy[yy]) d = min(d,slack[yy]);
for(int xx : lft) if(visx[xx]) lx[xx] += d;
for(int yy : rht){
if(!visy[yy]) slack[yy] -= d;
else ly[yy] -= d;
}
}
}
int tot = 0;
for(int y : rht) tot += dist[y][match[y]];
return tot;
}
vector<pair<int,int> > path;
void record(int id, int tag){
KM(id,tag);
for(int y : rht){
int x = match[y];
vector<int> route;
int u = y;
while(u!=-1){
route.push_back(u);
u = pre[x][u];
}
int sz = (int)route.size();
int l = 0;
while(l!=sz-1){
int r = l;
while(col[route[r]]==col[route[r+1]]) r++;
r++;
for(int i = r; i > l; i--) path.push_back(make_pair(route[i],route[i-1]));
swap(col[route[l]],col[route[r]]);
l = r;
}
}
}
void solve(){
scanf("%d %d %s",&n,&m,buf+1);
for(int i = 1; i <= n; i++) col[i] = buf[i] - '0';
for(int i = 1; i <= n; i++) G[i].clear();
for(int i = 1; i <= n; i++) tpcol[i] = -1;
for(int i = 1; i <= m; i++){
int u, v; scanf("%d %d",&u,&v);
G[u].push_back(v);
G[v].push_back(u);
}
ID = 0;
bool flag = false;
for(int i = 1; i <= n; i++) getPath(i);
path.clear();
for(int i = 1; i <= n; i++) if(tpcol[i]==-1){
pt[++ID].clear();
if(!setColor(i, ID)){
flag = true;
break;
}
int c0 = 0, c1 = 0;
int tpc0 = 0, tpc1 = 0;
for(int x : pt[ID]){
if(col[x]==0) c0++;
else c1++;
if(tpcol[x]==0) tpc0++;
else tpc1++;
}
if(c0!=tpc0 and c0!=tpc1){
flag = true;
break;
}
int mt0 = INF, mt1 = INF;
if(c0==tpc0) mt0 = KM(ID,0);
if(c0==tpc1) mt1 = KM(ID,1);
if(mt0==INF and mt1==INF){
flag = true;
break;
}
if(mt0<mt1) record(ID,0);
else record(ID,1);
}
if(flag){
puts("-1");
return;
}
printf("%d\n",path.size());
for(auto p : path) printf("%d %d\n",p.first,p.second);
}
int main(){
int tt;
for(scanf("%d",&tt); tt; tt--) solve();
return 0;
}

HDU5740 Glorious Brilliance【最短路 KM匹配】的更多相关文章

  1. Glorious Brilliance (最短路 + 带权二分图匹配)

    这是一道代码大题.一开始读错题意了,然后理解成直接看上去的那种相邻,然后想不通好久!!! 把不同联通的图分离出来,然后先预处理一下形成之后的相邻图的状态,然后根据01确定哪一些是需要更换状态的,然后建 ...

  2. 【转】KM匹配题集

    转自:http://blog.csdn.net/shahdza/article/details/7779324 [HDU]2255 奔小康赚大钱 模板题★1533 Going Home 模板题★242 ...

  3. HDU 5740 - Glorious Brilliance

    题意: 给出已0 1染色的无向图(不一定联通),一次操作为一对相邻点颜色互换. 问使任意相邻点颜色不同,最少需要多少次操作 分析: 交换两点的代价即为两点间最短路. 故用BFS找出所有点到任意点的最短 ...

  4. KM匹配板子

    /* gyt Live up to every day */ #include<cstdio> #include<cmath> #include<iostream> ...

  5. 【转载】图论 500题——主要为hdu/poj/zoj

    转自——http://blog.csdn.net/qwe20060514/article/details/8112550 =============================以下是最小生成树+并 ...

  6. 【HDOJ图论题集】【转】

    =============================以下是最小生成树+并查集====================================== [HDU] How Many Table ...

  7. hdu图论题目分类

    =============================以下是最小生成树+并查集====================================== [HDU] 1213 How Many ...

  8. HDU图论题单

    =============================以下是最小生成树+并查集====================================== [HDU] 1213 How Many ...

  9. POJ 2516 Minimum Cost (KM最优匹配)

    题意:有N家家店,每家店都对K种货物有需求:同时有M家仓库,对K钟货物有供应.对于每种货物,每个仓库送至每家店都有自己的单位费用.求满足所有店所有货物的最小费用 分析:对于每一种货物,如果总需求大于总 ...

随机推荐

  1. #3使用html+css+js制作网页 番外篇 制作接收php

    使用html+css+js制作网页 番外篇 制作接收php 本系列链接 基础 php语法 例子 本系列链接 #1使用html+css+js制作网站教程 准备 #2使用html+css+js制作网站教程 ...

  2. 解放双手,markdown文章神器,Typora+PicGo+七牛云图床实现自动上传图片

    本文主要分享使用Typora作为Markdown编辑器,PicGo为上传图片工具,使用七牛云做存储来解放双手实现图片的自动化上传与管理.提高写作效率,提升逼格.用过 Markdown 的朋友一定会深深 ...

  3. Kali Debian 修改时区

    1.KaliDebian默认读取时区文件是/etc/localtime(重建软连接,达到修改时区目的) 2./etc/localtime是个软连接,连接到时区文件/usr/share/zoneinfo ...

  4. 【EXP】WINDOWS下如何导出

    有些时候需要在windows下通过远程来导出数据 那么windows下怎么导出呢 例子: exp hr/hr@192.168.1.222:1521/zhang file=d:backup.dmp lo ...

  5. CTFshow-萌新赛杂项_劝退警告

    下载附件 https://www.lanzous.com/i9wocah 下载后得到一个劝退警告.zip 解压得到一张gif图片 使用binwalk分析发现包含zip 于是拿到了一个压缩包 打开后发现 ...

  6. rename命令和批量重命名

    本文为转载文章,转发自 https://blog.csdn.net/GGxiaobai/article/details/53507454 早期版本的rename是C语言版本,如今新的Ubuntu中采用 ...

  7. Django前后端分离项目部署

    vue+drf的前后端分离部署笔记 前端部署过程 端口划分: vue+nginx的端口 是81 vue向后台发请求,首先发给的是代理服务器,这里模拟是nginx的 9000 drf后台运行在 9005 ...

  8. 从synchronized和lock区别入手聊聊java锁机制

    写这篇文章之前,我去百度了一下啥叫锁,百度百科上写道:置于可启闭的器物上,以钥匙或暗码开启.确实我们一般理解的锁就是门锁,密码锁,但是在计算机科学中,锁又是啥,说实话,这个问题我也思考了很久,也没法很 ...

  9. (10)-Python3之--引入

    1.什么是模块 .py文件就是模块 模块名有命名要求: 1.不要以数字.下划线开头.特殊符号.也不要以中文开头. 2.通常来说,都是以字母开头. 3.不要以关键字来命名.内置函数.内置模块.不要以第三 ...

  10. Maven 基础详解

    一.编写pom.xml文件     Maven项目的核心是pom.xml.POM(Project Object Model,项目对象模型)定义了项目的基本信息,用于描述项目如何构建,声明项目依赖等等. ...