Minimal k-covering CodeForces - 976F

2024-02-22 18:58

本文主要是介绍Minimal k-covering CodeForces - 976F,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

点击打开链接

具体建图看官方题解点击打开链接

因为要求最小的边覆盖 所以要去掉最多的边 所以可以用网络流来解决

先从最大值到最小值来枚举 源点到二分图第一部分以及二分图第二部分到汇点所连的边的权值代表当前情况下每个点需要去掉几条边 然后跑完一遍最大流后残量网络中的边即为所求 然后加大源汇点的流量 表示下一次要多去掉一些边

#include <bits/stdc++.h>
using namespace std;
#define N 0x3f3f3f3fstruct node
{int v;int w;int id;int next;
};queue <int> que;
vector <int> ans[4010];
node edge[12010];
int first[4010],degree[4010],dis[4010],gap[4010],cur[4010],pre[4010];
int n1,n2,m,num,ss,ee;void addedge(int u,int v,int w,int id)
{edge[num].v=v;edge[num].w=w;edge[num].id=id;edge[num].next=first[u];first[u]=num++;return;
}void bfs()
{int i,u,v;while(!que.empty()) que.pop();memset(dis,-1,sizeof(dis));memset(gap,0,sizeof(gap));que.push(ee);dis[ee]=0;while(!que.empty()){u=que.front();que.pop();gap[dis[u]]++;for(i=first[u];i!=-1;i=edge[i].next){v=edge[i].v;if(dis[v]==-1){que.push(v);dis[v]=dis[u]+1;}}}return;
}int isap()
{int j,u,v,w,flow,minn,ans;bfs();memcpy(cur,first,sizeof(first));memset(pre,-1,sizeof(pre));u=ss,flow=N,ans=0;while(dis[ss]<num){int &i=cur[u];for(;i!=-1;i=edge[i].next){v=edge[i].v,w=edge[i].w;if(dis[v]+1==dis[u]&&w>0){pre[v]=i;u=v,flow=min(flow,w);if(u==ee){while(u!=ss){edge[pre[u]].w-=flow;edge[pre[u]^1].w+=flow;u=edge[pre[u]^1].v;}ans+=flow,flow=N;}break;}}if(i==-1){if(--gap[dis[u]]==0) break;cur[u]=first[u];minn=num-1;for(j=first[u];j!=-1;j=edge[j].next){v=edge[j].v,w=edge[j].w;if(w>0){minn=min(minn,dis[v]);}}dis[u]=minn+1;gap[dis[u]]++;if(u!=ss){u=edge[pre[u]^1].v;}}}return ans;
}int main()
{int gou[4010];int i,j,u,v,w,id,k;scanf("%d%d%d",&n1,&n2,&m);memset(first,-1,sizeof(first));memset(degree,0,sizeof(degree));num=0;for(i=1;i<=m;i++){scanf("%d%d",&u,&v);addedge(u,n1+v,1,i);addedge(n1+v,u,0,0);degree[u]++,degree[n1+v]++;}k=N;for(i=1;i<=n1+n2;i++){k=min(k,degree[i]);}ss=n1+n2+1,ee=n1+n2+2;for(i=1;i<=n1;i++){addedge(ss,i,degree[i]-k,0);addedge(i,ss,0,0);}for(i=n1+1;i<=n1+n2;i++){gou[i]=num;addedge(i,ee,degree[i]-k,0);addedge(ee,i,0,0);}num=n1+n2+2;for(i=k;i>=1;i--){isap();ans[i].clear();for(u=1;u<=n1;u++){for(j=first[u];j!=-1;j=edge[j].next){v=edge[j].v;if(v!=ss){w=edge[j].w,id=edge[j].id;if(w) ans[i].push_back(id);}}}for(j=first[ss];j!=-1;j=edge[j].next){edge[j].w++;}for(u=n1+1;u<=n1+n2;u++){edge[gou[u]].w++;}}for(i=0;i<=k;i++){printf("%d",ans[i].size());for(j=0;j<ans[i].size();j++){printf(" %d",ans[i][j]);}printf("\n");}return 0;
}

 

