二分圖

二分圖判定:

題目鏈接:二分圖判定

dfs:

const int __=10005;

bool flag;
int col[__];
vector<int>G[__];

void dfs(int x,int y)
{
    if(col[x])
    {
        if(col[x]!=y)flag=false;
        return;
    }
    col[x]=y;
    for(int i=0; i<G[x].size() && flag; i++)
        dfs(G[x][i],3-y);
}

void init(void)
{
    memset(col,0,sizeof(col));
    flag=true;
}

int main()
{
    int _;
    scanf("%d",&_);
    while(_--)
    {
        init();
        int n,m;
        scanf("%d%d",&n,&m);
        for(int i=1; i<=m; i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            G[x].push_back(y);
            G[y].push_back(x);
        }
        for(int i=1; i<=n; i++)
            if(!col[i])dfs(i,1);
        if(flag)printf("Correct\n");
        else printf("Wrong\n");
        for(int i=1;i<=n;i++)G[i].clear();
    }
    return 0;
}

最大匹配:

題目鏈接:最大匹配-匈牙利算法

dfs:

vector<int>G[1005];
bool vis[1005];
int m,n;
int match[1005];

bool dfs(int x)
{
    int len=G[x].size();
    for(int i=0; i<len; i++)
    {
        int nex=G[x][i];
        if(!vis[nex])
        {
            vis[nex]=true;
            if(!match[nex]||dfs(match[nex]))
            {
                match[nex]=x;
                match[x]=nex;
                return true;
            }
        }
    }
    return false;
}

int main()
{
    int x,y,ans=0;
    scanf("%d%d",&n,&m);
    for(int i=0;i<m;i++)
    {
        scanf("%d%d",&x,&y);
        G[x].push_back(y);
        G[y].push_back(x);
    }
    for(int i=1;i<=n;i++)
    {
        memset(vis,false,sizeof(vis));
        if(!match[i]&&dfs(i))ans++;
    }
    printf("%d\n",ans);
    for(int i=1;i<=n;i++)G[i].clear();
    return 0;
}

二維最大匹配:

題目鏈接:Antenna Placement

struct point
{
    int x,y;
};
char maze[50][50];
bool vis[50][50];
int m,n;
vector<point>vec[50][50];
point match[50][50];

void judgepoint(int x,int y)
{
    point t;
    if(x>1&&maze[x-1][y]=='*')
    {
        t.x=x-1;
        t.y=y;
        vec[x][y].push_back(t);
    }
    if(x<m&&maze[x+1][y]=='*')
    {
        t.x=x+1;
        t.y=y;
        vec[x][y].push_back(t);
    }
    if(y>1&&maze[x][y-1]=='*')
    {
        t.x=x;
        t.y=y-1;
        vec[x][y].push_back(t);
    }
    if(y<n&&maze[x][y+1]=='*')
    {
        t.x=x;
        t.y=y+1;
        vec[x][y].push_back(t);
    }
}

bool dfs(int x,int y)
{
    int len=vec[x][y].size();
    for(int i=0; i<len; i++)
    {
        point temp=vec[x][y][i];
        if(!vis[temp.x][temp.y])
        {
            vis[temp.x][temp.y]=1;
            if((!match[temp.x][temp.y].x&&!match[temp.x][temp.y].y)
                    ||dfs(match[temp.x][temp.y].x,match[temp.x][temp.y].y))
            {
                match[temp.x][temp.y].x=x;
                match[temp.x][temp.y].y=y;
                match[x][y].x=temp.x;
                match[x][y].y=temp.y;
                return true;
            }
        }
    }
    return false;
}

void init(void)
{
    memset(maze,0,sizeof(maze));
    memset(match,0,sizeof(match));
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        init();
        scanf("%d%d",&m,&n);
        for(int i=1; i<=m; i++)
            scanf("%s",maze[i]+1);
        for(int i=1; i<=m; i++)
            for(int j=1; j<=n; j++)
                if(maze[i][j]=='*')
                    judgepoint(i,j);
        int ans=0;
        for(int i=1; i<=m; i++)
            for(int j=1; j<=n; j++)
                if(maze[i][j]=='*'&&!match[i][j].x&&!match[i][j].y)
                {
                    memset(vis,0,sizeof(vis));
                    if(dfs(i,j))ans++;
                }

        for(int i=1; i<=m; i++)
            for(int j=1; j<=n; j++)
                if(maze[i][j]=='*'&&!match[i][j].x&&!match[i][j].y)
                    ans++;
        printf("%d\n",ans);
        for(int i=1; i<=m; i++)
            for(int j=1; j<=n; j++)
                if(maze[i][j]=='*')vec[i][j].clear();
    }
    return 0;
}

題目鏈接:最小點覆蓋和最大獨立集

最小點覆蓋=最大匹配數

