BZOJ 2244 [SDOI2011]拦截导弹 (三维偏序CDQ+线段树)
题目大意:
不愧为SDOI的duliu题
第一问?二元组的最长不上升子序列长度?裸的三维偏序问题,直接上$CDQ$
由于是不上升,需要查询某一范围的最大值,并不是前缀最大值,建议用线段树实现
第二问是个什么玩意??
画画图发现需要正反各做一次$CDQ$来统计
如果某个位置正反的答案$-1$就是最长长度
那么它被选择的次数就是 正着统计作为末尾的次数*反着统计作为末尾的次数
概率就是这个值/总次数
又发现某个位置作为末尾的次数可能非常非常大!
比如$1\;1\;2\;2\;3\;3\;4\;4\;5\;5....$这个值甚至达到了$2^{n/2}$
而题目又是让我们求概率,所以这个次数必须要用$double$存
蒟蒻的代码写得比较恶心..
另外听一些神犇说$sort$中的$cmp$里不能写$<=$或$>=$,不然会$RE$
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N1 50100
#define ll long long
#define dd double
#define inf 0x3f3f3f3f3f3f3f3fll
using namespace std; int gint()
{
int ret=,fh=;char c=getchar();
while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
while(c>=''&&c<=''){ret=ret*+c-'';c=getchar();}
return ret*fh;
}
int n,m,mh,mv,K;
inline void chk(int &ma,dd &sum,int w,dd num)
{
if(w>ma) ma=w,sum=num;
else if(w==ma) sum+=num;
} struct SEG{
int ma[N1<<]; dd sum[N1<<];
void pushup(int rt)
{
ma[rt]=; sum[rt]=;
chk(ma[rt],sum[rt],ma[rt<<],sum[rt<<]);
chk(ma[rt],sum[rt],ma[rt<<|],sum[rt<<|]);
}
void update(int x,int l,int r,int rt,int w,dd num)
{
if(l==r) { chk(ma[rt],sum[rt],w,num); return; }
int mid=(l+r)>>;
if(x<=mid) update(x,l,mid,rt<<,w,num);
else update(x,mid+,r,rt<<|,w,num);
pushup(rt);
}
void query(int L,int R,int l,int r,int rt,int &w,dd &num)
{
if(L<=l&&r<=R) { chk(w,num,ma[rt],sum[rt]); return; }
int mid=(l+r)>>;
if(L<=mid) query(L,R,l,mid,rt<<,w,num);
if(R>mid) query(L,R,mid+,r,rt<<|,w,num);
}
void clr(int x,int l,int r,int rt)
{
ma[rt]=; sum[rt]=;
if(l==r) return; int mid=(l+r)>>;
if(x<=mid) clr(x,l,mid,rt<<);
else clr(x,mid+,r,rt<<|);
}
}s;
struct node{int h,v,t,ans; dd sum;}a[N1],c[N1],tmp[N1];
int h[N1],v[N1],que[N1],tl; int cmp1(node s1,node s2){ if(s1.h!=s2.h) return s1.h>s2.h; return s1.v>=s2.v; }
int cmp2(node s1,node s2){ if(s1.h!=s2.h) return s1.h<s2.h; return s1.v<=s2.v; }
int cmpp(node s1,node s2){ if(s1.h!=s2.h) return s1.h>s2.h; return s1.v>s2.v; }
int cmpn(node s1,node s2){ if(s1.h!=s2.h) return s1.h<s2.h; return s1.v<s2.v; } void CDQ1(int L,int R)
{
if(R-L<=) return;
int M=(L+R)>>,i,j,k,cnt;
for(i=L,j=M,k=L;k<R;k++)
{
if(a[k].t<M) tmp[i++]=a[k];
else tmp[j++]=a[k];
}
for(k=L;k<R;k++) a[k]=tmp[k];
CDQ1(L,M);
for(i=L,j=M;i<M&&j<R;)
{
if(cmp1(a[i],a[j])) { s.update(a[i].v,,mv,,a[i].ans+,a[i].sum); que[++tl]=i; i++; }
else { s.query(a[j].v,mv,,mv,,a[j].ans,a[j].sum); j++; }
}
while(j<R) { s.query(a[j].v,mv,,mv,,a[j].ans,a[j].sum); j++;}
while(tl) { s.clr(a[que[tl--]].v,,mv,); }
CDQ1(M,R);
for(i=L,j=M,cnt=;i<M&&j<R;)
{
if(cmp1(a[i],a[j])) { tmp[++cnt]=a[i]; i++; }
else { tmp[++cnt]=a[j]; j++;}
}
while(i<M) tmp[++cnt]=a[i++];
while(j<R) tmp[++cnt]=a[j++];
for(k=L;k<R;k++) a[k]=tmp[k-L+];
}
void CDQ2(int L,int R)
{
if(R-L<=) return;
int M=(L+R)>>,i,j,k,cnt;
for(i=L,j=M,k=L;k<R;k++)
{
if(a[k].t<M) tmp[i++]=a[k];
else tmp[j++]=a[k];
}
for(k=L;k<R;k++) a[k]=tmp[k];
CDQ2(L,M);
for(i=L,j=M;i<M&&j<R;)
{
if(cmp2(a[i],a[j])) { s.update(a[i].v,,mv,,a[i].ans+,a[i].sum); que[++tl]=i; i++; }
else { s.query(,a[j].v,,mv,,a[j].ans,a[j].sum); j++; }
}
while(j<R) { s.query(,a[j].v,,mv,,a[j].ans,a[j].sum); j++;}
while(tl) { s.clr(a[que[tl--]].v,,mv,); }
CDQ2(M,R);
for(i=L,j=M,cnt=;i<M&&j<R;)
{
if(cmp2(a[i],a[j])) { tmp[++cnt]=a[i]; i++; }
else { tmp[++cnt]=a[j]; j++;}
}
while(i<M) tmp[++cnt]=a[i++];
while(j<R) tmp[++cnt]=a[j++];
for(k=L;k<R;k++) a[k]=tmp[k-L+];
}
dd ret[N1];
int p[N1];
int de; int main()
{
scanf("%d",&n);
int i,j,x,y,ans=; dd sum=;
for(i=;i<=n;i++) { h[i]=a[i].h=gint(); v[i]=a[i].v=gint(); a[i].t=i; a[i].ans=; a[i].sum=; }
sort(h+,h+n+); mh=unique(h+,h+n+)-(h+); sort(v+,v+n+); mv=unique(v+,v+n+)-(v+);
for(i=;i<=n;i++) { a[i].h=lower_bound(h+,h+mh+,a[i].h)-h; a[i].v=lower_bound(v+,v+mv+,a[i].v)-v;}
sort(a+,a+n+,cmpp);
CDQ1(,n+);
for(i=;i<=n;i++) { c[i]=a[i]; p[a[i].t]=i; chk(ans,sum,a[i].ans,a[i].sum); }
sort(a+,a+n+,cmpn); for(i=;i<=n;i++) { a[i].ans=; a[i].sum=; a[i].t=n-a[i].t+; }
CDQ2(,n+);
for(i=;i<=n;i++)
{
j=p[n-a[i].t+];
if(a[i].ans+c[j].ans-<ans) ret[c[j].t]=;
else ret[c[j].t]=1.0*c[j].sum*a[i].sum/sum;
}
printf("%d\n",ans);
for(i=;i<=n;i++) printf("%.5lf ",ret[i]);
puts(""); return ;
}
BZOJ 2244 [SDOI2011]拦截导弹 (三维偏序CDQ+线段树)的更多相关文章
- bzoj 2244 [SDOI2011]拦截导弹(DP+CDQ分治+BIT)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2244 [题意] 给定n个二元组,求出最长不上升子序列和各颗导弹被拦截的概率. [思路] ...
- bzoj 2244 [SDOI2011]拦截导弹(dp+CDQ+树状数组)
传送门 题解 看了半天完全没发现这东西和CDQ有什么关系…… 先把原序列翻转,求起来方便 然后把每一个位置表示成$(a,b,c)$其中$a$表示位置,$b$表示高度,$c$表示速度,求有多少个位置$a ...
- BZOJ 2244: [SDOI2011]拦截导弹 DP+CDQ分治
2244: [SDOI2011]拦截导弹 Description 某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度.并且能够拦截 ...
- bzoj 2244: [SDOI2011]拦截导弹 cdq分治
2244: [SDOI2011]拦截导弹 Time Limit: 30 Sec Memory Limit: 512 MBSec Special JudgeSubmit: 237 Solved: ...
- BZOJ 2244: [SDOI2011]拦截导弹 (CDQ分治 三维偏序 DP)
题意 略- 分析 就是求最长不上升子序列,坐标取一下反就是求最长不下降子序列,比较大小是二维(h,v)(h,v)(h,v)的比较.我们不看概率,先看第一问怎么求最长不降子序列.设f[i]f[i]f[i ...
- BZOJ 2244 [SDOI2011]拦截导弹 ——CDQ分治
三维偏序,直接CDQ硬上. 正反两次CDQ统计结尾的方案数,最后统计即可. #include <cstdio> #include <cstring> #include < ...
- BZOJ 2244: [SDOI2011]拦截导弹 [CDQ分治 树状数组]
传送门 题意:三维最长不上升子序列以及每个元素出现在最长不上升子序列的概率 $1A$了好开心 首先需要从左右各求一遍,长度就是$F[0][i]+F[1][i]-1$,次数就是$G[0][i]*G[1] ...
- bzoj 2244: [SDOI2011]拦截导弹
#include<cstdio> #include<iostream> #include<algorithm> #define M 100009 using nam ...
- BZOJ:2244: [SDOI2011]拦截导弹
问题: printf("%.5f ",0):为什么错了? 注意: 初始值很重要 题解: 三维偏序问题: 记录从前往后最长上升子序列长度pref,条数preg 从后往前suff,su ...
随机推荐
- C#常用 API函数大全
常用Windows API1. API之网络函数WNetAddConnection 创建同一个网络资源的永久性连接WNetAddConnection2 创建同一个网络资源的连接WNetAddConne ...
- tomcat 内存溢出问题(OutOfMemoryError: PermGen space)
导入公司项目的时候出现的问题,在此记录处理方法. tomcat在启动的时候报错:OutOfMemoryError: PermGen space PermGen space的全称是Permanent G ...
- Jquery Math ceil()、floor()、round()比较与用法
Math.ceil():向上取值 如:Math.ceil(2.1) -- 结果为 3 Math.ceil(-2.1) -- 结果为-2 结论:正入 负舍 Math.floor(): 先下取值 入 ...
- JAVA基础总结【面试】
前言 近间陆续面试了不少的求职的前(JAVA).后(WEB)端开发人员,包括实习生.应届毕业生.一两年工作经验的.也有三四年工作经验的,也算见过了比较多的开发人员,想在这里做个总结,本次主要讲一讲面试 ...
- POJ 3709
简单的单调队列优化,注意是哪些点加入队列即可. #include <iostream> #include <cstdio> #include <algorithm> ...
- POJ 2407
裸 的求欧拉函数 #include <iostream> #include <cstdio> #include <cstring> #include <alg ...
- HBase使用flush命令之后存储的位置
HBase使用flush命令之后存储的位置 根据系统安装位置的不一样而不一样,当前是在: hadoop fs -ls /apps/hbase/data/data/default/t1 下面: 使用ha ...
- 【iOS开发系列】UIDevice设备信息
[1] 推断是否是横向屏: BOOL b=UIDeviceOrientationIsLandscape([UIDevice currentDevice].orientation); 获取设备uniqu ...
- iOS CoreImage图片处理动态渲染(滤镜)
// // ViewController.m // CoreImageOfDong // // Created by Dong on 15/6/30. // Copyright (c) 201 ...
- NumPy和Pandas常用库
NumPy和Pandas常用库 1.NumPy NumPy是高性能科学计算和数据分析的基础包.部分功能如下: ndarray, 具有矢量算术运算和复杂广播能力的快速且节省空间的多维数组. 用于对整组数 ...