原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ191.html

题目传送门 - UOJ191

题意

自行移步集训队论文2016中罗哲正的论文。

题解

自行移步集训队论文2016中罗哲正的论文。

代码

所以说我就是来存代码的?

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
LL read(){
LL x=0,f=0;
char ch=getchar();
while (!isdigit(ch))
f|=ch=='-',ch=getchar();
while (isdigit(ch))
x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
return f?-x:x;
}
const int N=300015,mod=998244353;
int Type,n,m,Q;
vector <int> e[N];
struct Gragh{
static const int M=N;
int cnt,y[M],nxt[M],fst[N];
void clear(){
cnt=1;
memset(fst,0,sizeof fst);
}
void add(int a,int b){
y[++cnt]=b,nxt[cnt]=fst[a],fst[a]=cnt;
}
}qs;
struct Point{
int x,y;
Point(){}
Point(int _x,int _y){
x=_x,y=_y;
}
Point operator - (Point a){
return Point(x-a.x,y-a.y);
}
};
LL Cross(Point a,Point b){
return (LL)a.x*b.y-(LL)b.x*a.y;
}
struct Query{
Point a;//a for cross
int f,x,tmp;
LL ans;
}q[N];
int depth[N],fa[N];
Point P[N];
int new_node(int f,int d){
n++;
e[n].clear(),fa[n]=f,depth[n]=d;
e[f].push_back(n);
e[n].push_back(f);
return n;
}
void build_tree(){
static int st[N],d;
st[d=Q=0]=n=1;
e[n].clear(),fa[n]=0,depth[n]=d;
while (m--){
int opt=read(),L,R,x,y;
if (opt==1){
x=read(),y=read();
d++;
st[d]=new_node(st[d-1],d);
P[n]=Point(x,y);
}
else if (opt==2)
d--;
else if (opt==3){
L=read(),R=read(),x=read(),y=read();
Q++;
q[Q].a=Point(x,y);
q[Q].f=st[L],q[Q].x=st[R];
q[Q].ans=-5e18;
}
}
qs.clear();
for (int i=1;i<=Q;i++)
qs.add(q[i].x,i);
}
int Size,RT;
int size[N],Max[N],vis[N];
void get_root(int x,int pre){
size[x]=1,Max[x]=0;
for (auto y : e[x])
if (y!=pre&&!vis[y]){
get_root(y,x);
size[x]+=size[y];
Max[x]=max(Max[x],size[y]);
}
Max[x]=max(Max[x],Size-size[x]);
if (Max[x]<Max[RT])
RT=x;
}
void get_size(int x,int pre){
size[x]=1;
for (auto y : e[x])
if (y!=pre&&!vis[y])
get_root(y,x),size[x]+=size[y];
}
Point A[N];
int qid[N],qtot;
void get_query(int x){
for (int id=qs.fst[x];id;id=qs.nxt[id])
if (depth[q[qs.y[id]].f]<=depth[RT])
qid[++qtot]=qs.y[id];
for (auto y : e[x])
if (y!=fa[x]&&!vis[y])
get_query(y);
}
bool cmp(int a,int b){
return q[a].tmp<q[b].tmp;
}
bool cmp2(int a,int b){
return Cross(q[a].a,q[b].a)<0;
}
Point tmp[N];
int Tmp[N];
void Unknown(int *qid,int qt,Point *P,int &pt){
if (pt<=1)
return (void)sort(qid+1,qid+qt+1,cmp2);
int mid=pt>>1,ls=mid,rs=pt-mid,kk=0;
while (kk<qt&&q[qid[kk+1]].tmp<mid)
kk++;
for (int i=kk+1;i<=qt;i++)
q[qid[i]].tmp-=mid;
Unknown(qid,kk,P,ls);
Unknown(qid+kk,qt-kk,P+mid,rs);
for (int i=kk+1,j=1;i<=qt;i++){
while (j<ls&&Cross(P[j+1]-P[j],q[qid[i]].a)<0)
j++;
q[qid[i]].ans=max(q[qid[i]].ans,Cross(q[qid[i]].a,P[j]));
}
int tot=0,j=1;
for (int i=1;i<=ls;){
Point now=j<=rs&&P[j+mid].x<=P[i].x?P[mid+(j++)]:P[i++];
while (tot>1&&Cross(tmp[tot]-tmp[tot-1],now-tmp[tot-1])>=0)
tot--;
tmp[++tot]=now;
}
for (;j<=rs;j++){
Point now=P[j+mid];
while (tot>1&&Cross(tmp[tot]-tmp[tot-1],now-tmp[tot-1])>=0)
tot--;
tmp[++tot]=now;
}
while (tot>1&&tmp[tot].x==tmp[tot-1].x&&tmp[tot].y<=tmp[tot-1].y)
tot--;
int h=1;
while (h<tot&&tmp[h].x==tmp[h+1].x&&tmp[h].y<=tmp[h+1].y)
h++;
for (int i=h;i<=tot;i++)
P[i-h+1]=tmp[i];
pt=tot-h+1;
tot=0,j=kk+1;
for (int i=1;i<=kk;)
if (j<=qt&&Cross(q[qid[j]].a,q[qid[i]].a)<0)
Tmp[++tot]=qid[j++];
else
Tmp[++tot]=qid[i++];
for (;j<=qt;j++)
Tmp[++tot]=qid[j];
for (int i=1;i<=qt;i++)
qid[i]=Tmp[i];
}
void solve(int x){
Max[RT=0]=N,get_root(x,0),x=RT,get_size(x,0);
int t=0;
for (int i=x;!vis[i];i=fa[i])
A[++t]=P[i];
vis[x]=1,RT=x,qtot=0,get_query(x);
for (int i=1;i<=qtot;i++)
q[qid[i]].tmp=min(t,depth[x]-depth[q[qid[i]].f]+1);
sort(qid+1,qid+qtot+1,cmp);
int k=1;
while (k<=qtot&&q[qid[k]].tmp<t)
k++;
Unknown(qid,k-1,A,t);
sort(qid+k,qid+qtot+1,cmp2);
for (int i=k,j=1;i<=qtot;i++){
while (j<t&&Cross(A[j+1]-A[j],q[qid[i]].a)<0)
j++;
q[qid[i]].ans=max(q[qid[i]].ans,Cross(q[qid[i]].a,A[j]));
}
for (auto y : e[x])
if (!vis[y])
Size=size[y],solve(y);
}
int main(){
Type=read();
while (m=read()){
build_tree(),Size=n;
memset(vis,0,sizeof vis);
vis[0]=1,solve(1);
int ans=0;
for (int i=1;i<=Q;i++)
ans^=(int)((q[i].ans%mod+mod)%mod);
cout << ans << endl;
}
return 0;
}

  

