关闭

BZOJ 2115 Xor(dfs&独立回路&异或消元)

标签: 异或消元数学
321人阅读 评论(0) 收藏 举报
分类:

学习的秦总的写法写的dfs找独立回路和异或消元:秦总博客

独立回路的定义和独立回路个数=MN+1的证明,博客里面都有,但是dfs找的过程好像需要2倍的空间,不能理解。。。

分析:这题要明白一个东西,u>v的所有路径可以由u>v的一条任意路径+该图的所有独立回路的线性组合完成,其实想一想还是很好理解的,有一点很重要:比如1>7是一条路,然后怎么进过回路2>3,3>4,4>2,因为整个图是连通的,所以17点必然有一个点在一个回路里面,这样这个回路走两遍就不影响异或值了,这个时候回路上面的所有点又可以走一个回路,所以如果一直这样拓展一定能够走到2>3,3>4,4>2这条回路,并且其他回路都自己抵消了。

附上代码:

#include <bits/stdc++.h>
#define LL long long
#define FOR(i,x,y)  for(int i = x;i < y;++ i)
#define IFOR(i,x,y) for(int i = x;i > y;-- i)

using namespace std;

const int maxn = 50050;
const int maxm = 100010;

int head[maxn],edge_cnt;

struct Edge{
    int u,v,nt;
    LL val;
}edge[maxm<<1];

void add_edge(int u,int v,LL val){
    edge[edge_cnt].u = u;
    edge[edge_cnt].v = v;
    edge[edge_cnt].val = val;
    edge[edge_cnt].nt = head[u];
    head[u] = edge_cnt++;
}

int n,m,tot;

LL val[maxm<<1];
LL d[maxn];
bool vis[maxn];

void Build_Graph(){
    memset(vis,false,sizeof(vis));
    edge_cnt = 0;
    memset(head,-1,sizeof(head));
    int u,v;
    LL val;
    FOR(i,0,m){
        scanf("%d%d%lld",&u,&v,&val);
        add_edge(u,v,val);
        add_edge(v,u,val);
    }
}

void dfs(int u,int fa){
    vis[u] = true;
    for(int i = head[u];i != -1;i = edge[i].nt){
        int v = edge[i].v;
        if(v == fa) continue;
        if(!vis[v]){
            d[v] = d[u] ^ edge[i].val;
            dfs(v,u);
        }
        else{
            val[tot++] = d[u]^d[v]^edge[i].val;
        }
    }
}

int xorguass(){
    int row = 0;
    for(int i = 62;i >= 0;-- i){
        int j;
        for(j = row;j < tot;++ j) {
            if(val[j] & (1LL<<i)) break;
        }
        if(j != tot){
            swap(val[j],val[row]);
            for(j = 0;j < tot;++ j){
                if(j == row)    continue;
                if(val[j] & (1LL<<i)){
                    val[j] ^= val[row];
                }
            }
            ++ row;
        }
    }
    return row;
}

void work(){
    LL ans = d[n];
    n = xorguass();
    FOR(i,0,n)  ans = max(ans,ans^val[i]);
    printf("%lld\n",ans);
}
int main()
{
    //freopen("test.in","r",stdin);
    while(~scanf("%d%d",&n,&m)){
        Build_Graph();
        tot = 0;
        d[1] = 0;
        dfs(1,-1);
        work();
    }
    return 0;
}

0
0

猜你在找
【直播】机器学习&数据挖掘7周实训--韦玮
【套餐】系统集成项目管理工程师顺利通关--徐朋
【直播】3小时掌握Docker最佳实战-徐西宁
【套餐】机器学习系列套餐(算法+实战)--唐宇迪
【直播】计算机视觉原理及实战--屈教授
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之矩阵--黄博士
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之凸优化--马博士
【套餐】Javascript 设计模式实战--曾亮
查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:52835次
    • 积分:1968
    • 等级:
    • 排名:第19858名
    • 原创:150篇
    • 转载:0篇
    • 译文:0篇
    • 评论:8条
    最新评论