• Evanyou Blog 彩带


      先上题目传送门

      最小路径覆盖其实是一类题目,一般用二分图匹配或者网络流都可以做。相关定理:最小路径覆盖数=顶点数-最大割(即最大匹配)

      这道题的要求除了求出最小路径覆盖数之外还要求输出每一条路径,这里蒟蒻用的是最简单粗暴的深搜,从每一个路径覆盖的起点开始一个个搜索直到输出所有点。

      奉上用Dinic算法做的正解:

    #include<cstdio>
    #include<cstring>
    #include<cstdlib>
    #include<cmath>
    #include<queue>
    #include<iostream>
    #include<algorithm>
    #define inf 1e9
    using namespace std;
    const int N=10010;
    int n,m,dis[N<<1],ans;
    int head[N<<1],size=1;
    int fa[N<<1],start,endd;
    struct Node{
      int from,to,val,next;
    }edge[N<<1];
    queue<int>team;
    inline int read()
    {
      char ch=getchar();int num=0;bool flag=false;
      while(ch<'0'||ch>'9'){if(ch=='-')flag=true;ch=getchar();}
      while(ch>='0'&&ch<='9'){num=num*10+ch-'0';ch=getchar();}
      return flag?-num:num;
    }
    inline void add(int x,int y,int z)
    {
      edge[++size].from=x;
      edge[size].to=y;
      edge[size].val=z;
      edge[size].next=head[x];
      head[x]=size;
    }
    inline void add_edge(int x,int y)
    {add(x,y,1);add(y,x,0);}
    inline bool bfs()
    {
      memset(dis,-1,sizeof(dis));
      while(!team.empty())team.pop();
      team.push(start);dis[start]=0;
      while(!team.empty()){
        int x=team.front();team.pop();
        for(int i=head[x];i!=-1;i=edge[i].next){
          int y=edge[i].to;
          if(edge[i].val>0&&dis[y]==-1){
        dis[y]=dis[x]+1;team.push(y);
          }
        }
      }
      return dis[endd]>-1;
    }
    inline int dinic(int u,int flow)
    {
      if(u==endd||flow==0)return flow;
      int now=0;
      for(int i=head[u];i!=-1;i=edge[i].next){
        int v=edge[i].to;
        if(edge[i].val>0&&dis[v]==dis[u]+1){
          int ka=dinic(v,min(flow-now,edge[i].val));
          edge[i].val-=ka;
          edge[i^1].val+=ka;
          now+=ka;
          if(now==flow)break;
        }
      }
      return now;
    }
    inline int find(int x)
    {return fa[x]==x?x:fa[x]=find(fa[x]);}
    inline void print(int u)
    {
      printf("%d ",u);
      for(int i=head[u];i!=-1;i=edge[i].next){
        int v=edge[i].to;
        if(edge[i].val==0&&v>n)
          print(v-n);
      }
    }
    void ready()
    {
      memset(head,-1,sizeof(head));
      n=read();m=read();start=0;endd=n<<1|1;
      for(int i=1;i<=m;i++){
        int x=read();int y=read();
        add_edge(x,y+n);}
      for(int i=1;i<=n;i++){
        add_edge(start,i);
        add_edge(i+n,endd);
      }
    }
    void work()
    {
      while(bfs())
        ans+=dinic(start,inf);
      for(int i=1;i<=n;i++)fa[i]=i;
      for(int i=2;i<=size;i++){
        int u=edge[i].from,v=edge[i].to;
        if(u>start&&u<=n&&v>n&&v<endd&&edge[i].val==0)
          fa[find(v-n)]=find(u);}
        for(int i=1;i<=n;i++)
          if(find(i)==i)
        print(i),printf("
    ");
        printf("%d",n-ans);
    }
    int main()
    {
      ready();
      work();
      return 0;
    }

      这里再上一份用二分图匹配算法做的最小路径覆盖模板(不输出路径):

    //It is made by HolseLee on 30th Dec 2017
    #include<cstdio>
    #include<cstring>
    #include<cstdlib>
    #include<cmath>
    #include<iostream>
    #include<algorithm>
    #include<vector>
    using namespace std;
    const int N=10010;
    int n,m,ans,ma[N];
    bool vis[N];
    vector<int>team[N];
    inline int read()
    {
      char ch=getchar();int num=0;bool flag=false;
      while(ch<'0'||ch>'9'){if(ch=='-')flag=true;ch=getchar();}
      while(ch>='0'&&ch<='9'){num=num*10+ch-'0';ch=getchar();}
      return flag?-num:num;
    }
    inline bool dfs(int u)
    {
      for(int i=0;i<team[u].size();i++){
        int v=team[u][i];
        if(vis[v])continue;
        vis[v]=true;
        if(ma[v]==-1||dfs(ma[v])){
          ma[v]=u;
          return true;}
      }
      return false;
    }
    void ready()
    {
      memset(ma,-1,sizeof(ma));
      n=read();m=read();
      for(int i=1;i<=n;i++)
        team[i].clear();
      for(int i=1;i<=m;i++){
        int x=read();int y=read();
        team[x].push_back(y);
      }
    }
    void work()
    {
      int ret=0;
      for(int i=1;i<=n;i++){
        memset(vis,false,sizeof(vis));
        ret+=dfs(i);
      }
      ans=n-ret;
      printf("%d",ans);
    }
    int main()
    {
      ready();
      work();
      return 0;
    }
  • 相关阅读:
    java 调用可执行文件时,ProcessBuilder异常CreateProcess error=2
    easyUI行内编辑与jdbc批量更新
    Oracle中merge into应用举例
    mybatis插入List<Map<String, String>>批量数据到Oracle数据库
    插入数据库值大于数据库字段设置的长度导致的mapper执行中断,控制台未报错
    ORA-25156:旧样式的外部联接(+)不能与 ANSI链接一起使用
    Oracle列转行函数listagg使用演示
    Oracle查询序列和函数方法封装
    js取url问号后的参数方法封装
    linux查看文件具体时间和大小
  • 原文地址:https://www.cnblogs.com/cytus/p/8149032.html
Copyright © 2020-2023  润新知