题解:
真是一道神题!!!
大家还是围观JZP的题解吧(网址找不到了。。。)
代码:
#include<cstdio> #include<cstdlib> #include<cmath> #include<cstring> #include<algorithm> #include<iostream> #include<vector> #include<map> #include<set> #include<queue> #include<string> #define inf 1000000000 #define maxn 200000+5 #define maxm 200000+5 #define eps 1e-10 #define ll long long #define pa pair<int,int> #define for0(i,n) for(int i=0;i<=(n);i++) #define for1(i,n) for(int i=1;i<=(n);i++) #define for2(i,x,y) for(int i=(x);i<=(y);i++) #define for3(i,x,y) for(int i=(x);i>=(y);i--) #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go) #define for5(n,m) for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) #define mod 1000000007 using namespace std; inline int read() { int x=,f=;char ch=getchar(); while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();} while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();} return x*f; }
int n,m,s,t,cut[][],a[maxm],b[maxm],c[maxm],maxflow,tot=,head[maxn],cur[maxn],h[maxn]; queue<int>q;
bool v[maxn]; struct edge{int go,next,v;}e[maxm]; inline void add(int x,int y,int v) { e[++tot]=(edge){y,head[x],v};head[x]=tot; e[++tot]=(edge){x,head[y],};head[y]=tot; } bool bfs() { for(int i=;i<=n;i++)h[i]=-; q.push(s);h[s]=; while(!q.empty()) { int x=q.front();q.pop(); for(int i=head[x];i;i=e[i].next) if(e[i].v&&h[e[i].go]==-) { h[e[i].go]=h[x]+;q.push(e[i].go); } } return h[t]!=-; } int dfs(int x,int f) { if(x==t) return f; int tmp,used=; for(int i=cur[x];i;i=e[i].next) if(e[i].v&&h[e[i].go]==h[x]+) { tmp=dfs(e[i].go,min(e[i].v,f-used)); e[i].v-=tmp;if(e[i].v)cur[x]=i; e[i^].v+=tmp;used+=tmp; if(used==f)return f; } if(!used) h[x]=-; return used; }
inline void dfs(int x)
{
v[x]=;
for4(i,x)if(e[i].v&&!v[y])dfs(y);
} void dinic() {
maxflow=; while(bfs()) { for (int i=;i<=n;i++)cur[i]=head[i];maxflow+=dfs(s,inf); } }
inline void solve(int l,int r)
{
if(l==r)return;
s=b[l];t=b[r];
for1(i,tot)e[i].v=a[i];
dinic();
for1(i,n)v[i]=;
dfs(s);
for1(i,n)if(v[i])for1(j,n)if(!v[j])cut[i][j]=cut[j][i]=min(cut[i][j],maxflow);
int t1=l-,t2=r+;
for2(i,l,r)if(v[b[i]])c[++t1]=b[i];else c[--t2]=b[i];
for2(i,l,r)b[i]=c[i];
solve(l,t1);solve(t2,r);
} int main() {
int T=read();
while(T--)
{
n=read();m=read();
for1(i,n)head[i]=;tot=;
for1(i,m)
{
int x=read(),y=read(),w=read();
add(x,y,w);add(y,x,w);
}
for1(i,tot)a[i]=e[i].v;
for1(i,n)b[i]=i;
for1(i,n)for1(j,n)cut[i][j]=inf;
solve(,n);
int q=read();
while(q--)
{
int x=read(),ans=;
for1(i,n)for2(j,i+,n)if(cut[i][j]<=x)ans++;
printf("%d\n",ans);
}
printf("\n");
} return ; }
2229: [Zjoi2011]最小割
Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 685 Solved: 263
[Submit][Status]
Description
小
白在图论课上学到了一个新的概念——最小割,下课后小白在笔记本上写下了如下这段话:
“对于一个图,某个对图中结点的划分将图中所有结点分成两个部分,如果结点s,t不在同一个部分中,则称这个划分是关于s,t的割。
对于带权图来说,将所有顶点处在不同部分的边的权值相加所得到的值定义为这个割的容量,而s,t的最小割指的是在关于s,t的割中容量最小的割”
现给定一张无向图,小白有若干个形如“图中有多少对点它们的最小割的容量不超过x呢”的疑问,小蓝虽然很想回答这些问题,但小蓝最近忙着挖木块,于是作为
仍然是小蓝的好友,你又有任务了。
Input
输入文件第一行有
且只有一个正整数T,表示测试数据的组数。 对于每组测试数据, 第一行包含两个整数n,m,表示图的点数和边数。
下面m行,每行3个正整数u,v,c(1<=u,v<=n,0<=c<=106),表示有一条权为c的无向边(u,v)
接下来一行,包含一个整数q,表示询问的个数 下面q行,每行一个整数x,其含义同题目描述。
Output
对于每组测试数据,输出应包括q行,第i行表示第i个问题的答案。对于点对(p,q)和(q,p),只统计一次(见样例)。
两组测试数据之间用空行隔开。
Sample Input
1
5 0
1
0
5 0
1
0
Sample Output
10
【数据范围】
对于100%的数据 T<=10,n<=150,m<=3000,q<=30,x在32位有符号整数类型范围内。
图中两个点之间可能有多条边