题目链接

Description

给定一个有 \(n\) 个点和 \(m\) 条边的无向图,其中每一条边 \(e_i\) 都有一个权值记为 \(w_i\) 。

对于给出的两个点 \(a\) 和 \(b\) ,求一条 \(a\) 到 \(b\) 的路径,使得路径上的边权的 \(OR\)(位或)和最小,输出这个值。(也就是说,如果将路径看做边的集合 \(\{e_1,e_2,…,e_k\}\),那么这条路径的代价为 \(w_1\ OR\ w_2\ OR\ …\ OR\ w_k\),现在求一条路径使得其代价最小,输出这个代价。如果不存在这样的路径,输出 \(-1\)。

Input

Easy

\(2\le n \le 10^4, 0 \le m \le 10^6, 0 \le c_i \le 2^{62}-1\)

Hard

\(2\le n \le 10^4, 0 \le m \le 10^6, 0 \le c_i \le 2^{62}-1\)

题解 By zerol

Easy

\(dp[u][k]\) 表示到结点 \(u\) 的代价为 \(k\) 的方案是否存在,然后在图上转移(dfs 一下就好了)。

Hard

假设答案的二进制位全是 \(1\),然后从高位到低位考虑,如果将该位置为 \(0\) 不破坏连通性的话就置为 \(0\),这样肯定最优。

判断连通性可以用 并查集 或者 搜索,反正 \(O(M)\) 就行。

复杂度 \(O(63M)\) 。

对比

稍微回顾一下之前的 bzoj 2115 [Wc2011] Xor 路径最大异或和 线性基

共同点:从高位向低位做,判断当前位能否置为1或0(毕竟都是位运算)

不同点:

Xor那道题是预处理出来一条路径,以及所有可以补充于其上的路径。所以,判断能否置为1即是看能否添加这条路径。

而这道题,因为OR运算的性质,直观想来,添加的路径越少OR和就越小。所以,判断当前位能否置为0即是通过连通性来判断。

Code

Easy

#include <bits/stdc++.h>
#define maxn 1100
#define maxm 10010
using namespace std;
bool vis[maxn][maxn];
struct Edge { int to, ne, w; } edge[maxm << 1];
int tot, ne[maxn];
void add(int u, int v, int w) {
edge[tot] = {v, ne[u], w};
ne[u] = tot++;
}
typedef long long LL;
void dfs(int u, int ors) {
if (vis[u][ors]) return;
vis[u][ors] = true;
for (int i = ne[u]; ~i; i = edge[i].ne) {
int v = edge[i].to;
dfs(v, ors | edge[i].w);
}
}
int main() {
memset(ne, -1, sizeof ne);
int n, m;
scanf("%d%d", &n,&m);
int u, v, w;
for (int i = 0; i < m; ++i) {
scanf("%d%d%d", &u,&v,&w);
add(u,v,w); add(v, u, w);
}
scanf("%d%d", &u, &v);
dfs(u, 0);
for (int i = 0; i <= 1024; ++i) {
if (vis[v][i]) { printf("%d\n", i); return 0; }
}
puts("-1");
return 0;
}

Hard

#include <bits/stdc++.h>
#define maxn 10010
#define maxm 1000010
using namespace std;
typedef long long LL;
vector<int> a[maxn];
struct Edge { int u, v; } edge[maxm];
void add(LL w, int id) {
int cnt = 0;
while (w) {
if (w & 1) a[cnt].push_back(id);
w >>= 1, ++cnt;
}
}
int s, t, n, m, fa[maxn], sz[maxn];
bool exc[maxm], flag[64];
int find(int x) { return fa[x] == x ? x : fa[x] = find(fa[x]); }
void unionn(int u, int v) {
u = find(u), v = find(v);
if (sz[u] < sz[v]) swap(u, v);
sz[u] += sz[v], fa[v] = u;
}
bool ok(int id) {
memset(exc, 0, sizeof exc);
for (int i = 62; i >= id; --i) {
if (!flag[i]) {
for (auto x : a[i]) exc[x] = true;
}
}
for (int i = 1; i <= n; ++i) fa[i] = i, sz[i] = 1;
for (int i = 0; i < m; ++i) {
if (!exc[i]) unionn(edge[i].u, edge[i].v);
}
return find(s) == find(t);
}
int main() {
scanf("%d%d", &n,&m);
for (int i = 0; i < m; ++i) {
int u, v; LL w;
scanf("%d%d%lld", &u,&v,&w);
edge[i] = {u, v};
add(w, i);
}
scanf("%d%d", &s, &t);
if (!ok(64)) { puts("-1"); return 0; }
LL ans = 0;
for (int i = 62; i >= 0; --i) {
ans <<= 1;
if (!ok(i)) flag[i] = 1, ans |= 1;
}
printf("%lld\n", ans);
return 0;
}

EOJ Monthly 2018.1 F 最小OR路径的更多相关文章

  1. EOJ Monthly 2018.8 D. Delivery Service-树上差分(边权/边覆盖)(边权转点权)(模板题)

    D. Delivery Service 单测试点时限: 2.5 秒 内存限制: 512 MB EOJ Delivery Service Company handles a massive amount ...

  2. EOJ Monthly 2018.2

    A. 坑爹的售票机 题意 用\(1,5,10,25,50,100\)的纸币买\(n\)张单价为\(p\)的船票,且一次性最多买\(k\)张,求钱数恰好时最少需要多少张纸币. Hard: \(n,k,p ...

  3. EOJ Monthly 2018.1

    985月赛,当时鸽了,现在想补一补 A. 石头剪刀布的套路 Time limit per test: 1.0 seconds Memory limit: 256 megabytes 现在有一种石头剪刀 ...

  4. EOJ Monthly 2018.11 D. 猜价格

    猜价格 分两种情况讨论: k≤n,先猜至多 k 次 1,由于回答 <1 肯定是假的,所以可以把剩余系下是哪次错试出来,然后用至多 n 次搞定. k>n,每个数都猜两次,如果两次结果不一样, ...

  5. 【EOJ Monthly 2018.7】【D数蝌蚪】

    https://acm.ecnu.edu.cn/contest/92/problem/D/ D. 数蝌蚪 Time limit per test: 2.0 seconds Memory limit:  ...

  6. EOJ Monthly 2018.7

    准备继续大学acm啦 又要开始愉快的码码码啦 第一次在华东师大OJ上面做题 看来EOJ上的积分体质是假的,我怎么一把上红??? A.数三角形 神tm的防AK题放在A,出题人很不友好啊... 先写了个暴 ...

  7. EOJ Monthly 2018.3

    985月赛我只喜欢ECNU.jpg A. 打工时不可能打工的 Time limit per test: 2.0 seconds Memory limit: 256 megabytes 我 Ayano ...

  8. EOJ Monthly 2018.4

    A. ultmaster 的小迷妹们 Time limit per test: 2.0 seconds Memory limit: 256 megabytes ultmaster 男神和他的小迷妹们准 ...

  9. 【EOJ Monthly 2018.2 (Good bye 2017)】

    23333333333333333 由于情人节要回家,所以就先只放代码了. 此题是与我胖虎过不去. [E. 出老千的 xjj] #include<cstdio> #include<c ...

随机推荐

  1. 10分钟了解 react 引入的 Hooks

    "大家好,我是谷阿莫,今天要将的是一个...",哈哈哈,看到这个题我就想到这个开头.最近react 官方在 2018 ReactConf 大会上宣布 React v16.7.0-a ...

  2. myeclipse10.5 crack(2012-12-27-bd 写的日志迁移

    首先去网上下一个破解文件如图: 解压过后打开的文件夹如图: 再打开crack文件夹如图: 运行run.bat如果点击它没反应就是你没有安装jdk,它如果运行就如图所示: 到这一步就在第一个方框user ...

  3. 致敬wusir懒孩子自有懒孩子的生存之道之二

    https://www.cnblogs.com/wupeiqi/ https://www.cnblogs.com/Eva-J/ https://www.cnblogs.com/wupeiqi/p/90 ...

  4. 04vim的使用

    linux常用命令 workon 查看已经安装的虚拟环境 deactivate 退出虚拟环境 whoami 查看用户 sudo bash install.sh 添加权限 pwd 查看在那个路径下 cd ...

  5. 用python给图片添加半透明水印

    # coding:utf-8 from PIL import Image, ImageDraw, ImageFont def add_text_to_image(image, text): font ...

  6. Find a path HDU - 5492 (dp)

    Find a path Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total ...

  7. 【树状数组】CF961E Tufurama

    挺巧妙的数据结构题(不过据说这是一种套路? E. Tufurama One day Polycarp decided to rewatch his absolute favourite episode ...

  8. MySQL查询数据库中表名或字段名

    查询数据库中所有表名 select table_name from information_schema.tables where table_schema='csdb' and table_type ...

  9. Java面向对象---抽象类与接口

    final关键字 1.final关键字在java中被称为完结器,表示最终的意思: 2.final能声明类.方法.属性: 3.使用final声明的类不能被继承,使用final声明的方法不能被重写,使用f ...

  10. tensorboard在cmd运行成功但在浏览器中不能正常显示的问题解决

     我是配置了两个python环境,python3.5和anconda3.5,强烈建议使用python3.5版本,算是比较稳定的! cmd在运行时是按顺序查找的文件,如果说是python3.6这个版本问 ...