前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >2019 牛客暑期多校训练营 第五场 F maximum clique 1 最大独立集

2019 牛客暑期多校训练营 第五场 F maximum clique 1 最大独立集

作者头像
用户2965768
发布2019-08-14 14:33:32
3590
发布2019-08-14 14:33:32
举报
文章被收录于专栏:wym

题意:求n个点中bit至少有两位不同的数的最大个数并输出这些点

解:至少两位 补集是 至多一位不同也即恰好一位不同,因为这些数相异 ,不存在0为不同。

则原来求最大独立集 = 点数 - 最大匹配。 用dinc求最大匹配,按照顶层和非顶层点输出

代码语言:javascript
复制
#include <iostream>
#include <cstring>
#include <queue>
#include <stdio.h>
using namespace std;
#define INF 0x7ffffff
const int M = 1000005;
struct E{
    int nxt,v,w;
};
E edge[M];
int dpth[M],head[M],cnt;
int n,m,s,t;
void add(int u,int v,int w){
    cnt++;
    edge[cnt].v = v;
    edge[cnt].w = w;
    edge[cnt].nxt = head[u];
    head[u] = cnt;
    cnt++;
    edge[cnt].v = u;
    edge[cnt].w = 0;
    edge[cnt].nxt = head[v];
    head[v] = cnt;
}
//初始化从0开始存边 i^1是i的反向边
void init(){
    cnt=-1;
    memset(head,-1,sizeof(head));
}
//初始化分层图
bool bfs(){
    memset(dpth,0,sizeof(dpth));
    queue<int> q; q.push(s); dpth[s] = 1;
    while(!q.empty()){
        int u = q.front(); q.pop();
        for(int v,i=head[u];i!=-1;i=edge[i].nxt){
            {
                v = edge[i].v; 
                if(edge[i].w>0&&dpth[v]==0){
                    dpth[v] = dpth[u]+1;
                    q.push(v);
                }
            }
        }
    }
    if(dpth[t]!=0)
    return true;
    else return false;
}
int dfs(int u,int dist){
    if(u==t)
     return dist;
    for(int v,i=head[u];i!=-1;i=edge[i].nxt){
        v = edge[i].v;
        if((dpth[v]==dpth[u]+1)&&edge[i].w!=0)
        {
            int di = dfs(v,min(dist,edge[i].w));
            if(di>0){
                edge[i].w-=di;
                edge[i^1].w+=di;
                return di;
            }
        }
    }
    return 0;
}
int dinic(){
    int ans=0;
    while(bfs())
        while(int di = dfs(s,INF))
            ans+=di;
    return ans;
}
bool link[M];
int a[M];
int main()
{
    int n,m;
    scanf("%d",&n);
    init();
    s = n+1,t = n+2;
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        if(__builtin_popcount(a[i])&1)link[i] = true;
        if(link[i])add(s,i,1);
        else add(i,t,1);
    }
    for(int i=1;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            if(__builtin_popcount(a[i]^a[j])==1){
                if(link[i])add(i,j,1);
                else add(j,i,1);
            }
        }
    }
    int ans = n - dinic();
    printf("%d\n",ans);
    for(int i=1;i<=n;i++){
        if(dpth[i]!=0&&link[i])printf("%d ",a[i]);
        else if(dpth[i]==0&&!link[i])printf("%d ",a[i]);
    }
    return 0;
}
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2019年08月03日,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档