bzoj1227 P2154 [SDOI2009]虔诚的墓主人
组合数学+离散化+树状数组
先看题,结合样例分析,易得每个墓地的虔诚度=C(正左几棵,k)*C(正右几棵,k)*C(正上几棵,k)*C(正下几棵,k),如果任意一遍的棵树<k,则虔诚度=0。
所以我们可以预处理出C(w,k)。
再看数据范围:“对于100%的数据,满足1 ≤ N, M ≤ 1,000,000,000,0 ≤ xi ≤ N,0 ≤ yi ≤ M,1 ≤ W ≤ 100,000,1 ≤ k ≤ 10。”
对比一下 n,m 和 w 的大小,肯定要离散化。
排序离散化后,我们可以枚举。
按 x 坐标排序后,我们在枚举时可以方便的得出C(正上几棵,k)*C(正下几棵,k)
但是 C(正左几棵,k)*C(正右几棵,k) 和 y 坐标有关系,是动态变化的。
所以我们可以用树状数组维护正左几棵,正右几棵。
然后注意一些细节:2147483648 在我的编译器有警告qaq,所以只能直接用 2147483648ll 了qwq
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
using namespace std;
inline int Int(){
char c=getchar(); int x=;
while(!isdigit(c)) c=getchar();
while(isdigit(c)) x=(x<<)+(x<<)+(c^),c=getchar();
return x;
}
inline int min(int &a,int &b) {return a<b ?a:b;}
const int maxw=1e5+;
struct data{int x,y;} a[maxw];
inline bool cmp1(const data &x0,const data &x1){return x0.x<x1.x||(x0.x==x1.x&&x0.y<x1.y);}
inline bool cmp2(const data &x0,const data &x1){return x0.y<x1.y||(x0.y==x1.y&&x0.x<x1.x);}
int n,m,w,k,ed,ans,c[maxw][],tmp[maxw],totx[maxw],toty[maxw],p[maxw],sum[maxw],last[maxw];
inline void add(int x,int y){for(;x<=ed;x+=(x&-x)) sum[x]+=y;}
inline int query(int x) {int ans=; for(;x;x-=(x&-x)) ans+=sum[x]; return ans;}
int main(){
n=Int(); m=Int(); w=Int();
for(int i=;i<=w;++i) a[i].x=Int(),a[i].y=Int();
k=Int(); int cnt; for(int i=;i<=w;++i) c[i][]=;
for(int i=;i<=w;++i)
for(int j=;j<=min(i,k);++j)
c[i][j]=c[i-][j]+c[i-][j-];
//预处理组合数
sort(a+,a+w+,cmp1); cnt=; tmp[]=;
for(int i=;i<=w;++i) cnt= a[i].x==a[i-].x ? cnt:cnt+,tmp[i]=cnt;
for(int i=;i<=w;++i) ++totx[a[i].x=tmp[i]];
sort(a+,a+w+,cmp2); cnt=; tmp[]=;
for(int i=;i<=w;++i) cnt= a[i].y==a[i-].y ? cnt:cnt+,tmp[i]=cnt;
for(int i=;i<=w;++i) ++toty[a[i].y=tmp[i]];
ed=a[w].y;
//横纵坐标离散化
sort(a+,a+w+,cmp1);
for(int i=;i<=w;++i){
cnt= i==||a[i].x!=a[i-].x ? :cnt+;
int u=a[i].y,v= (++p[u])>=k&&toty[u]-p[u]>=k ? 1LL*c[p[u]][k]*c[toty[u]-p[u]][k]%2147483648ll:;
add(u,v-last[u]); last[u]=v;
if(i==w||a[i].x!=a[i+].x||a[i+].y-a[i].y<=||cnt<k||totx[a[i].x]-cnt<k) continue;
ans+= 1LL*c[cnt][k]*c[totx[a[i].x]-cnt][k]%2147483648ll*(query(a[i+].y-)-query(a[i].y))%2147483648ll;
}//树状数组维护
printf("%d",ans< ?(ans+2147483648ll)%2147483648ll:ans);
return ;
}
bzoj1227 P2154 [SDOI2009]虔诚的墓主人的更多相关文章
- 【BZOJ1227】[SDOI2009]虔诚的墓主人(线段树)
[BZOJ1227][SDOI2009]虔诚的墓主人(线段树) 题面 BZOJ 洛谷 题解 显然发现答案就是对于每一个空位置,考虑上下左右各有多少棵树,然后就是这四个方向上树的数量中选\(K\)棵出来 ...
- 【BZOJ1227】[SDOI2009]虔诚的墓主人
E. 虔诚的墓主人 题目描述 小W 是一片新造公墓的管理人.公墓可以看成一块N×M 的矩形,矩形的每个格点,要么种着一棵常青树,要么是一块还没有归属的墓地.当地的居民都是非常虔诚的基督徒,他们愿意提前 ...
- 【bzoj1227】 SDOI2009—虔诚的墓主人
http://www.lydsy.com/JudgeOnline/problem.php?id=1227 (题目链接) 题意 一个n*m的公墓,一个点上要么是墓地,要么是常青树,给出一个数K,并规定每 ...
- P2154 [SDOI2009]虔诚的墓主人 树状数组
https://www.luogu.org/problemnew/show/P2154 题意 在一个坐标系中,有w(1e5)个点,这个图中空点的权值是正上,正下,正左,正右各取k个的排列组合情况.计算 ...
- [洛谷P2154] SDOI2009 虔诚的墓主人
问题描述 小W是一片新造公墓的管理人.公墓可以看成一块N×M的矩形,矩形的每个格点,要么种着一棵常青树,要么是一块还没有归属的墓地. 当地的居民都是非常虔诚的基督徒,他们愿意提前为自己找一块合适墓地. ...
- P2154 [SDOI2009]虔诚的墓主人
略有一点点思维的题. 首先,如果一个点上,下,左,右分别有\(a,b,c,d\)棵树,那这个点的十字架方案为\(C_{a}^{k}C_{b}^{k}C_{c}^{k}C_{d}^{k}\). 按x坐标 ...
- luogu P2154 [SDOI2009]虔诚的墓主人
luogu 下面记一个点上下左右点数分别为\(u_i,d_i,l_i,r_i\) 枚举每个中间点太慢了,考虑枚举两个点之间横的一条线段,这里面的点左边点数目都相同,右边点数目都相同,然后只要查一下区间 ...
- BZOJ1227 SDOI2009 虔诚的墓主人【树状数组+组合数】【好题】*
BZOJ1227 SDOI2009 虔诚的墓主人 Description 小W 是一片新造公墓的管理人.公墓可以看成一块N×M 的矩形,矩形的每个格点,要么种着一棵常青树,要么是一块还没有归属的墓地. ...
- bzoj1227 [SDOI2009]虔诚的墓主人(组合公式+离散化+线段树)
1227: [SDOI2009]虔诚的墓主人 Time Limit: 5 Sec Memory Limit: 259 MBSubmit: 803 Solved: 372[Submit][Statu ...
随机推荐
- linux 执行远程linux上的shell脚本或者命令以及scp 上传文件到ftp--免密码登陆
场景:在linux A 上执行Linux B上的shell脚本和命令 步骤1.设置ssh免登陆 1.SSH无密码登录 # 本地服务器执行(A机器):生成密钥对 ssh-keygen -t dsa -P ...
- Css3动画属性总汇
http://css3lib.alloyteam.com/uilib/animation/demo1/#cta Css3动画属性总汇 Stay hungry. Stay foolish. Attent ...
- Hybrid设计--账号体系的建设
前后端分离:开发效率高,没有SEO 现在是重客户端设计:交互和业务逻辑是前端来写,适合做前后端分离.对前端更友好,提高了效率. 传统模式开发:整个业务逻辑是server端写,不适合做前后端分离.ser ...
- 解决mysql的内存表“table is full”错误
最后参考http://blog.sina.com.cn/s/blog_6942a1590101429h.html 来解决,摘录下核心 后来GOOGLE得知,需要重建该表才可以. 1. 设置新的参数 m ...
- Go断后,Dart冲前,Google的野心
今天,我要讲讲新的前端web语言-Dart 我是网上看到这段话,激起了兴趣,虽然我不能算是真正的web前端从业者!哈哈 [ Dart语言由谷歌制造,用来替代Javascript,弥补Javascrip ...
- django模型(增删改查等)
Django提供的模型操作起来数据库非常方便 以自定义模型名字Demo为例: 获取所有数据:all() demo.objects.all() #返回queryset对象 #这个语句和Select * ...
- 安装FusionInsight
1.在华为平台上下载整体客户端,不建议下载单个组件客户端,后期关联测试还是要装上的. 2.下载后需要将服务器上的客户端拷贝到本地.打开xShell,新建会话,登陆本地虚拟机上的Linux系统(19 ...
- install apache-activemq
进入bin/linux-x86-64下 cd apache-activemq-5.14.1/bin/linux-x86-64/ 启动 ./activemq start 五.启动成功后,访问 activ ...
- poj2114 寻找树上存在长度为k点对,树上的分治
寻找树上存在长度为k点对,树上的分治 代码和 这个 差不多 ,改一下判断的就好 #include <iostream> #include <algorithm> #inc ...
- 20155228 获取技能的成功经验和关于C语言学习的调查
内容提要 你有什么技能比大多人(超过90%以上)更好?针对这个技能的获取你有什么成功的经验?与老师博客中的学习经验有什么共通之处? 有关C语言学习的调查 你是怎么学习C语言的?(作业,实验,教材,其他 ...