【BZOJ1227】[SDOI2009]虔诚的墓主人(线段树)

题面

BZOJ

洛谷

题解

显然发现答案就是对于每一个空位置,考虑上下左右各有多少棵树,然后就是这四个方向上树的数量中选\(K\)棵出来的方案数的乘积。显然离散化之后对于答案没有任何影响,所以直接离散化。

然而这样的点数还是\(O(n^2)\)级别,我们把行列拆开考虑。如果我们钦定一行,从左往右看,对于一段连续的空地而言,左右的组合数的乘积是不会变化的,只有上下的乘积会改变,所以可以考虑用一个什么东西维护上下乘积,而左右乘积改变的次数之和恰好等于树的个数,这个是可以接受的。

然而在换行的时候上下乘积是会改变的,然而发现这个的改变次数也恰好是树的个数次,所以总的改变次数就是\(O(n)\)级别的。发现钦定行之后需要维护上下组合数乘积的结果的区间和,用线段树维护即可。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;
#define ll long long
#define MAX 100100
#define lson (now<<1)
#define rson (now<<1|1)
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
int n,m,W,K,X[MAX],Y[MAX],ans;
int Sx[MAX],Sy[MAX],tx,ty;
vector<int> L[MAX];
int R[MAX],sR[MAX];
int C[MAX][22];
bool cmp(int a,int b){return Y[a]<Y[b];}
int t[MAX<<2];
void Modify(int now,int l,int r,int p,int w)
{
if(l==r){t[now]=w;return;}
int mid=(l+r)>>1;
if(p<=mid)Modify(lson,l,mid,p,w);
else Modify(rson,mid+1,r,p,w);
t[now]=t[lson]+t[rson];
}
int Query(int now,int l,int r,int L,int R)
{
if(L>R)return 0;if(L<=l&&r<=R)return t[now];
int mid=(l+r)>>1,ret=0;
if(L<=mid)ret+=Query(lson,l,mid,L,R);
if(R>mid)ret+=Query(rson,mid+1,r,L,R);
return ret;
}
int main()
{
n=read();m=read();W=read();
for(int i=1;i<=W;++i)Sx[++tx]=X[i]=read(),Sy[++ty]=Y[i]=read();
K=read();
sort(&Sx[1],&Sx[tx+1]);sort(&Sy[1],&Sy[ty+1]);
tx=unique(&Sx[1],&Sx[tx+1])-Sx-1;ty=unique(&Sy[1],&Sy[ty+1])-Sy-1;
for(int i=1;i<=W;++i)X[i]=lower_bound(&Sx[1],&Sx[tx+1],X[i])-Sx;
for(int i=1;i<=W;++i)Y[i]=lower_bound(&Sy[1],&Sy[ty+1],Y[i])-Sy;
for(int i=1;i<=W;++i)L[X[i]].push_back(i),R[Y[i]]+=1;
for(int i=1;i<=tx;++i)sort(L[i].begin(),L[i].end(),cmp);
for(int i=1;i<=tx;++i)L[i].push_back(W+1);Y[W+1]=ty+1;
for(int i=0;i<=W;++i)C[i][0]=1;
for(int i=1;i<=W;++i)
for(int j=1;j<=i&&j<=K;++j)
C[i][j]=C[i-1][j]+C[i-1][j-1];
for(int i=1;i<=tx;++i)
{
int ss=0,sum=L[i].size(),l,r;
for(int j=0;j<sum;++j)
{
l=j?Y[L[i][j-1]]:0;r=Y[L[i][j]];
ans+=ss*Query(1,1,ty,l+1,r-1);
if(j==sum-1)break;
ss=C[j+1][K]*C[sum-j-2][K];++sR[r];
Modify(1,1,ty,r,C[sR[r]][K]*C[R[r]-sR[r]][K]);
}
}
if(ans<0)ans+=2147483648ll;
printf("%d\n",ans);
return 0;
}

