http://www.lydsy.com/JudgeOnline/problem.php?id=4519

最小割树

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm> using namespace std; #define N 900
#define M 9000
const int inf=2e9; int n; int tot=;
int front[N],nxt[M<<],to[M<<],val[M<<],from[M<<];
int lev[N],num[N];
int path[N];
int cur[N]; int src,decc; int a[N],tmp[N];
bool use[N]; int dis[N][N]; int ans[]; void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
} void add(int u,int v,int w)
{
to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; from[tot]=u; val[tot]=w;
} bool bfs()
{
queue<int>q;
for(int i=;i<=n;++i) lev[i]=n;
q.push(decc);
lev[decc]=;
int now,t;
while(!q.empty())
{
now=q.front();
q.pop();
for(int i=front[now];i;i=nxt[i])
{
t=to[i];
if(lev[t]==n && val[i^])
{
lev[t]=lev[now]+;
q.push(t);
}
}
}
return lev[src]!=n;
} int augment()
{
int now=decc,flow=inf;
int i;
while(now!=src)
{
i=path[now];
flow=min(flow,val[i]);
now=from[i];
}
now=decc;
while(now!=src)
{
i=path[now];
val[i]-=flow;
val[i^]+=flow;
now=from[i];
}
return flow;
} int isap()
{
int flow=;
if(!bfs()) return ;
memset(num,,sizeof(num));
for(int i=;i<=n;++i) num[lev[i]]++,cur[i]=front[i];
int now=src,t;
while(lev[src]<n)
{
if(now==decc)
{
flow+=augment();
now=src;
}
bool advanced=false;
for(int i=cur[now];i;i=nxt[i])
{
t=to[i];
if(lev[t]==lev[now]- && val[i])
{
advanced=true;
path[t]=i;
cur[now]=i;
now=t;
break;
}
}
if(!advanced)
{
int mi=n;
for(int i=front[now];i;i=nxt[i])
if(val[i]) mi=min(mi,lev[to[i]]);
if(!--num[lev[now]]) break;
num[lev[now]=mi+]++;
cur[now]=front[now];
if(now!=src) now=from[path[now]];
}
}
return flow;
} void dfs(int x)
{
use[x]=true;
for(int i=front[x];i;i=nxt[i])
if(!use[to[i]] && val[i]) dfs(to[i]);
} void solve(int l,int r)
{
if(l==r) return;
for(int i=;i<=tot;i+=) val[i]=val[i+]=val[i]+val[i+]>>;
src=a[l]; decc=a[r];
int flow=isap();
memset(use,false,sizeof(use));
dfs(src);
for(int i=;i<=n;++i)
if(use[i])
for(int j=;j<=n;++j)
if(!use[j])
dis[i][j]=dis[j][i]=min(dis[i][j],flow);
int i=l,j=r;
for(int k=l;k<=r;++k)
if(use[a[k]]) tmp[i++]=a[k];
else tmp[j--]=a[k];
for(int k=l;k<=r;++k) a[k]=tmp[k];
solve(l,i-);
solve(j+,r);
} int main()
{
int m;
int u,v,w;
memset(dis,,sizeof(dis));
read(n); read(m);
for(int i=;i<=n;++i) a[i]=i;
while(m--)
{
read(u); read(v); read(w);
add(u,v,w);
add(v,u,w);
}
solve(,n);
int sum=;
for(int i=;i<=n;++i)
for(int j=i+;j<=n;++j)
ans[++sum]=dis[i][j];
sort(ans+,ans+sum+);
sum=unique(ans+,ans+sum+)-ans-;
cout<<sum;
return ;
}

4519: [Cqoi2016]不同的最小割

Time Limit: 20 Sec  Memory Limit: 512 MB
Submit: 866  Solved: 499
[Submit][Status][Discuss]

Description

学过图论的同学都知道最小割的概念:对于一个图,某个对图中结点的划分将图中所有结点分成
两个部分,如果结点s,t不在同一个部分中,则称这个划分是关于s,t的割。对于带权图来说,将
所有顶点处在不同部分的边的权值相加所得到的值定义为这个割的容量,而s,t的最小割指的是在
关于s,t的割中容量最小的割。
而对冲刺NOI竞赛的选手而言,求带权图中两点的最小割已经不是什么难事了。我们可以把
视野放宽,考虑有N个点的无向连通图中所有点对的最小割的容量,共能得到N(N−1)
2个数值。
这些数值中互不相同的有多少个呢?这似乎是个有趣的问题。
 

Input

输入文件第一行包含两个数N,M,表示点数和边数。接下来M行,每行三个数u,v,w,
表示点u和点v(从1开始标号)之间有条边权值是w。
1<=N<=850 1<=M<=8500 1<=W<=100000
 

Output

输出文件第一行为一个整数,表示个数。

 

Sample Input

4 4
1 2 3
1 3 6
2 4 5
3 4 4

Sample Output

3

最新文章

  1. 前端开发之走进Vue.js
  2. jbuilder的set!方法重构接口
  3. 扩展easyui 的表单验证 (转)
  4. DataTable添加行和列数据
  5. 如果设置http.get超时控制
  6. VLC命令行的应用
  7. 表单提交中get与post的区别
  8. css 问题总结
  9. 一种最坏情况线性运行时间的选择算法 - The missing worst-case linear-time Select algorithm in CLRS.
  10. 3)Javascript设计模式:Observer模式
  11. VS2013 F12无法转到函数的定义处,总是从“元数据”获取的问题 ——解决方法
  12. 用 CSS3 做一个流星雨动画
  13. 【BZOJ2152】聪聪可可(点分治)
  14. 网站开发进阶(三十四)编码中的setCharacterEncoding 理解
  15. java中利用dom4j解析XML文件
  16. java 中的迭代
  17. eclipse c++11 cmake gnuradio
  18. CY7C68013 USB接口相机开发记录 - 第一天:资料下载
  19. 设计一个BCD码计数器。
  20. [转]VB 读写ini 配置文件

热门文章

  1. 【CS231N】1、图像分类
  2. java异常处理及自定义异常的使用
  3. Alpha版本冲刺(九)
  4. [并查集] 1107. Social Clusters (30)
  5. Keil C51 的printf
  6. js 算法
  7. 安恒杯2月月赛-应该不是xss
  8. UFLDL学习笔记 ---- 主成分分析与白化
  9. 微信小程序与java后台交互
  10. P4071 [SDOI2016]排列计数