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,以下是我的解题报告索引,每一题几乎都有详细的说明,供各位码农参考.根据我自己做的进度持续更新中...... ...
随机推荐
- Microsoft SQL Server 混合云博客系列
Microsoft 云操作系统愿景的核心支柱之一就是借助我们的混合云基础结构改造数据中心.在 Windows Azure 基础结构服务正式发布后的几个月里,我们一直在发布博客,介绍 Windows A ...
- leetCode 24. Swap Nodes in Pairs (双数交换节点) 解题思路和方法
Swap Nodes in Pairs Given a linked list, swap every two adjacent nodes and return its head. For exa ...
- BI系统的应用组织思路与数据分析模式
BI商业智能软件一般都会提供若干数据整合.数据查询.分析与评价.数据可视化及数据分享的手段,但是在BI项目的构建与实施过程中,如果不按照一定的应用组织思路.数据分析模式及分析流程使用这些工具或手段,呈 ...
- c 对某个整数做因式分解
1 #include <stdio.h> int main(void) { int n,i; scanf("%d",&n); printf("%d=& ...
- 调用AnimateWindow API来实现弹出效果
下面是实例的cs代码 public partial class frm_Main : Form { //使用Windows Api AnimateWindow [DllImport("use ...
- Informatica 9.5.1 安装配置
Informatica 结构 1个或多个资源库(Respository) PowerCenter数据整合引擎是基于元数据驱动的,提供了基于数据驱动的元数据知识库(Repository),该元数据知识 ...
- Swift和OC 混编
1.首先创建一个Swift工程 2.导入或者创建一个OC文件(.h和.m) 3.再创建一个桥连接文件 4.然后文件样子为 5.在桥接链接里面导入头文件 6.通过targets->->bui ...
- Web性能优化方案
第一章 打开网站慢现状分析 在公司访问部署在IDC机房的VIP网站时会感觉很慢.是什么原因造成的?为了缩短页面的响应时间,改进我们的用户体验,我们需要知道用户的时间花在等待什么东西上. 可以跟踪一下我 ...
- java.util.concurrent BlockingQueue
BlockingQueue 它实现了Queue接口.它是A BlockingQueue with one thread putting into it, and another thread taki ...
- QStringLiteral的两篇外文解释(编译期转换成QString)
http://blog.qt.io/blog/2014/06/13/qt-weekly-13-qstringliteral/ https://woboq.com/blog/qstringliteral ...