hdu 4612 Warm up
http://acm.hdu.edu.cn/showproblem.php?pid=4612
将原图进行缩点 变成一个树
树上每条边都是一个桥 然后加一条边要加在树的直径两端才最优
代码:
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<algorithm>
#include<queue>
#include<bitset>
#include<deque>
#include<numeric> #pragma comment(linker, "/STACK:1024000000,1024000000") using namespace std; typedef long long ll;
typedef unsigned int uint;
typedef pair<int,int> pp;
const double eps=1e-9;
const int INF=0x3f3f3f3f;
const ll MOD=1000000007;
const int N=210000;
const int M=2100000;
int head1[N],I1;
int head2[N],I2;
struct node
{
int j,next;
}edge1[M],edge2[M];
bool ev[M];
int dist[N];
int low[N],dfn[N],f[N],deep;
bool in[N],visited[N];
int num[N],sel[N];
stack<int>st;
queue<int>qt;
vector<int>vt[N];
void add1(int i,int j)
{
edge1[I1].j=j;
edge1[I1].next=head1[i];
head1[i]=I1++;
}
void add2(int i,int j)
{
edge2[I2].j=j;
edge2[I2].next=head2[i];
head2[i]=I2++;
}
void tarjan(int x)//将环缩点
{
visited[x]=true;
in[x]=true;
st.push(x);
low[x]=dfn[x]=deep++;
for(int t=head1[x];t!=-1;t=edge1[t].next)
if(!ev[t])
{
ev[t]=true;
ev[t^1]=true;
int j=edge1[t].j;
if(visited[j]==false)
{
tarjan(j);
low[x]=min(low[x],low[j]); }else if(in[j]==true)
{
low[x]=min(low[x],dfn[j]);
}
}
if(low[x]==dfn[x])
{
while(st.top()!=x)
{
int k=st.top(); st.pop();
in[k]=false;
f[k]=x;
vt[x].push_back(k);
}
int k=st.top(); st.pop();
in[k]=false;
f[k]=x;
}
}
void buildNewGraph(int k)
{
memset(head2,-1,sizeof(head2));I2=0;
for(int i=0;i<k;++i)
{
for(int t=head1[i];t!=-1;t=edge1[t].next)
{
int j=edge1[t].j;
if(f[i]!=f[j])
{
add2(f[j],f[i]);
++num[f[i]];
}
}
}
}
void init(int n,int m)
{
memset(head1,-1,sizeof(head1));
I1=0;
while(m--)
{
int l,r;
scanf("%d %d",&l,&r);
--l;--r;
add1(l,r);
add1(r,l);
}
}
int bfs(int x1,int n)
{
memset(dist,-1,sizeof(dist));
queue<int>qt;
dist[x1]=0;
qt.push(x1);
while(!qt.empty())
{
int x=qt.front();
qt.pop();
for(int t=head2[x];t!=-1;t=edge2[t].next)
{
int j=edge2[t].j;
if(dist[j]==-1)
{
dist[j]=dist[x]+1;
qt.push(j);
}
}
}
int k=0;
for(int i=0;i<n;++i)
if(dist[i]>dist[k])
k=i;
return k;
}
int solve(int n)
{
int s=0,sum=0;
for(int i=0;i<n;++i)
if(f[i]==i)
{s=i;++sum;}
int k=bfs(s,n);
k=bfs(k,n);
return (sum-1-dist[k]);
}
int main()
{
//freopen("data.in","r",stdin);
int n,m;
while(scanf("%d %d",&n,&m)!=EOF)
{
if(n==0&&m==0)
break;
init(n,m);
while(!st.empty()) st.pop();
for(int i=0;i<n;++i)
{vt[i].clear();f[i]=i;}
memset(in,false,sizeof(in));
memset(visited,false,sizeof(visited));
memset(ev,false,sizeof(ev));
deep=0;
for(int i=0;i<n;++i)
if(!visited[i])
tarjan(i);
buildNewGraph(n);
printf("%d\n",solve(n));
}
return 0;
}
hdu 4612 Warm up的更多相关文章
- Hdu 4612 Warm up (双连通分支+树的直径)
题目链接: Hdu 4612 Warm up 题目描述: 给一个无向连通图,问加上一条边后,桥的数目最少会有几个? 解题思路: 题目描述很清楚,题目也很裸,就是一眼看穿怎么做的,先求出来双连通分量,然 ...
- HDU 4612——Warm up——————【边双连通分量、树的直径】
Warm up Time Limit:5000MS Memory Limit:65535KB 64bit IO Format:%I64d & %I64u Submit Stat ...
- HDU 4612 Warm up tarjan 树的直径
Warm up 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=4612 Description N planets are connected by ...
- 【HDU 4612 Warm up】BCC 树的直径
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=4612 题意:一个包含n个节点m条边的无向连通图(无自环,可能有重边).求添加一条边后最少剩余的桥的数 ...
- hdu 4612 Warm up 双连通+树形dp思想
Warm up Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others) Total S ...
- HDU 4612 Warm up tarjan缩环+求最长链
Warm up Problem Description N planets are connected by M bidirectional channels that allow instant ...
- HDU 4612 Warm up(2013多校2 1002 双连通分量)
Warm up Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others)Total Su ...
- hdu 4612 Warm up 桥缩点
4612Warm hdu up 题目:给出一个图,添加一条边之后,问能够在新图中得到的最少的桥的数量. 分析:我们可以双联通分量进行缩点,原图变成了一棵树.问题变成了:求树中添加一条边之后,使得不在圈 ...
- hdu 4612 Warm up(缩点+树上最长链)
本来就是自己负责图论,结果水了= = 题目其实很裸,就是求桥的数量,只是要新加上一条边罢了.做法:先缩点.再在树上搜最长链(第一场多校的hdu 4607Park Visit就考了最长链,小样,套个马甲 ...
随机推荐
- zoj2589Circles(平面图的欧拉定理)
链接 连通图中: 设一个平面图形的顶点数为n,划分区域数为r,一笔画笔数为也就是边数m,则有: n+r-m=2 那么不算外面的那个大区域的话 就可以写为 n+r-m = 1 那么这个题就可以依次求出每 ...
- git tag之后如何修改
先 git clone 整个仓库,然后 git checkout tag_name 就可以取得 tag 对应的代码了. 但是这时候 git 可能会提示你当前处于一个“detached HEAD&quo ...
- C++ Primer 第二章 引用 指针 const限定符
1.引用: 为对象起了另外一个名字,引用类型引用另外一种类型,通过将声明符写成&d的形式来定义引用类型,其中d也就是声明的变量名(声明符就是变量名). PS:1.通过图片中编译所提示的报错信息 ...
- Sql Server中通配符的使用
Sql Server中通配符的使用 通配符_ "_"号表示任意单个字符,该符号只能匹配一个字符."_"可以放在查询条件的任意位置,且只能代表一个字符.一个汉字只 ...
- div滚动条演示
<!DOCTYPE html> <html> <head> <meta charset="gb2312" /> <title& ...
- 【bzoj3160】【xsy1726】万径人踪灭
[bzoj3160]万径人踪灭 题意 给定一个由'a'和'b'构成的字符串,求不连续回文子序列的个数. \(n\leq 100000\) 分析 还是蛮不错的. 这道题基本上是自己想到的. 除了没有利用 ...
- 能源项目xml文件 -- app-init.xml
<?xml version="1.0" encoding="UTF-8"?> <beans xmlns="http://www.sp ...
- python 练习 7
#!/usr/bin/python # -*- coding: utf-8 -*- def gcd(x,y): #最大公因子 if x>y:x%=y while x: x,y=y%x,x ret ...
- hdu----(2084)数塔(dp)
数塔 Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submissi ...
- ZigZag Conversion [LeetCode]
The string "PAYPALISHIRING" is written in a zigzag pattern on a given number of rows like ...