最大獨立集=點數-最大匹配數

vector<int>G[1005];
bool vis[1005];
int m,n;
int match[1005];

bool dfs(int x)
{
    int len=G[x].size();
    for(int i=0; i<len; i++)
    {
        int nex=G[x][i];
        if(!vis[nex])
        {
            vis[nex]=true;
            if(!match[nex]||dfs(match[nex]))
            {
                match[nex]=x;
                match[x]=nex;
                return true;
            }
        }
    }
    return false;
}

int main()
{
    int x,y,ans=0;
    scanf("%d%d",&n,&m);
    for(int i=0;i<m;i++)
    {
        scanf("%d%d",&x,&y);
        G[x].push_back(y);
        G[y].push_back(x);
    }
    for(int i=1;i<=n;i++)
    {
        memset(vis,false,sizeof(vis));
        if(!match[i]&&dfs(i))ans++;
    }
    printf("%d\n%d\n",ans,n-ans);
    for(int i=1;i<=n;i++)G[i].clear();
    return 0;
}

題目鏈接:最小路徑覆蓋

最小路徑覆蓋=點數-最大匹配數

vector<int>G[1005];
bool vis[1005];
int match[1005];

bool dfs(int x)
{
    int len=G[x].size();
    for(int i=0; i<len; i++)
    {
        int nex=G[x][i];
        if(!vis[nex])
        {
            vis[nex]=true;
            if(!match[nex]||dfs(match[nex]))
            {
                match[nex]=x;
                match[x]=nex;
                return true;
            }
        }
    }
    return false;
}

int main()
{
    int n,m,x,y,ans=0;
    scanf("%d%d",&n,&m);
    for(int i=0;i<m;i++)
    {
        scanf("%d%d",&x,&y);
        G[x].push_back(n+y);
    }
    for(int i=1;i<=n;i++)
    {
        memset(vis,false,sizeof(vis));
        if(!match[i]&&dfs(i))ans++;
    }
    printf("%d\n",n-ans);
    return 0;
}

題目鏈接:帶權二分匹配

km:

int cost[305][305],lx[305],ly[305];
int match[305],d;
bool visx[305],visy[305];

bool dfs(int x,int ny)
{
    visx[x]=true;
    for(int y=1; y<=ny; y++)
    {
        if(visy[y])continue;
        int t=lx[x]+ly[y]-cost[x][y];
        if(!t)
        {
            visy[y]=true;
            if(!match[y] || dfs(match[y],ny))
                return match[y]=x;
        }
        else d=min(d,t);
    }
    return false;
}

int km(int nx,int ny)
{
    memset(lx,0,sizeof(lx));
    memset(ly,0,sizeof(ly));
    memset(match,0,sizeof(match));
    for(int i=1; i<=nx; i++)
        for(int j=1; j<=ny; j++)
            lx[i]=max(lx[i],cost[i][j]);
    for(int x=1; x<=nx; x++)
        while(1)
        {
            d=inf;
            memset(visx,false,sizeof(visx));
            memset(visy,false,sizeof(visy));
            if(dfs(x,ny))break;
            if(d==inf)return -1;
            for(int x=1; x<=nx; x++)
                if(visx[x])lx[x]-=d;
            for(int y=1; y<=ny; y++)
                if(visy[y])ly[y]+=d;
        }
    int res=0;
    for(int y=1; y<=ny; y++)
        if(match[y])
            res+=cost[match[y]][y];
    return res;
}


int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        for(int i=1; i<=n; i++)
            for(int j=1; j<=n; j++)
                scanf("%d",&cost[i][j]);
        printf("%d\n",km(n,n));
    }
    return 0;
}
最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。

推薦閱讀更多精彩內容

  • 說說二分圖,其實圖論的題難點不在用算法,難在如何建圖,只有圖建好了,剩下的就簡單了,在這說說求二分圖的算法,即匈牙...
    Anxdada閱讀 333評論 0 0
  • 二分圖又稱作二部圖,是圖論中的一種特殊模型。 設G=(V,E)是一個無向圖,如果頂點V可分割為兩個互不相交的子集(...
    Myth52125閱讀 10,307評論 0 4
  • 定義## 記圖 G = (V, E) 最大流問題##### 記每條邊所能傳輸的最大流量為 c(e) 記每條邊所能傳...
    MrGopher閱讀 985評論 0 1
  • 題目描述 給定一個二分圖,結點個數分別為n,m,邊數為e,求二分圖最大匹配數 輸入輸出格式 輸入格式第一行,n,m...
    Ricardo_Y_Li閱讀 335評論 0 1
  • 離結婚的期限還有兩個月零三天,但結婚照還沒拍,婚前任何準備的東西都還沒有準備,婚期越近我自己越焦慮,也許大家都...
    命運的弄潮兒閱讀 279評論 0 0