loj3161「NOI2019」I 君的探险(随机化,整体二分)

loj

Luogu

题解时间

对于 $ N \le 500 $ 的点,毫无疑问可以直接 $ O(n^2) $ 暴力询问解决。

考虑看起来最好做的 $ B $ 类。

由于有每个点的父亲编号小于该点的优良特性,很容易想到整体二分。

考虑用整体二分求出每个点的父亲:

对于一个分治区间,毫无疑问 $ [l,mid] $ 的节点的父亲在左区间。

而对于另外一半节点,考虑将左半节点全部modify,此时右半某个节点亮起则说明左半节点至少有一个与之相邻,则其父亲在左半区间内。

而对于其他数据,我们依然可以大胆尝试类似的做法,只不过要维护已经确定的边的影响,

并且注意每次将序列随机打乱,同时已经加完边的节点要从序列中删除。

#include<bits/stdc++.h>
// #include"explore.h"
using namespace std;
typedef long long lint;
struct pat{int x,y;pat(int x=0,int y=0):x(x),y(y){}bool operator<(const pat &p)const{return x==p.x?y<p.y:x<p.x;}};
void modify(int x);
int query(int x);
void report(int x,int y);
int check(int x);
const int N=400011;
int n,m;
namespace task500
{
bool check(){return n<=500;}
int tag[N];
void explore()
{
for(int i=0;i<n-1;i++){modify(i);for(int j=i+1;j<n;j++)if(query(j)!=tag[j]) tag[j]^=1,report(i,j);}
}
}
vector<int> ve[N];int tcnt;
namespace taskB
{
bool check(){return n==99997||n==199997;}
void fuck(int l,int r,int px)
{
if(l==r){for(auto x:ve[px])if(x!=l) report(l,x);return;}
int mid=l+r>>1,ls=++tcnt,rs=++tcnt;
for(int i=l;i<=mid;i++) modify(i);
for(auto x:ve[px])
if(x<=mid||query(x)) ve[ls].push_back(x);
else ve[rs].push_back(x);
for(int i=l;i<=mid;i++) modify(i);
fuck(l,mid,ls),fuck(mid+1,r,rs);
}
void explore()
{
tcnt++;for(int i=0;i<n;i++) ve[tcnt].push_back(i);
fuck(0,n-1,tcnt);
}
}
namespace tasknormal
{
int cnt;
int a[N],on[N];
vector<int> to[N];
void report(int x,int y){::report(x,y);to[x].push_back(y),to[y].push_back(x),cnt++;}
int query(int x){int ret=::query(x);for(auto y:to[x]) ret^=on[y];return ret;}
void fuck(int l,int r,int px)
{
if(l==r){for(auto x:ve[px])if(x!=l) report(a[l],a[x]);ve[px].clear();return;}
int mid=l+r>>1,ls=++tcnt,rs=++tcnt;
for(int i=l;i<=mid;i++) modify(a[i]),on[a[i]]=1;
for(auto x:ve[px])
if(x<=mid||query(a[x])) ve[ls].push_back(x);
else ve[rs].push_back(x);
for(int i=l;i<=mid;i++) modify(a[i]),on[a[i]]=0;
fuck(l,mid,ls),fuck(mid+1,r,rs);
ve[px].clear();
}
void explore()
{
for(int i=0;i<n;i++) a[i]=i;
while(cnt<m)
{
random_shuffle(a,a+n);
tcnt=1;for(int i=0;i<n;i++) ve[tcnt].push_back(i);
fuck(0,n-1,tcnt);if(cnt<m)for(int i=0;i<n;i++)if(check(a[i])){swap(a[i],a[n-1]),n--,i--;}
}
}
}
void explore(int n,int m)
{
::n=n,::m=m;
if(task500::check()) task500::explore();
else if(taskB::check()) taskB::explore();
else tasknormal::explore();
}

