HDU5740 Glorious Brilliance【最短路 KM匹配】
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匹配】的更多相关文章
- Glorious Brilliance (最短路 + 带权二分图匹配)
这是一道代码大题.一开始读错题意了,然后理解成直接看上去的那种相邻,然后想不通好久!!! 把不同联通的图分离出来,然后先预处理一下形成之后的相邻图的状态,然后根据01确定哪一些是需要更换状态的,然后建 ...
- 【转】KM匹配题集
转自:http://blog.csdn.net/shahdza/article/details/7779324 [HDU]2255 奔小康赚大钱 模板题★1533 Going Home 模板题★242 ...
- HDU 5740 - Glorious Brilliance
题意: 给出已0 1染色的无向图(不一定联通),一次操作为一对相邻点颜色互换. 问使任意相邻点颜色不同,最少需要多少次操作 分析: 交换两点的代价即为两点间最短路. 故用BFS找出所有点到任意点的最短 ...
- KM匹配板子
/* gyt Live up to every day */ #include<cstdio> #include<cmath> #include<iostream> ...
- 【转载】图论 500题——主要为hdu/poj/zoj
转自——http://blog.csdn.net/qwe20060514/article/details/8112550 =============================以下是最小生成树+并 ...
- 【HDOJ图论题集】【转】
=============================以下是最小生成树+并查集====================================== [HDU] How Many Table ...
- hdu图论题目分类
=============================以下是最小生成树+并查集====================================== [HDU] 1213 How Many ...
- HDU图论题单
=============================以下是最小生成树+并查集====================================== [HDU] 1213 How Many ...
- POJ 2516 Minimum Cost (KM最优匹配)
题意:有N家家店,每家店都对K种货物有需求:同时有M家仓库,对K钟货物有供应.对于每种货物,每个仓库送至每家店都有自己的单位费用.求满足所有店所有货物的最小费用 分析:对于每一种货物,如果总需求大于总 ...
随机推荐
- 利用PHP递归 获取所有的上级栏目
/** * 获取所有的上级栏目 * @param $category_id * @param array $array * @return array * @author 宁佳兵 <meilij ...
- tf.argmax(vector,axis)函数的使用
1.返回值 vector为向量,返回行或列的最大值的索引号: vector为矩阵,返回值是向量,返回每行或每列的最大值的索引号. 2.参数 vector为向量或者矩阵 axis = 0 或1 0:返回 ...
- Deep Learn I'm back.
Intorduction: 时隔好几个月,我准备重新进入Deep Learning 的领域.昨天和老师聊了很多,之前觉得我做的工作就是排列组合,在水论文,灌水.但老师却说:这也是为将来的研究打基础. ...
- Java开发手册之异常日志
1.捕获异常的时候,有一种特殊情况,就是方法体内部所抛出的并不是Exception而是Error,这个时候,上层方法捕获Exception就会失败.所以在某些场合需要捕获更高一级别的Throwable ...
- 【Oracle】10g 11g下载路径
现在直接点击不能下载了 要经过oracle许可才可以下载 如果嫌麻烦可以用迅雷直接下载密码是这个 一般不会动了 大家也不用帮我找回密码了 每次都改 也很麻烦的用迅雷下就不用密码了 下载也不会卡到最后 ...
- playwright自动化项目搭建
这是关于playwright系列介绍的最后一篇.搭建基于 playwright 的自动化项目. GitHub地址: https://github.com/defnngj/playwright-pro ...
- 零基础怎么学Python编程,新手常犯哪些错误?
Python是人工智能时代最佳的编程语言,入门简单.功能强大,深获初学者的喜爱. 很多零基础学习Python开发的人都会忽视一些小细节,进而导致整个程序出现错误.下面就给大家介绍一下Python开发者 ...
- .net code+vue 文件上传
后端技术 .net code 官方文档 https://docs.microsoft.com/zh-cn/aspnet/core/mvc/models/file-uploads?view=aspnet ...
- 网络基础知识之Cisco
2021-01-2313:35:21 大家好,我是新手,刚学完Cisco协议,以后,我会每天都会给大家一些学习笔记 和一些学习心得,这方面不会的可以找我.qq2934896930. 网络的定义: 计算 ...
- jquery 数据查询
jquery 数据查询 <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> & ...