UOJ #119. 【UR #8】决战圆锥曲线 线段树 暴力 剪枝 卡常数

最新推荐文章于 2021-09-06 23:20:23 发布
最新推荐文章于 2021-09-06 23:20:23 发布
阅读量610 收藏
点赞数 1
分类专栏： 线段树/树状数组 —————————中级数据结构 暴力
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/BlackJack_/article/details/78271775
版权



线段树暴力搞 剪剪枝就好了


#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;

typedef double db;
typedef long long ll;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(ll x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

//(100000005x+20150609) mod 998244353			 ˆy mod 100001
const int N=100100,A=100000005,B=20150609,C=998244353,D=100001;

int n,x0;

inline int rnd()
{x0=(1ll*A*x0+B)%C;return x0/100;}

struct EDGE{int mx,mn;bool rev;}tr[N<<2];

inline void pushup(int k)
{
	tr[k].mx=max(tr[k<<1].mx,tr[k<<1|1].mx);
	tr[k].mn=min(tr[k<<1].mn,tr[k<<1|1].mn);
}

inline void getrever(int k)
{
	int tmp=tr[k].mx;
	tr[k].mx=D-tr[k].mn-1;
	tr[k].mn=D-tmp-1;
}

inline void pushdown(int k)
{
	if(tr[k].rev)
	{
		getrever(k<<1);getrever(k<<1|1);
		tr[k].rev=0;
		tr[k<<1].rev^=1;tr[k<<1|1].rev^=1;
	}
}

void build(int k,int l,int r)
{
	if(l==r)
	{
		tr[k].mx=tr[k].mn=rnd()%D;
		return ;
	}
	int mid=(l+r)>>1;
	build(k<<1,l,mid);build(k<<1|1,mid+1,r);
	pushup(k);
}

void modify(int k,int l,int r,int x,int val)
{
	if(l^r)pushdown(k);
	if(l==r)
	{
		tr[k].mx=tr[k].mn=val;
		return ;
	}
	int mid=(l+r)>>1;
	x<=mid ? modify(k<<1,l,mid,x,val) : modify(k<<1|1,mid+1,r,x,val);
	pushup(k);
}

void rever(int k,int l,int r,int x,int y)
{
	if(l^r)pushdown(k);
	if(l>=x&&r<=y)
	{
		getrever(k);
		tr[k].rev^=1;
		return ;
	}
	int mid=(l+r)>>1;
	if(x<=mid)	rever(k<<1,l,mid,x,y);
	if(y>mid)	rever(k<<1|1,mid+1,r,x,y);
	pushup(k);
}

ll ans;

inline ll col(int x,int y,int a,int b,int c)
{return 1ll*x*a+1ll*y*b+1ll*x*y*c;}

void query(int k,int l,int r,int x,int y,int a,int b,int c)
{
	if(l^r)pushdown(k);
	int mid=(l+r)>>1;
	if(l>=x&&r<=y)
	{
		if(l==r)
		{
			ans=max(ans,col(l,tr[k].mx,a,b,c));
			return ;
		}
		if(col(r,tr[k<<1|1].mx,a,b,c)>ans)	query(k<<1|1,mid+1,r,x,y,a,b,c);
		if(col(mid,tr[k<<1].mx,a,b,c)>ans)	query(k<<1,l,mid,x,y,a,b,c);
		return ;
	}
	if(x<=mid)	query(k<<1,l,mid,x,y,a,b,c);
	if(y>mid)	query(k<<1|1,mid+1,r,x,y,a,b,c);
}

int main()
{
	n=read();int m=read();x0=read();
	build(1,1,n);
	char opt[2];
	register int x,y,a,b,c;
	while(m--)
	{
		scanf("%s",opt);
		switch(opt[0])
		{
			case 'C':
				x=rnd()%n+1;y=rnd()%D;
				modify(1,1,n,x,y);
				break;
			case 'R':
				x=rnd()%n+1;y=rnd()%n+1;if(x>y)swap(x,y);
				rever(1,1,n,x,y);
				break;
			case 'Q':
				x=rnd()%n+1;y=rnd()%n+1;if(x>y)swap(x,y);
				a=read();b=read();c=read();
				ans=0;
				query(1,1,n,x,y,a,b,c);
				print(ans);puts("");
				break;
		}
	}
}
/*
3 3 2705443
C
R
Q 872784 195599 7

13035048532
*/



疯狂卡常版

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define fastcall __attribute__((optimize("-O3")))
#define IL __inline__ __attribute__((always_inline))

typedef double db;
typedef long long ll;

#include<bits/stdc++.h>
using namespace std;
#ifdef ONLINE_JUDGE
char *TT,*mo,but[(1<<15)+2];
#define getchar() ((TT==mo&&(mo=(TT=but)+fread(but,1,1<<15,stdin)),TT==mo)?0:*TT++)
#endif
inline int read()
{
    int x=0,c=0;
    while(c<'0'||c>'9')c=getchar();
	while(c>='0'&&c<='9')x=x*10+c-'0',c=getchar();
    return x;
}
void print(ll x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

//(100000005x+20150609) mod 998244353			 ˆy mod 100001
const int N=100100,A=100000005,B=20150609,C=998244353,D=100001;

int n,x0;

fastcall inline int rnd()
{x0=(1ll*A*x0+B)%C;return x0/100;}

struct EDGE{int mx,mn;bool rev;}tr[N<<2];

fastcall inline void pushup(int k)
{
	tr[k].mx=max(tr[k<<1].mx,tr[k<<1|1].mx);
	tr[k].mn=min(tr[k<<1].mn,tr[k<<1|1].mn);
}

fastcall inline void getrever(int k)
{
	int tmp=tr[k].mx;
	tr[k].mx=D-tr[k].mn-1;
	tr[k].mn=D-tmp-1;
}

fastcall inline void pushdown(int k)
{
	if(tr[k].rev)
	{
		getrever(k<<1);getrever(k<<1|1);
		tr[k].rev=0;
		tr[k<<1].rev^=1;tr[k<<1|1].rev^=1;
	}
}

fastcall void build(int k,int l,int r)
{
	if(l==r)
	{
		tr[k].mx=tr[k].mn=rnd()%D;
		return ;
	}
	int mid=(l+r)>>1;
	build(k<<1,l,mid);build(k<<1|1,mid+1,r);
	pushup(k);
}

fastcall void modify(int k,int l,int r,int x,int val)
{
	if(l^r)pushdown(k);
	if(l==r)
	{
		tr[k].mx=tr[k].mn=val;
		return ;
	}
	int mid=(l+r)>>1;
	x<=mid ? modify(k<<1,l,mid,x,val) : modify(k<<1|1,mid+1,r,x,val);
	pushup(k);
}

fastcall void rever(int k,int l,int r,int x,int y)
{
	if(l^r)pushdown(k);
	if(l>=x&&r<=y)
	{
		getrever(k);
		tr[k].rev^=1;
		return ;
	}
	int mid=(l+r)>>1;
	if(x<=mid)	rever(k<<1,l,mid,x,y);
	if(y>mid)	rever(k<<1|1,mid+1,r,x,y);
	pushup(k);
}

ll ans;

fastcall inline ll col(int x,int y,int a,int b,int c)
{return 1ll*x*a+1ll*y*b+1ll*x*y*c;}

fastcall void query(int k,int l,int r,int x,int y,int a,int b,int c)
{
	if(l^r)pushdown(k);
	int mid=(l+r)>>1;
	if(l>=x&&r<=y)
	{
		if(l==r)
		{
			ans=max(ans,col(l,tr[k].mx,a,b,c));
			return ;
		}
		if(col(r,tr[k<<1|1].mx,a,b,c)>ans)	query(k<<1|1,mid+1,r,x,y,a,b,c);
		if(col(mid,tr[k<<1].mx,a,b,c)>ans)	query(k<<1,l,mid,x,y,a,b,c);
		return ;
	}
	if(x<=mid)	query(k<<1,l,mid,x,y,a,b,c);
	if(y>mid)	query(k<<1|1,mid+1,r,x,y,a,b,c);
}

int main()
{
	n=read();int m=read();x0=read();
	build(1,1,n);
	char opt;
	register int x,y,a,b,c;
	while(m--)
	{
		opt=getchar();
		while(opt!='C'&&opt!='R'&&opt!='Q')opt=getchar();
		switch(opt)
		{
			case 'C':
				x=rnd()%n+1;y=rnd()%D;
				modify(1,1,n,x,y);
				break;
			case 'R':
				x=rnd()%n+1;y=rnd()%n+1;if(x>y)swap(x,y);
				rever(1,1,n,x,y);
				break;
			case 'Q':
				x=rnd()%n+1;y=rnd()%n+1;if(x>y)swap(x,y);
				a=read();b=read();c=read();
				ans=0;
				query(1,1,n,x,y,a,b,c);
				print(ans);puts("");
				break;
		}
	}
}
/*
3 3 2705443
C
R
Q 872784 195599 7

13035048532
*/


BlackJack_
关注
专栏目录
unh6uoj.github.io
03-07
欢迎来到GitHub Pages 您可以使用的来维护和预览Markdown文件中网站的内容。 每当您提交到该存储库时，GitHub Pages都将运行从Markdown文件中的内容重建站点中的页面。 降价促销 Markdown是一种轻巧且易于使用的...
【UR#8】决战圆锥曲线（线段树）
zxyoi_dreamer的博客（不定期诈尸）
04-16 320
传送门 题解： 首先考虑暴力扫一遍出答案。 然后注意到可能成为答案的 i,ji,ji,j 必然满足 i<j,yi>yji < j ,y_i>y_ji<j,yi​>yj​ 由于 yyy 是随机生成，期望的后缀最大值只有 O(log⁡n)O(\log n)O(logn) 个。 于是我们直接用线段树找到后缀最大值即可，维护可以直接维护区间最大值，然后从右儿子开始df...
uoj#119. 【UR #8】决战圆锥曲线（线段树+复杂度分析）
weixin_34061042的博客
02-28 114
题解 传送门 题解 然而要我来说我感觉只是个爆搜啊…… //minamoto #include<bits/stdc++.h> #define R register #define ll long long #define ls (p<<1) #define rs (p<<1|1) #define fp(i,a,b) for(R int i=a,I=b+1;i&...
UOJ #119. 【UR #8】决战圆锥曲线
weixin_30425949的博客
07-15 93
Description Solution 此题解题关键在于数据随机,根据这个进行复杂度分析 不想写题解了 #include<bits/stdc++.h> #define ls (o<<1) #define rs (o<<1|1) using namespace std; typedef long long ll; const int N=2e5+10; in...
[UOJ 118][UOJ Round #8]赴京赶考
某沙茶的OI代码库
06-10 760
题目链接http://uoj.ac/problem/118思路比赛的时候没仔细想就打了个30分的暴力QAQ 比赛完了以后这题一堆满分，我才发现这是个水题对于(ai,bi)(a_i,b_i)，若ai+1≠aia_{i+1}\neq a_i，则无论bib_i取何值，显然(ai,bi)(a_i,b_i)到(ai+1,bi)(a_{i+1},b_i)都是需要花费1单位时间的。因此我们可以在x维度和y维度分
【2020省选模拟】题解
Stargazer的博客
04-17 146
T1： 原题uoj119uoj119uoj119 考虑对于一个区间yyy最大的值(i,yi)(i,y_i)(i,yi​) 显然如果某一个位置比他更优一定是j>ij>ij>i的(j,yj)(j,y_j)(j,yj​) 每次找区间最大值即可 由于数据随机期望只找logloglog次即可 #include<bits/stdc++.h> using namespace std...
uoj-data-converter.zip
09-21
"UOJ Data Converter"是一个专门针对UOJ（Unreal Online Judge）平台的数据转换工具。这个工具的主要目的是帮助用户将特定格式的数据转换为UOJ系统所接受的输入输出格式，以便于在该在线判题系统上进行编程题目的...
UOJ228基础数据结构练习题.cpp
01-11
UOJ228基础数据结构练习题.cpp
南京大学服务器上搭建UOJ系统的保姆级步骤
09-30
从申请南京大学服务器开始，然后搭建docker系统，之后创建UOJ系统用于发布课堂作业，用于帮助助教进行UOJ系统搭建，帮助助教更好的完成老师布置的任务
基于UOJ开发的在线评测系统
最新发布
05-19
【作品名称】：基于UOJ开发的在线评测系统 【适用人群】：适用于希望学习不同技术领域的小白或进阶学习者。可作为毕设项目、课程设计、大作业、工程实训或初期项目立项。 【项目介绍】：更人性化的后台界面，更...
HDU7059 +HDU7116 (线段树暴力修改)
矮纸斜行闲作草
09-06 208
HDU7059 传送门 给了两种修改： 操作1就是区间每个数减去自己的lowbit,由于a[i]是不超过int型范围的，所以操作1最多进行32次，可以单点暴力修改，同时记录这个区间是否都是0，用于剪枝 操作2是区间每个数加上highbit，其实也就是最高位的贡献乘2，我们考虑把最高位与其他位分开维护区间和。 这样操作2就乘了区间乘某个数了，经典的线段树维护操作。 具体要存的信息就是： 1、二进制下最高位之和 ，涉及区间乘2 2、其余位之和，涉及单点暴力修改 3、区间是否都是0，用于剪枝，如果都
UOJ#169 元旦老人与数列(线段树)
hnust_Derker的博客
06-29 437
题目链接 思路： &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;\ \ \ \ 16年的论文题，这里没有办法再用标记来做了， 但是除了历史最值之外都可以利用论文前面提到的方法来做，所以对于历史最值用那个方法额外再维护一下历史最值即可。push_downpush_downpush\_down的时候，当前的记录的历史最值只能传给有最小值的孩子，没有...
UOJ #118. 【UR #8】赴京赶考
sxb_201的博客
07-08 895
先是找规律 然后 贪心 #include #include #include #include using namespace std; const int MAXN = 1e5+10; int a[MAXN*2],b[MAXN*2]; int sum_a[MAXN*2],sum_b[MAXN*2]; int n,m,q; int f(int x1,int y1,int x2,int
线段树暴力修改HDU5930
qq_36876305的博客
01-26 277
先占坑，这题实在不会。。。。 #include #include #include using namespace __gnu_pbds; using namespace std; #define SZ(X) ((int)X.size()) #define mp make_pair #define pb push_back #define RALL(X) X.rbe
BZOJ 2733: [HNOI2012]永无乡
BlackJack
03-15 5019
[HNOI2012]永无乡 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 3211  Solved: 1699 [Submit][Status][Discuss] Description 永无乡包含 n 座岛，编号从 1 到 n，每座岛都有自己的独一无二的重要度，按照重要度可 以将这 n 座岛排名，名次用 1 到 n 来表示。某些岛
BZOJ 3065: 带插入区间K小值 替罪羊树套权值线段树 详解
BlackJack
09-16 956
3065: 带插入区间K小值 Time Limit: 60 Sec  Memory Limit: 512 MB Submit: 3617  Solved: 1173 [Submit][Status][Discuss] Description 从前有n只跳蚤排成一行做早操，每只跳蚤都有自己的一个弹跳力a[i]。跳蚤国王看着这些跳蚤国欣欣向荣的情景，感到非常高兴。这时跳蚤国王决定理性愉悦
AtCoder Regular Contest 063 F : Snuke’s Coloring 2 利用性质 单调栈 线段树
BlackJack
05-25 920
F - すぬけ君の塗り絵 2 / Snuke's Coloring 2Time limit : 4sec / Memory limit : 256MBScore : 1600 pointsProblem StatementThere is a rectangle in the xy-plane, with its lower left corner at (0,0) and its upper r...
SPOJ 3734: Periodni 笛卡尔树 树形dp 组合数学
BlackJack
02-07 814
PERIODNI - Periodni Luka is bored in chemistry class so he is staring at a large periodic table of chemical elements hanging from a wall above the blackboard. To kill time, Luka decided to make his
BZOJ 5039: [Jsoi2014]序列维护 线段树
BlackJack
09-05 749
5039: [Jsoi2014]序列维护 Time Limit: 20 Sec  Memory Limit: 256 MB Submit: 56  Solved: 39 [Submit][Status][Discuss] Description JYY 有一个维护数列的任务。 他希望你能够来帮助他完成。 JYY 现在有一个长度为 N 的序列 a1,a2,…,aN，有如下三种操作：
UOJ#435题目“Simple Tree”
05-30
UOJ#435题目“Simple Tree”要求实现一种数据结构，支持以下两种操作： 1. 修改某个节点的权值； 2. 查询某个节点到根节点路径上权值的最大值。 因为需要查询某个节点到根节点的路径，所以可以考虑使用树链剖分...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值