loj3161「NOI2019」I 君的探险(随机化,整体二分)的更多相关文章

  1. 【LOJ】#2985. 「WC2019」I 君的商店

    LOJ#2985. 「WC2019」I 君的商店 一道很神仙的题啊QAQ 居然是智商题--不是乱搞或者是大数据 我们可以用2N问出一个最大值是1 然后对于任意两个值\(x + y\)和\(a\)比较 ...

  2. loj2985「WC2019」I 君的商店(二分,思维)

    loj2985「WC2019」I 君的商店(二分,思维) loj Luogu 题解时间 真的有点猛的思维题. 首先有一个十分简单的思路: 花费 $ 2N $ 确定一个为 $ 1 $ 的数. 之后每次随 ...

  3. luogu P5473 [NOI2019]I 君的探险 交互 随机 二分 分治 整体二分

    LINK:I 君的探险 神仙题! 考虑一个暴力的做法 每次点亮一个点 询问全部点 这样询问次数为 \(\frac{n\cdot (n-1)}{2}\) 可以通过前5个点. 考虑都为A的部分分 发现一个 ...

  4. LOJ #2985. 「WC2019」I 君的商店

    传送门 搬题解QwQ 首先最大值一定为 \(1\),直接扫一遍两两比较 \(O(2N)\) 求出最大值 设最大值位置为 \(a\),对于任意两个没有确定的位置 \(x,y\) 询问 \([a,x+y] ...

  5. LOJ 3158: 「NOI2019」序列

    题目传送门:LOJ #3158. 题意简述: 给定两个长度为 \(n\) 的正整数序列 \(a,b\),要求在每个序列中都选中 \(K\) 个下标,并且要保证同时在两个序列中都被选中的下标至少有 \( ...

  6. LOJ 3160: 「NOI2019」斗主地

    题目传送门:LOJ #3160. 简要题意: 有一个长度为 \(n\) 的序列 \(a\),初始时 \(a_i=i\) 或 \(a_i=i^2\),这取决于 \(\mathrm{type}\) 的值. ...

  7. LOJ 3159: 「NOI2019」弹跳

    题目传送门:LOJ #3159. 题意简述: 二维平面上有 \(n\) 个整点,给定每个整点的坐标 \((x_i,y_i)\). 有 \(m\) 种边,第 \(i\) 种边从 \(p_i\) 号点连向 ...

  8. LOJ 3156: 「NOI2019」回家路线

    题目传送门:LOJ #3156. 题意简述: 有一张 \(n\) 个点 \(m\) 条边的有向图,边有两个权值 \(p_i\) 和 \(q_i\)(\(p_i<q_i\))表示若 \(p_i\) ...

  9. 「NOI2019」弹跳(KD树)

    题意:w×h网格中有n个点,m条边.每条边可以从p点花费t时间到一个矩形中的任意点,求1号点到每个点的最少时间. \(1<=w,h<=n<=70000,1<=m<=150 ...

随机推荐

  1. Solution -「JOISC 2021」「LOJ #3491」道路建设

    \(\mathcal{Description}\)   Link.   平面上有 \(n\) 个互不重合的点 \((x_{1..n},y_{1..n})\),求其两两曼哈顿距离的前 \(m\) 小值. ...

  2. design PLL

    PLL PLL主要是NCO中通过相位比较器,控制压控振荡器,控制频率.

  3. c++ 拷贝构造函数、拷贝运算符、析构函数

    拷贝构造函数.拷贝运算符.析构函数 拷贝构造函数.拷贝运算符.析构函数 定义行为像值的类 class HasPtr{ public: HasPtr(const string &s = stri ...

  4. Python基础—装饰器(Day11)

    装饰器 1.装饰器是在不改变原函数的执行的情况下为原函数增额外的功能. 简单版装饰器import time def func1(): print('执行速度') def timmer(f): star ...

  5. NeurIPS 2017 | QSGD: Communication-Efficient SGD via Gradient Quantization and Encoding

    由于良好的可扩展性,随机梯度下降(SGD)的并行实现是最近研究的热点.实现并行化SGD的关键障碍就是节点间梯度更新时的高带宽开销.因此,研究者们提出了一些启发式的梯度压缩方法,使得节点间只传输压缩后的 ...

  6. 关于 Word2Vec 使用时遇到的一系列问题!!

    ​ 1  训练时 model = Word2Vec(x, size=250, window=5, min_count=5, workers=12, iter=10, sg=1) 这句代码一直报错 查了 ...

  7. Oracle账户被锁:the account is locked

    转至:https://blog.csdn.net/weixin_37615080/article/details/80400239?utm_medium=distribute.pc_relevant_ ...

  8. 60天shell脚本计划-2/12-渐入佳境

    --作者:飞翔的小胖猪 --创建时间:2021年2月1日 --修改时间:2021年2月5日 说明 每日上传更新一个shell脚本,周期为60天.如有需求的读者可根据自己实际情况选用合适的脚本,也可在评 ...

  9. ssh执行命令或者脚本里最好先加上 source /etc/profile

    记: 写一次shell脚本,发现总提示 找不到JAVA_HOME .明明已经配了环境变量,并且java -version也已经能输出信息,但就是提示. 后来在 脚本中加入第一句代码 source /e ...

  10. 一、ES6基础

    一.ECMAScript和JavaScript关系 JavaScript 的创造者 Netscape 公司,决定将 JavaScript 提交给标准化组织 ECMA,希望这种语言能够成为国际标 准,但 ...