【BZOJ1227】[SDOI2009]虔诚的墓主人(线段树)的更多相关文章

  1. BZOJ1227 SDOI2009 虔诚的墓主人【树状数组+组合数】【好题】*

    BZOJ1227 SDOI2009 虔诚的墓主人 Description 小W 是一片新造公墓的管理人.公墓可以看成一块N×M 的矩形,矩形的每个格点,要么种着一棵常青树,要么是一块还没有归属的墓地. ...

  2. [BZOJ1227][SDOI2009]虔诚的墓主人 组合数+树状数组

    1227: [SDOI2009]虔诚的墓主人 Time Limit: 5 Sec  Memory Limit: 259 MBSubmit: 1433  Solved: 672[Submit][Stat ...

  3. BZOJ1227 [SDOI2009]虔诚的墓主人 【树状数组】

    题目 小W 是一片新造公墓的管理人.公墓可以看成一块N×M 的矩形,矩形的每个格点,要么种着一棵常青树,要么是一块还没有归属的墓地.当地的居民都是非常虔诚的基督徒,他们愿意提前为自己找一块合适墓地.为 ...

  4. bzoj1227: [SDOI2009]虔诚的墓主人(树状数组,组合数)

    传送门 首先,对于每一块墓地,如果上下左右各有$a,b,c,d$棵树,那么总的虔诚度就是$C_k^a*C_k^b*C_k^c*C_k^d$ 那么我们先把所有的点都给离散,然后按$x$为第一关键字,$y ...

  5. bzoj1227 [SDOI2009]虔诚的墓主人(组合公式+离散化+线段树)

    1227: [SDOI2009]虔诚的墓主人 Time Limit: 5 Sec  Memory Limit: 259 MBSubmit: 803  Solved: 372[Submit][Statu ...

  6. [bzoj1227] [SDOI2009]虔诚的墓主人

    终于填上了这个万年巨坑....从初二的时候就听说过这题...然后一直不敢写QAQ 现在感觉也不是很烦(然而我还是写麻烦了 离散化一波,预处理出组合数什么的.. 要维护对于当前行,每列上方和下方节点凑出 ...

  7. luogu2154 [SDOI2009] 虔诚的墓主人 离散化 树状数组 扫描线

    题目大意 公墓可以看成一块N×M的矩形,矩形的每个格点,要么种着一棵常青树,要么是一块还没有归属的墓地.一块墓地的虔诚度是指以这块墓地为中心的十字架的数目,一个十字架可以看成中间是墓地,墓地的正上.正 ...

  8. bzoj1227 P2154 [SDOI2009]虔诚的墓主人

    P2154 [SDOI2009]虔诚的墓主人 组合数学+离散化+树状数组 先看题,结合样例分析,易得每个墓地的虔诚度=C(正左几棵,k)*C(正右几棵,k)*C(正上几棵,k)*C(正下几棵,k),如 ...

  9. Bzoj 1227: [SDOI2009]虔诚的墓主人 树状数组,离散化,组合数学

    1227: [SDOI2009]虔诚的墓主人 Time Limit: 5 Sec  Memory Limit: 259 MBSubmit: 895  Solved: 422[Submit][Statu ...

随机推荐

  1. 20155207 EXP7 EXP8 EXP9 实验补交

    20155207 EXP7 EXP8 EXP9 实验补交 20155207 EXP7 网络欺诈技术防范 20155207 EXP8 Web基础 20155207 <网络对抗> Exp9 W ...

  2. 网络对抗技术 2017-2018-2 20152515 Exp3 免杀原理与实践

    基础问题回答 (1)杀软是如何检测出恶意代码的? 答:分析恶意程序的行为特征,分析其代码流将其性质归类于恶意代码. (2)免杀是做什么? 答:一般是对恶意软件做处理,让它不被杀毒软件所检测,也是渗透测 ...

  3. 20155238 Java第13周课堂实践

    类定义 实验内容及要求 设计并实现一个Book类,定义义成Book.java,Book 包含书名,作者,出版社和出版日期,这些数据都要定义getter和setter.定义至少三个构造方法,接收并初始化 ...

  4. 2017-2018-2 20155333 《网络对抗技术》 Exp1 PC平台逆向破解

    2017-2018-2 20155333 <网络对抗技术> Exp1 PC平台逆向破解 1. 逆向及Bof基础实践说明 1.1 实践目标 本次实践的对象是一个名为pwn1的linux可执行 ...

  5. mfc CListBox

    通过ID操作对象 CListBox(列表框)控件 CListBox类常用成员 CListBox插入数据 CListBox删除数据 CListBox运用示例 一.CListBox类常用成员 CListB ...

  6. Promise 执行顺序

    加深印象 setTimeout(function() { console.log("timeout-start"); }, 200) // 改为100呢 console.log(& ...

  7. python高并发和多线程的关系

    “高并发和多线程”总是被一起提起,给人感觉两者好像相等,实则 高并发 ≠ 多线程 多线程是完成任务的一种方法,高并发是系统运行的一种状态,通过多线程有助于系统承受高并发状态的实现.   高并发是一种系 ...

  8. laraver框架学习

    最近开始学习laravel框架,这个框架在国外很流行,近些年开始在国内流行.自己而是刚开始学习这个框架. 使用composer 更新系统内的依赖包 在终端输入:composer update Entr ...

  9. C#_正则表达式

    概述 正则表达式,主要是用符号描述了一类特定的文本(模式).而正则表达式引擎则负责在给定的字符串中,查找到这一特定的文本. 本文主要是列出常用的正则表达式符号,加以归类说明.本文仅仅是快速理解了正则表 ...

  10. 软件测试_APP测试_兼容性测试

    APP的兼容测试主要就是测试APP的安装.启动.运行.卸载测试,以及安装时间 .启动时间.CPU占用.内存占用.流量耗用.电量耗用等性能上的测试. 兼容性测试点: 一.内部兼容性: 1.与本地和其他主 ...