传送门


IOI强行交互可还行,我Luogu的代码要改很多才能交到UOJ去……

发现问题是对边权做限制的连通块类问题,考虑\(Kruskal\)重构树进行解决。

对于图上的边\((u,v)(u<v)\),我们建两棵\(Kruskal\)重构树,一棵按照\(u\)从大到小加边表示人形时的活动区域,一棵按照\(v\)从小到大加边表示狼形时的活动区域。

那么我们的每组询问就变为了:给出两段区间,询问是否存在一个点同时覆盖这两个区间(即人形转换为狼形的地点)。这个是经典的二维数点问题,使用主席树计算即可。

话说感觉一堆namespace的码风看起来好舒服啊~

#include<bits/stdc++.h>
//This code is written by Itst
using namespace std;

inline int read(){
    int a = 0;
    char c = getchar();
    bool f = 0;
    while(!isdigit(c)){
        if(c == '-')
            f = 1;
        c = getchar();
    }
    while(isdigit(c)){
        a = (a << 3) + (a << 1) + (c ^ '0');
        c = getchar();
    }
    return f ? -a : a;
}

const int MAXN = 2e5 + 10;
int N , M , Q;
namespace ST{
#define mid ((l + r) >> 1)
    struct node{
        int l , r , sum;
    }Tree[MAXN * 50];
    int root[MAXN] , cnt;

    void insert(int& x , int p , int l , int r , int tar){
        x = ++cnt;
        Tree[x] = Tree[p];
        ++Tree[x].sum;
        if(l == r)
            return;
        if(mid >= tar)
            insert(Tree[x].l , Tree[p].l , l , mid , tar);
        else
            insert(Tree[x].r , Tree[p].r , mid + 1 , r , tar);
    }

    bool query(int x , int y , int l , int r , int L , int R){
        if(Tree[x].sum == Tree[y].sum)
            return 0;
        if(l >= L && r <= R)
            return 1;
        bool f = 0;
        if(mid >= L)
            f |= query(Tree[x].l , Tree[y].l , l , mid , L , R);
        if(!f && mid < R)
            f |= query(Tree[x].r , Tree[y].r , mid + 1 , r , L , R);
        return f;
    }
}

namespace Itst{
#define P pair < int , int >
    int fa[2][MAXN << 1] , ch[2][MAXN << 1][2] , jump[2][MAXN << 1][21] , que[2][MAXN << 1] , dfn[2][MAXN] , cnt[2] , ts[2] , rg[2][MAXN << 1][2];

    void init(){
        cnt[0] = cnt[1] = N;
        for(int i = 1 ; i <= N ; ++i)
            fa[0][i] = fa[1][i] = i;
    }

    int find(int ind , int x){
        return fa[ind][x] == x ? x : (fa[ind][x] = find(ind , fa[ind][x]));
    }

    inline void link(int ind , int x , int y){
        x = find(ind , x);
        y = find(ind , y);
        if(x == y)
            return;
        ++cnt[ind];
        fa[ind][x] = fa[ind][y] = fa[ind][cnt[ind]] = cnt[ind];
        ch[ind][cnt[ind]][0] = x;
        ch[ind][cnt[ind]][1] = y;
    }

    void dfs(int ind , int x , int p){
        jump[ind][x][0] = p;
        for(int i = 1 ; jump[ind][x][i - 1] ; ++i)
            jump[ind][x][i] = jump[ind][jump[ind][x][i - 1]][i - 1];
        if(x <= N){
            rg[ind][x][0] = rg[ind][x][1] = dfn[ind][x] = ++ts[ind];
            que[ind][x] = x;
            return;
        }
        dfs(ind , ch[ind][x][0] , x);
        dfs(ind , ch[ind][x][1] , x);
        que[ind][x] = ind ? min(que[ind][ch[ind][x][0]] , que[ind][ch[ind][x][1]]) : max(que[ind][ch[ind][x][0]] , que[ind][ch[ind][x][1]]);
        rg[ind][x][0] = rg[ind][ch[ind][x][0]][0];
        rg[ind][x][1] = rg[ind][ch[ind][x][1]][1];
    }

