二分図最小点オーバーライドセットテンプレート(C++版)
18060 ワード
最小カバー数を求めるだけなら、直接最大マッチングを求めるのがいいです.また、ポイントセットを上書きする必要がある場合は、テンプレートを写します.ここは古典的な最小カバーポイントセットのテーマです.だから、左右のポイントセットは確定しています.前のテンプレートと最大マッチングは違います.
#include
using namespace std;
const int N=100050;
const int M=100050;
int n,m;
char s[N];
struct Edge{
int v,next;
}edge[M];
int cnt,head[N];
void init(){
cnt=0;
memset(head,-1,sizeof(head));
}
void addEdge(int u,int v){
edge[cnt]=Edge{v,head[u]};
head[u]=cnt++;
}
//Left
int Left[N],Right[N];
//S T
bool S[N],T[N];
//X Y
vector<int> X,Y;
bool dfs(int u){
S[u]=true;
for(int i=head[u];i!=-1;i=edge[i].next){
int v=edge[i].v;
if(!T[v]){
T[v]=true;
if(Left[v]==-1 || dfs(Left[v])){
Left[v]=u;
Right[u]=v;
return true;
}
}
}
return false;
}
void solve(){
memset(S,0,sizeof(S));
memset(T,0,sizeof(T));
//
for(int i=1;i<=n;i++){
if(Right[i]==-1){
dfs(i);
}
}
//
for(int u=1;u<=n;u++){
if(!S[u]){
X.push_back(u);
}
}
//
for(int v=1;v<=m;v++){
if(T[v]){
Y.push_back(v);
}
}
}
int main(void){
scanf("%d%d",&n,&m);
init();
for(int i=1;i<=n;i++){
scanf("%s",s+1);
for(int j=1;j<=m;j++){
if(s[j]=='*'){
addEdge(i,j);
}
}
}
solve();
// ,
printf("%d",X.size());
for(int i=0;i<X.size();i++){
printf(" %d",X[i]);
}
printf("
");
printf("%d",Y.size());
for(int i=0;i<Y.size();i++){
printf(" %d",Y[i]);
}
printf("
");
return 0;
}