(2016北京集训十二)【xsy1542】疯狂求导
题解:
这题看起来很难。。。但是实际上并没有想象中的那么难
第一眼看上去不会求导公式怎么办?不要紧,题目背景非常良心的给出了题目中的导数计算公式
求完导合并同类项很恶心怎么办?不要紧,样例解释说明了不需要合并同类项(然后有许多人因为这个爆〇了)
一看这种题目形式明显就是大数据结构,外面的序列明显线段树维护,次数也可以用线段树,但是线段树套线段树容易MLE;
所以用树状数组套线段树实现
具体就是以1~n为下标建线段树,外面用树状数组维护次数,每次在树状数组上查询即可
写完过样例直接1A就是爽
代码:
1 #include<algorithm> 2 #include<iostream> 3 #include<cstring> 4 #include<cstdio> 5 #include<cmath> 6 #include<queue> 7 #define inf 2147483647 8 #define eps 1e-9 9 #define lb(x) (x&-x) 10 using namespace std; 11 typedef long long ll; 12 const int N=1000000,bit=524288; 13 struct node{ 14 int num,ls,rs; 15 ll v; 16 }t[10000001]; 17 int n,m,ans=0,op,l,r,tot=0,rts[1000001]; 18 int TheseNumbersDoNotMakeAnySense; 19 ll s; 20 void newn(int &u){ 21 if(!u)u=++tot; 22 } 23 void pd(int u,int l,int r){ 24 if(t[u].num){ 25 int mid=(l+r)/2; 26 newn(t[u].ls); 27 newn(t[u].rs); 28 t[t[u].ls].num+=t[u].num; 29 t[t[u].ls].v+=t[u].num*(mid-l+1); 30 t[t[u].rs].num+=t[u].num; 31 t[t[u].rs].v+=t[u].num*(r-mid); 32 t[u].num=0; 33 } 34 } 35 void updata(int l,int r,int &u,int L,int R,int v){ 36 newn(u); 37 if(L<=l&&r<=R){ 38 t[u].num+=v; 39 t[u].v+=v*(r-l+1); 40 return; 41 } 42 int mid=(l+r)/2; 43 pd(u,l,r); 44 if(L<=mid)updata(l,mid,t[u].ls,L,R,v); 45 if(mid<R)updata(mid+1,r,t[u].rs,L,R,v); 46 t[u].v=t[t[u].ls].v+t[t[u].rs].v; 47 } 48 int query(int l,int r,int u,int L,int R){ 49 if(!u)return 0; 50 if(L<=l&&r<=R){ 51 return t[u].v; 52 } 53 int mid=(l+r)/2,ret=0; 54 pd(u,l,r); 55 if(L<=mid)ret=query(l,mid,t[u].ls,L,R); 56 if(mid<R)ret+=query(mid+1,r,t[u].rs,L,R); 57 return ret; 58 } 59 void ins(int l,int r,int s){ 60 for(;s<=N;s+=lb(s)){ 61 updata(1,n,rts[s],l,r,1); 62 } 63 } 64 ll ask(int l,int r,int s){ 65 ll ret=0; 66 for(;s;s-=lb(s)){ 67 ret+=query(1,n,rts[s],l,r); 68 } 69 return ret; 70 } 71 void work(int l,int r,ll s){ 72 ll ret=ask(l,r,N),tmp=0; 73 if(ret<=s){ 74 printf("1 %lld\n",ret); 75 ans=1; 76 return; 77 } 78 s=ret-s; 79 for(int bt=bit;bt;bt/=2){ 80 if(tmp+bt<=N){ 81 int q=query(1,n,rts[tmp+bt],l,r); 82 //printf("%lld %lld %lld\n",tmp,bt,q); 83 if(s>q)s-=q,tmp+=bt; 84 } 85 } 86 ll rt=ask(l,r,tmp+1); 87 rt=ret-rt; 88 printf("%lld %lld\n",tmp+2,rt); 89 ans=tmp+2; 90 } 91 int main(){ 92 scanf("%d%d",&n,&m); 93 for(int i=1;i<=n;i++){ 94 scanf("%d",&TheseNumbersDoNotMakeAnySense); 95 } 96 for(int i=1;i<=m;i++){ 97 scanf("%d%d%d%lld",&op,&l,&r,&s); 98 if(op==0){ 99 s^=ans; 100 ins(l,r,s); 101 }else{ 102 work(l,r,s); 103 } 104 } 105 return 0; 106 }
转载于:https://www.cnblogs.com/dcdcbigbig/p/9696111.html
(2016北京集训十二)【xsy1542】疯狂求导相关推荐
- (2016北京集训十)【xsy1529】小Q与进位制 - 分治FFT
题意很简单,就是求这个数... 其实场上我想出了分治fft的正解...然而不会打...然后打了个暴力fft挂了... 没啥好讲的,这题很恶心,卡常卡精度还爆int,要各种优化,有些dalao写的很复杂 ...
- 问题六十二:怎么求一元十次方程在区间内的所有不相等的实根(2)——修正“区间端点零值”问题
前续"问题六十二:怎么求一元十次方程在区间内的所有不相等的实根"和"问题五十九:怎么求一元六次方程在区间内的所有不相等的实根(3)--修正一个问题" 不管是求解 ...
- SLAM总结(二)-数学基础之求导和线性方程求解
SLAM总结(二)-数学基础之求导和线性方程求解 1.求导:高数中常见的是一个函数对一个自变量求导,属于标量对标量求导.在SLAM问题中,函数是目标函数(残差项,约束项),一般包含多个函数,用多维列向 ...
- PyTorch基础(二)-----自动求导Autograd
一.前言 上一篇文章中提到PyTorch提供了两个重要的高级功能,分别是: 具有强大的GPU加速的张量计算(如NumPy) 包含自动求导系统的的深度神经网络 第一个特性我们会在之后的学习中用到,这里暂 ...
- 机器学习中的矩阵向量求导(二) 矩阵向量求导之定义法
在机器学习中的矩阵向量求导(一) 求导定义与求导布局中,我们讨论了向量矩阵求导的9种定义与求导布局的概念.今天我们就讨论下其中的标量对向量求导,标量对矩阵求导, 以及向量对向量求导这三种场景的基本求解 ...
- LSGO软件技术团队2015~2016学年第十二周(1116~1122)总结
团队简述: LSGO软件技术团队成立于2010年10月,主要从事的应用方向为互联网与移动互联网(UI设计,前端开发,后台开发),地理信息系统:研究方向为大数据处理与机器学习.成立几年来为学校培养了一批 ...
- 八十二、归并排序求取复杂的逆序数
@Author:Runsen 逆序数,我在很多的面试题都见过,本质上来说难度是比较大,因为如果使用暴力法当数据量一大,必然就会爆掉.你现在就要记住逆序数就是考归并排序. 逆序数 给定一个数组array ...
- 【2016北京集训测试赛】river
HINT 注意是全程不能经过两个相同的景点,并且一天的开始和结束不能用同样的交通方式. [吐槽] 嗯..看到这题的想法的话..先想到了每个点的度为2,然后就有点不知所措了 隐隐约约想到了网络流,但并没 ...
- 2016北京集训测试赛(十三) Problem B: 网络战争
Solution KD tree + 最小割树 转载于:https://www.cnblogs.com/ZeonfaiHo/p/7420354.html
最新文章
- 教你打包Java程序,jar转exe随处可跑
- Fiddler抓包使用教程-过滤
- commons-lang包中我们常用的类的作用
- [svc]cfssl模拟https站点-探究浏览器如何校验证书
- php定义object数据类型,PHP数据类型(4):对象object
- centos中安装nginx
- win10远程登陆deepin xrdp
- grafana安装配置 influxdb安装配置
- js中如何获取对象的长度和名称
- 小米路由2+安装+php,小米路由器作为二级路由,完美使用APP和文件共享
- 微信小程序之海报生成
- c语言星空特效源代码,星空 - 网页特效代码|网页特效观止 - 让你的网页靓起来!...
- 考研刷题小程序的盈利模式分析
- 内德-米德方法——《数值计算方法》
- 容联云 PHP 实现短信发送验证码
- 在内网使用DBL落地网关,连接公网的freeswitch
- 3ds文件转换为.X文件
- zan php framework,zan框架入门(一)——协程
- linux usbnet驱动
- 专递课堂、名师课堂、名校网络课堂的建设方案
热门文章
- 《C#图解教程》读书笔记之六:接口和转换
- angular中父组件给子组件传值-@input
- spring注解方式注入bean
- css cursor 的可选值(鼠标的各种样式)
- Linux环境下安装Tigase XMPP Server
- Ubuntu 下安装adobe flash player
- sql中exists替换in的区别
- java版spring cloud+spring boot+redis多租户社交电子商务平台 (十一)docker部署spring cloud项目...
- MyBatis 几个常用的 jdbcType数据
- 无微不至:调整_lm_cache_res_cleanup解决Shared Pool 的4031问题