    int up(int ind , int x , int l , int r){
        for(int i = 19 ; i >= 0 ; --i)
            if(jump[ind][x][i] && que[ind][jump[ind][x][i]] >= l && que[ind][jump[ind][x][i]] <= r)
                x = jump[ind][x][i];
        return x;
    }

    void work(){
        dfs(0 , cnt[0] , 0);
        dfs(1 , cnt[1] , 0);
        vector < P > p;
        for(int i = 1 ; i <= N ; ++i)
            p.push_back(P(dfn[0][i] , dfn[1][i]));
        sort(p.begin() , p.end());
        for(int i = 1 ; i <= N ; ++i)
            ST::insert(ST::root[i] , ST::root[i - 1] , 1 , N , p[i - 1].second);

        while(Q--){
            int S = read() + 1 , E = read() + 1 , L = read() + 1 , R = read() + 1;
            E = up(0 , E , 1 , R);
            S = up(1 , S , L , N);
            puts(ST::query(ST::root[rg[0][E][1]] , ST::root[rg[0][E][0] - 1] , 1 , N , rg[1][S][0] , rg[1][S][1]) ? "1" : "0");
        }
    }
}

namespace init{
    struct Edge{
        int a , b;
    }Ed[MAXN << 1];

    bool cmp1(Edge a , Edge b){
        return a.b < b.b;
    }

    bool cmp2(Edge a , Edge b){
        return a.a > b.a;
    }

    void main(){
        N = read();
        M = read();
        Q = read();
        Itst::init();
        for(int i = 1 ; i <= M ; ++i){
            Ed[i].a = read() + 1;
            Ed[i].b = read() + 1;
            if(Ed[i].a > Ed[i].b)
                swap(Ed[i].a , Ed[i].b);
        }
        sort(Ed + 1 , Ed + M + 1 , cmp1);
        for(int i = 1 ; i <= M ; ++i)
            Itst::link(0 , Ed[i].a , Ed[i].b);
        sort(Ed + 1 , Ed + M + 1 , cmp2);
        for(int i = 1 ; i <= M ; ++i)
            Itst::link(1 , Ed[i].a , Ed[i].b);
        Itst::work();
    }
}

int main(){
#ifndef ONLINE_JUDGE
    freopen("in" , "r" , stdin);
    //freopen("out" , "w" , stdout);
#endif
    init::main();
    return 0;
}

