hdu1853解题报告
题意和解决回路匹配的思路如同hdu3488
(这里我第一次想到最短路,但是对于有回路这个不知道怎么处理,后来看了别人的解题报告才知道KM匹配,但是看到KM之后就自己想...想了很久....还是不知道回路这个地方怎么匹配......其实应该这样来想....总共有N个城市....如果是要形成回路..那么就是环,那么每一个城市都要和指向的城市匹配一次,也要被一个城市指向自己匹配一次...那么匹配的时候我们把所有的N个城市分成两拨,对于每一个城市都匹配一次就得到了一个完全匹配(既每个点都匹配一次,也就是每一个城市匹配一次、被匹配一次).....那么还有一个地方要注意的是:这里有重边...我们应该选择小的边更新)
这里唯一有一点点不同的就是可以存在无回路的情况,那就是有一个点或者多个点没有匹配到.....
那么在代码中体现出来:
// 31MS 272K
#include<stdio.h>
#include<string.h> #define MAX 101
#define INF 1<<30-1 int N,M;
int map[MAX][MAX];
int lx[MAX],ly[MAX],link[MAX],slar[MAX];
bool visx[MAX],visy[MAX]; bool dfs(int x)
{
visx[x] = true;
for(int y = 1; y <= N; y ++)
{
if(visy[y]) continue;// || x == y
int t = lx[x] + ly[y] - map[x][y];
if(t == 0)
{
visy[y] = true;
if(link[y] == -1 || dfs(link[y]))
{
link[y] = x;return true;
}
}
else if(slar[y] > t) slar[y] = t;
}
return false;
} int KM()
{
int i,j;
memset(ly,0,sizeof(ly));
memset(link,-1,sizeof(link));
for(i = 1 ;i <= N;i ++)
{
lx[i] = -INF;
for(j = 1; j <= N; j ++)
if(lx[i] < map[i][j]) lx[i] = map[i][j];
}
for(int x = 1; x <=N; x ++)
{
for(i = 1; i <= N ; i ++) slar[i] = INF;
while(1)
{
memset(visx,false,sizeof(visx));
memset(visy,false,sizeof(visy));
if(dfs(x)) break; int d = INF;
for(i = 1; i <= N;i ++)
if(!visy[i] && d > slar[i]) d = slar[i];
for(i = 1; i <= N; i ++)
if(visx[i]) lx[i] -= d;
for(i = 1; i <= N; i ++)
if(visy[i]) ly[i] += d;
else slar[i] -= d;
}
}
int ans = 0;bool flag = false;
for(i = 1; i <= N; i ++)
if(link[i] == -1 || map[link[i]][i] == -INF)//这里判断匹配不到的情况
{
flag=true;break;
}
if(flag) ans = 1;
return -ans;
} int main()
{
int i,j;
int a,b,c;
while(~scanf("%d%d",&N,&M))
{
for(i = 1; i <= N; i ++)
for(j = 1; j <= N; j ++)
map[i][j]=-INF;
while(M -- )
{
scanf("%d%d%d",&a,&b,&c);
map[a][b] = map[a][b] > -c ? map[a][b] : -c;//重边选择
}
printf("%d\n",KM());
}
return 0;
}
个人愚昧观点..欢迎指正与讨论
hdu1853解题报告的更多相关文章
- CH Round #56 - 国庆节欢乐赛解题报告
最近CH上的比赛很多,在此会全部写出解题报告,与大家交流一下解题方法与技巧. T1 魔幻森林 描述 Cortana来到了一片魔幻森林,这片森林可以被视作一个N*M的矩阵,矩阵中的每个位置上都长着一棵树 ...
- 二模13day1解题报告
二模13day1解题报告 T1.发射站(station) N个发射站,每个发射站有高度hi,发射信号强度vi,每个发射站的信号只会被左和右第一个比他高的收到.现在求收到信号最强的发射站. 我用了时间复 ...
- BZOJ 1051 最受欢迎的牛 解题报告
题目直接摆在这里! 1051: [HAOI2006]受欢迎的牛 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4438 Solved: 2353[S ...
- 习题:codevs 2822 爱在心中 解题报告
这次的解题报告是有关tarjan算法的一道思维量比较大的题目(真的是原创文章,希望管理员不要再把文章移出首页). 这道题蒟蒻以前做过,但是今天由于要复习tarjan算法,于是就看到codevs分类强联 ...
- 习题:codevs 1035 火车停留解题报告
本蒟蒻又来写解题报告了.这次的题目是codevs 1035 火车停留. 题目大意就是给m个火车的到达时间.停留时间和车载货物的价值,车站有n个车道,而火车停留一次车站就会从车载货物价值中获得1%的利润 ...
- 习题: codevs 2492 上帝造题的七分钟2 解题报告
这道题是受到大犇MagHSK的启发我才得以想出来的,蒟蒻觉得自己的代码跟MagHSK大犇的代码完全比不上,所以这里蒟蒻就套用了MagHSK大犇的代码(大家可以关注下我的博客,友情链接就是大犇MagHS ...
- 习题:codevs 1519 过路费 解题报告
今天拿了这道题目练练手,感觉自己代码能力又增强了不少: 我的思路跟别人可能不一样. 首先我们很容易就能看出,我们需要的边就是最小生成树算法kruskal算法求出来的边,其余的边都可以删掉,于是就有了这 ...
- NOIP2016提高组解题报告
NOIP2016提高组解题报告 更正:NOIP day1 T2天天爱跑步 解题思路见代码. NOIP2016代码整合
- LeetCode 解题报告索引
最近在准备找工作的算法题,刷刷LeetCode,以下是我的解题报告索引,每一题几乎都有详细的说明,供各位码农参考.根据我自己做的进度持续更新中...... ...
随机推荐
- php如何在原来的时间上加一天?一小时
php如何在原来的时间上加一天?一小时? <?phpecho "今天:",date('Y-m-d H:i:s'),"<br>";echo &q ...
- Sicily-1156
一. 模仿树的的先序遍历.范围是1000个节点.用数组存储节点的信息. 二. 要注意的是,头结点是不确定的,所以在前序遍历之前要找出头结点,除了头结点的下标值出现一次之外,其他结点 ...
- 大一C语言结课设计之《简单计算器》
/*===============================================*\ ** 设计目的:简单计算器,计算形如10*(20.2-30.6)+5.0/2的表达式值 ** 简 ...
- ASP.NET Calendar 控件
ASP.NET Calendar 控件 http://www.w3school.com.cn/aspnet/control_calendar.asp
- awk参数详解
wk是行处理器: 相比较屏幕处理的优点,在处理庞大文件时不会出现内存溢出或是处理缓慢的问题,通常用来格式化文本信息 awk处理过程: 依次对每一行进行处理,然后输出 awk命令形式: awk [-F| ...
- Chrome Firefox 自定义背景色
Chrome有个自定义背景色的文件 Custom.css 默认里面什么字都没写 html, body {background-color: #e0dcc0!important;} 这个颜色 ...
- PDF417码
PDF417码是由留美华人王寅敬(音)博士发明的.PDF是取英文Portable Data File三个单词的首字母的缩写,意为“便携数据文件”.因为组成条码的每一符号字符都是由4个条和4个空构成,如 ...
- 学校作业-Dp练习
题目 ★Stringsobits01 串 考虑排好序的 N(N<=31)位二进制数. 你会发现,这很有趣.因为他们是排列好的,而且包含所有可能的长度为 N 且含有 1 的个数小于等于 L(L ...
- [每日一题] 11gOCP 1z0-052 :2013-09-14 repeated parsing activity.................................A70
转载请注明出处:http://blog.csdn.net/guoyjoe/article/details/11699299 正确答案:D SQL语句的执行过程: 1.客户端输入sql语句update ...
- HDU 4720 Naive and Silly Muggles (外切圆心)
Naive and Silly Muggles Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Oth ...