UOJ#191. 【集训队互测2016】Unknown 点分治 分治 整体二分 凸包 计算几何的更多相关文章

  1. UOJ#191. 【集训队互测2016】Unknown

    题意:维护一个数列,每个元素是个二维向量,每次可以在后面加一个元素或者删除一个元素.给定P(x,y),询问对于[l,r]区间内的元素$S_i$,$S_i \times P$的最大值是多少. 首先简单地 ...

  2. uoj #190. 【集训队互测2016】消失的源代码 提交答案题

    Test 1: 发现是一个字母表的映射 把 \('a' \to 'z'\) 打进去找出映射就好了QAQ . Test 2: 求助 \(dalao\) 得知的点.. 答案是 : \(2016x^2 + ...

  3. 【loj2461】【2018集训队互测Day 1】完美的队列

    #2461. 「2018 集训队互测 Day 1」完美的队列 传送门: https://loj.ac/problem/2461 题解: 直接做可能一次操作加入队列同时会弹出很多数字,无法维护:一个操作 ...

  4. 【2018集训队互测】【XSY3372】取石子

    题目来源:2018集训队互测 Round17 T2 题意: 题解: 显然我是不可能想出来的……但是觉得这题题解太神了就来搬(chao)一下……Orzpyz! 显然不会无解…… 为了方便计算石子个数,在 ...

  5. 洛谷 P4463 - [集训队互测 2012] calc(多项式)

    题面传送门 & 加强版题面传送门 竟然能独立做出 jxd 互测的题(及其加强版),震撼震撼(((故写题解以祭之 首先由于 \(a_1,a_2,\cdots,a_n\) 互不相同,故可以考虑求出 ...

  6. 【uoj#94】【集训队互测2015】胡策的统计(集合幂级数)

    题目传送门:http://uoj.ac/problem/94 这是一道集合幂级数的入门题目.我们先考虑求出每个点集的连通生成子图个数,记为$g_S$,再记$h_S$为点集$S$的生成子图个数,容易发现 ...

  7. 【集训队互测2015】Robot

    题目描述 http://uoj.ac/problem/88 题解 维护两颗线段树,维护最大值和最小值,因为每次只有单点查询,所以可以直接在区间插入线段就可以了. 注意卡常,不要写STL,用链表把同类修 ...

  8. LOJ3069. 「2019 集训队互测 Day 1」整点计数(min_25筛)

    题目链接 https://loj.ac/problem/3069 题解 复数真神奇. 一句话题意:令 \(f(x)\) 表示以原点 \((0, 0)\) 为圆心,半径为 \(x\) 的圆上的整点数量, ...

  9. BZOJ3938 & UOJ88:[集训队互测2015]Robot——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=3938 http://uoj.ac/problem/88 小q有n只机器人,一开始他把机器人放在了一 ...

随机推荐

  1. thyemleaf:禁用JS缓存(原创)

    在开发时经常需要调整JS,但是调整后由于页面缓存的原因,看不到实时效果. 开发环境:springboot+thymeleaf 1.配置文件多模式 2.获得当前的激活的模式和随机数 import org ...

  2. Ex3_28 在2SAT问题中,给定一个字句的集合..._第十二次作业

    参考答案 ----------------------------------------------------------------------------------------------- ...

  3. tcp协议简单了解

    HTTP简介 HTTP协议是Hyper Text Transfer Protocol(超文本传输协议)的缩写,是用于从万维网(WWW:World Wide Web )服务器传输超文本到本地浏览器的传送 ...

  4. chkconfig: command not found

    问题描述 Ubuntu 16.04 下安装 Nginx 服务器,在添加 nginx 服务时出现如下信息 # chkconfig --add nginx chkconfig: command not f ...

  5. Java编程之前的复习和练习

    日期:2018.7.14 星期六 博客期:001 今天先是试着写一下博客,最近去青海旅游了,学习时间有点少,但空余时间还是有学习的,不管怎么样吧!先说一下我的这几天的成果——“Bignum”类,虽然很 ...

  6. SpringMVC文件下载与JSON格式

    点击查看上一章 现在JSON这种数据格式是被使用的非常的广泛的,SpringMVC作为目前最受欢迎的框架,它对JSON这种数据格式提供了非常友好的支持,可以说是简单到爆. 在我们SpringMVC中只 ...

  7. vue 之webpack打包工具的使用

    一.什么是webpack? webpack是一个模块打包工具.用vue项目来举例:浏览器它是只认识js,不认识vue的.而我们写的代码后缀大多是.vue的,在每个.vue文件中都可能html.js.c ...

  8. 浅谈FastJson的TypeReference用法

    简单描述:看同事提交的代码,发现有一行代码 似曾相识,但却朦朦胧胧,ε=(´ο`*)))唉很明显自己没掌握呗,于是乎,就百度了一下 干货:对进行泛型的反序列化,使用TypeReference可以明确的 ...

  9. eclipse java ee ide 配置

    配置WEB服务器 windows->preferences->server->Runtime-Environments->Add->选择对应版本->选择安装路径-& ...

  10. easyui 布局之window和panel一起使用时,拉动window宽高时panel不跟随一起变化

    项目开发中布局是每一个组件都由最外层的window和内部的至少一个panel组成,其他的细小组件再依次放到panel中. 问题:当拉动外部的window时我们希望内部的panel的宽高也跟着变化,但是 ...