Luogu4899 IOI2018 Werewolf 主席树、Kruskal重构树的更多相关文章

  1. Kruskal重构树入门

    这个知识点好像咕咕咕了好长了..趁还没退役赶紧补一下吧.. 讲的非常简略,十分抱歉.. 前置知识 Kruskal算法 一定的数据结构基础(如主席树) Kruskal重构树 直接bb好像不是很好讲,那就 ...

  2. [IOI2018] werewolf 狼人 kruskal重构树,主席树

    [IOI2018] werewolf 狼人 LG传送门 kruskal重构树好题. 日常安利博客文章 这题需要搞两棵重构树出来,这两棵重构树和我们平时见过的重构树有点不同(据说叫做点权重构树?),根据 ...

  3. [IOI2018] werewolf 狼人 [kruskal重构树+主席树]

    题意: 当你是人形的时候你只能走 \([L,N-1]\) 的编号的点(即大于等于L的点) 当你是狼形的时候你只能走 \([1,R]\) 的编号的点(即小于等于R的点) 然后问题转化成人形和狼形能到的点 ...

  4. UOJ#407. 【IOI2018】狼人 Kruskal,kruskal重构树,主席树

    原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ407.html 题解 套路啊. 先按照两个节点顺序各搞一个kruskal重构树,然后问题转化成两棵krus ...

  5. LOJ.2865.[IOI2018]狼人(Kruskal重构树 主席树)

    LOJ 洛谷 这题不就是Peaks(加强版)或者归程么..这算是\(IOI2018\)撞上\(NOI2018\)的题了? \(Kruskal\)重构树(具体是所有点按从小到大/从大到小的顺序,依次加入 ...

  6. 【BZOJ-3545&3551】Peaks&加强版 Kruskal重构树 + 主席树 + DFS序 + 倍增

    3545: [ONTAK2010]Peaks Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1202  Solved: 321[Submit][Sta ...

  7. BZOJ 3551: [ONTAK2010]Peaks加强版 [Kruskal重构树 dfs序 主席树]

    3551: [ONTAK2010]Peaks加强版 题意:带权图,多组询问与一个点通过边权\(\le lim\)的边连通的点中点权k大值,强制在线 PoPoQQQ大爷题解传送门 说一下感受: 容易发现 ...

  8. bzoj 3551 kruskal重构树dfs序上的主席树

    强制在线 kruskal重构树,每两点间的最大边权即为其lca的点权. 倍增找,dfs序对应区间搞主席树 #include<cstdio> #include<cstring> ...

  9. [IOI2018]狼人——kruskal重构树+可持久化线段树

    题目链接: IOI2018werewolf 题目大意:给出一张$n$个点$m$条边的无向图,点和边可重复经过,一个狼人初始为人形,有$q$次询问,每次询问要求人形态只能处于编号不小于$L$的点,狼形态 ...

随机推荐

  1. SD从零开始62-63,不完全日志,业务伙伴及业务伙伴确定

    [原创] SD从零开始62 不完全日志 不完全日志Incompletion log 一个不完全日志是销售凭证中对你公司重要的而还没有在系统中输入的所有数据的清单: 你可以在配置中为不完全日志定义这些数 ...

  2. Linux 下修改网卡MAC地址

    Linux下修改网卡MAC地址 by:授客 QQ:1033553122 例子:修改网卡接口eth0的mac地址 #停用网卡接口,比如eth0 # ifconfig eth0 down #编辑对应的网卡 ...

  3. loadrunner 场景设计-目标场景设计

    场景设计-目标场景设计 by:授客 QQ:1033553122 A.   概述 Goals Types for Goal-Oriented Scenarios 在以目标为向导的场景中,定义你想实现的测 ...

  4. loadrunner 运行脚本-命令行运行脚本

    Loadrunner 运行脚本-命令行运行脚本   by:授客 QQ:1033553122 脚本所在目录 Run-time Settings->Additional Attributes设置   ...

  5. spring容器、BeanFactory、ApplicatContext、WebApplicationContext

    1.spring容器 BeanFactory提供了IoC的功能:ApplicatContext建立在BeanFactory的基础上 在初始化BeanFactory的实现类时,并没有实例化Bean,直到 ...

  6. python常用模块json

    python jons模块 json模块 主要是解决数据格式的转换问题,比如python接收到json对象需要转换为python对象,供python处理,亦或者python数据需要发送到其给其他客户端 ...

  7. Java ——基础语法

    package myhello; // 本类所在的包的路径 import af.util.AfMath; // 导入对应的类 import java.util.Random; // 导入随机数的类 p ...

  8. MySQL表名区分大小写设置

    关闭MySQL服务 在服务运行目录找到my.ini或者my.cnf文件 打开文件,找到[mysqld]在下面增加一行 lower_case_table_names=0 (0:大小写敏感;1:大小写不敏 ...

  9. 理解LSTM

    本文基于Understanding-LSTMs进行概括整理,对LSTM进行一个简单的介绍 什么是LSTM LSTM(Long Short Term Memory networks)可以解决传统RNN的 ...

  10. centos7升级内核版本

    本文转载http://blog.csdn.net/nciasd/article/details/51490146,大神非常厉害!!!!! 查看当前系统的内核版本 # uname -r 1.导入key ...