这篇关于Minimal k-covering CodeForces - 976F的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/736191

相关文章

Codeforces Round #240 (Div. 2) E分治算法探究1

Codeforces Round #240 (Div. 2) E  http://codeforces.com/contest/415/problem/E 2^n个数,每次操作将其分成2^q份,对于每一份内部的数进行翻转(逆序),每次操作完后输出操作后新序列的逆序对数。 图一:  划分子问题。 图二: 分而治之,=>  合并 。 图三: 回溯:

Codeforces Round #261 (Div. 2)小记

A  XX注意最后输出满足条件,我也不知道为什么写的这么长。 #define X first#define Y secondvector<pair<int , int> > a ;int can(pair<int , int> c){return -1000 <= c.X && c.X <= 1000&& -1000 <= c.Y && c.Y <= 1000 ;}int m

Codeforces Beta Round #47 C凸包 (最终写法)

题意慢慢看。 typedef long long LL ;int cmp(double x){if(fabs(x) < 1e-8) return 0 ;return x > 0 ? 1 : -1 ;}struct point{double x , y ;point(){}point(double _x , double _y):x(_x) , y(_y){}point op

Codeforces Round #113 (Div. 2) B 判断多边形是否在凸包内

题目点击打开链接 凸多边形A, 多边形B, 判断B是否严格在A内。  注意AB有重点 。  将A,B上的点合在一起求凸包,如果凸包上的点是B的某个点,则B肯定不在A内。 或者说B上的某点在凸包的边上则也说明B不严格在A里面。 这个处理有个巧妙的方法,只需在求凸包的时候, <=  改成< 也就是说凸包一条边上的所有点都重复点都记录在凸包里面了。 另外不能去重点。 int

Codeforces 482B 线段树

求是否存在这样的n个数; m次操作,每次操作就是三个数 l ,r,val          a[l] & a[l+1] &......&a[r] = val 就是区间l---r上的与的值为val 。 也就是意味着区间[L , R] 每个数要执行 | val 操作  最后判断  a[l] & a[l+1] &......&a[r] 是否= val import ja

Codeforces Round 971 (Div. 4) (A~G1)

A、B题太简单,不做解释 C 对于 x y 两个方向,每一个方向至少需要 x / k 向上取整的步数,取最大值。 由于 x 方向先移动,假如 x 方向需要的步数多于 y 方向的步数,那么最后 y 方向的那一步就不需要了,答案减 1 代码 #include <iostream>#include <algorithm>#include <vector>#include <string>

Minimal coverage -uva 覆盖线段,贪心

一道经典的贪心问题,具体方法就是将(an,bn)区间,按照an从小到大的顺序进行排序,之后从0开始, 取最大的有效区间,这里用到了结构体的快排,否则可能会超时. #include<stdio.h>#include<stdlib.h>#include<string.h>#define MAX_SIZE 100000 + 10#define BOTTOM -50000 - 10str

Codeforces#295(Div.2)A、B(模拟+BFS)

解题报告链接:点击打开链接 C. 题目链接:点击打开链接 解题思路: 对于给定的字符串,取出现次数最多的字母(可以同时有多个)。由这些字母组成长度为n的字符串,求有多少种组合。最后用数学知识即可。 完整代码: #include <algorithm>#include <iostream>#include <cstring>#include <climits>

Codeforces Round #281 (Div. 2)A(构造+暴力模拟)

题目链接:http://codeforces.com/problemset/problem/493/A 解题思路: 暴力的判断,分三种情况去判断即可。注意如果之前已经被罚下场后,那么在后面的罚下情况不应该算在输出结果内。 完整代码: #include <algorithm>#include <iostream>#include <cstring>#include <co

Codeforces Round #182 (Div. 2)A(水题)

题目链接:http://codeforces.com/contest/302/problem/A 解题思路: 只要通过重新排列使区间内和为0即是1,否则是0. 完整代码: #include <algorithm>#include <iostream>#include <cstring>#include <complex>#include <cstdio>#inc