1227: [SDOI2009]虔诚的墓主人

Time Limit: 5 Sec  Memory Limit: 259 MB
Submit: 1083  Solved: 514
[Submit][Status][Discuss]

Description

小W 是一片新造公墓的管理人。公墓可以看成一块N×M 的矩形,矩形的每个格点,要么种着一棵常青树,要么是一块还没有归属的墓地。当地的居民都是非常虔诚的基督徒,他们愿意提前为自己找一块合适墓地。为了体现自己对主的真诚,他们希望自己的墓地拥有着较高的虔诚度。一块墓地的虔诚度是指以这块墓地为中心的十字架的数目。一个十字架可以看成中间是墓地,墓地的正上、正下、正左、正右都有恰好k 棵常青树。小W 希望知道他所管理的这片公墓中所有墓地的虔诚度总和是多少

Input

第一行包含两个用空格分隔的正整数N 和M,表示公墓的宽和长,因此这个矩形公墓共有(N+1) ×(M+1)个格点,左下角的坐标为(0, 0),右上角的坐标为(N, M)。第二行包含一个正整数W,表示公墓中常青树的个数。第三行起共W 行,每行包含两个用空格分隔的非负整数xi和yi,表示一棵常青树的坐标。输入保证没有两棵常青树拥有相同的坐标。最后一行包含一个正整数k,意义如题目所示。

Output

包含一个非负整数,表示这片公墓中所有墓地的虔诚度总和。为了方便起见,答案对2,147,483,648 取模。

Sample Input

5 6
13
0 2
0 3
1 2
1 3
2 0
2 1
2 4
2 5
2 6
3 2
3 3
4 3
5 2
2

Sample Output

6

HINT

图中,以墓地(2, 2)和(2, 3)为中心的十字架各有3个,即它们的虔诚度均为3。其他墓地的虔诚度为0。

所有数据满足1 ≤ N, M ≤ 1,000,000,000,0 ≤ xi ≤ N,0 ≤ yi ≤ M,1 ≤ W ≤ 100,000, 1 ≤ k ≤ 10。存在50%的数据,满足1 ≤ k ≤ 2。存在25%的数据,满足1 ≤ W ≤ 10000。

注意:”恰好有k颗树“,这里的恰好不是有且只有,而是从>=k的树中恰好选k棵

Source

 

O(nm)暴力:sigma(C(l(x,y),k)*C(r(x,y),k)*C(u(x,y),k)*C(d(x,y),k))((x,y)是整个地图上任意一点,l,r,u,d分别表示点(x,y)正左/右/上/下方的点数,C表示组合数)。

优化:
先将所有常青树以y坐标为第一关键字,x坐标为第二关键字,从小到大排序。对于第i棵常青树,如果第i-1棵常青树的y坐标与它相同,那么两棵常青树之间的空地的虔诚度之和为:

C(l(i-1),k)*C(r(i),k)*sigma(C(u(x),k)*C(d(x),k))(x介于第i-1~第i棵常青树之间。其中l(i-1)表示与第i-1棵常青树y坐标相同,且x坐标小于等于它的常青树数量,r(i)则表示与第i棵常青树y坐标相同,且x坐标大于等于它的常青树数量;u(x),d(x)表示点(x,y)正上/下方有多少常青树,C表示组合数)。

于是我们从左往右处理某行的某一个点时,要将树状数组中该点横坐标位置上的数进行修改

修改的值为就是现在的c(u[i],k)*c[d[i],k]减去原来的,也就是c(u[i],k)*c[d[i],k]-c(u[i]+1,k)*c[d[i]-1,k]

ps:网上有些人让ans自然溢出,可以加快速度。

  我只想说NO ZUO NO DIE。

别为了速度,WA了本可以过的点,那就得不偿失了!(因为我这样干过。。。)

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll mod=2147483648LL;
const int N=2e5+;
struct node{
int x,y;
}a[N];
int n,m,w,K,cnt,H[N],h[N],l[N],now[N];
ll ans,f[N],C[N][];//WA*1 看成50%的数据了
void getc(){
C[][]=;
for(int i=;i<=w;i++){
C[i][]=;
for(int j=;j<=min(K,i);j++){
C[i][j]=(C[i-][j]+C[i-][j-])%mod;
}
}
}
int find(int x){
int l=,r=cnt,mid;
while(l<=r){
mid=l+r>>;
if(H[mid]==x) return mid;
if(H[mid]>x) r=mid-;
else l=mid+;
}
}
int lowbit(int x){
return x&-x;
}
void updata(int x,ll val){
for(int i=x;i<=cnt;i+=lowbit(i)) f[i]=(f[i]+val)%mod;
}
ll query(int x){
ll res=;
for(int i=x;i;i-=lowbit(i)) res=(res+f[i])%mod;
return res;
}
bool cmp(const node &a,const node &b){
return a.y==b.y?a.x<b.x:a.y<b.y;
}
int main(){
scanf("%d%d%d",&m,&n,&w);
for(int i=;i<=w;i++){
scanf("%d%d",&a[i].x,&a[i].y);
H[++cnt]=a[i].x;H[++cnt]=a[i].y;
}
sort(H+,H+cnt+);
scanf("%d",&K);getc();
sort(a+,a+w+,cmp);
for(int i=;i<=w;i++) h[find(a[i].y)]++,l[find(a[i].x)]++;
int lf=;
for(int i=,d;i<=w;i++){
if(i>&&a[i-].y==a[i].y){
lf++;
ll t1=query(find(a[i].x)-)-query(find(a[i-].x));
ll t2=C[lf][K]*C[h[find(a[i].y)]-lf][K];//WA*2 忘了-1
ans+=t1*t2;
ans%=mod;
}
else lf=;
d=find(a[i].x);now[d]++;
ll change=(C[now[d]][K]*C[l[d]-now[d]][K]-C[now[d]-][K]*C[l[d]-now[d]+][K])%mod;
//c(u[i],k)*c[d[i],k]-c(u[i]+1,k)*c[d[i]-1,k]
updata(d,change);
}
if(ans<) ans+=mod;
printf("%d",(int)ans);
return ;
}

1227: [SDOI2009]虔诚的墓主人的更多相关文章

  1. BZOJ 1227: [SDOI2009]虔诚的墓主人

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

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

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

  3. BZOJ 1227 [SDOI2009]虔诚的墓主人 - 扫描线

    Solution 离散化 扫描线, 并用 $rest[i]$ 和 $cnt[i]$ 记录 第$i$列 总共有 $cnt[i]$棵常青树, 还有$rest[i]$ 没有被扫描到. 那么 第$i$ 列的方 ...

  4. 【以前的空间】bzoj 1227 [SDOI2009]虔诚的墓主人

    题解:hzw大神的博客说的很清楚嘛 http://hzwer.com/1941.html 朴素的做法就是每个点如果它不是墓地那么就可形成十字架的数量就是这个c(点左边的树的数量,k)*c(点右边的树的 ...

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

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

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

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

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

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

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

    [BZOJ1227][SDOI2009]虔诚的墓主人(线段树) 题面 BZOJ 洛谷 题解 显然发现答案就是对于每一个空位置,考虑上下左右各有多少棵树,然后就是这四个方向上树的数量中选\(K\)棵出来 ...

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

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

随机推荐

  1. vps安装wordpress遇到的问题(lnmp)

    1.要执行请求的操作,WordPress 需要访问您网页服务器的权限. 请输入您的 FTP 登录XXXX完美解决方法 因为在wordpress中新上传的插件的权限都是www用户的,而不是root或其他 ...

  2. Properties集合

    Map |--Hashtable |--Properties Properties集合特点: 1.该集合中的键和值都是字符串类型 2.集合中的数据可以保存在IO流中或者从IO流中获取数据. 通常该集合 ...

  3. 使用runloop阻塞线程的正确写法

    使用runloop阻塞线程的正确写法 runloop可以阻塞线程,等待其他线程执行后再执行. 比如: @implementation ViewController{    BOOL end;}…– ( ...

  4. Sea.Js的运行原理(转)

    1.CMD(Common Module Definition)规范 Sea.js采用了和Node相似的CMD规范,使用require.exports和module来组织模块.但Sea.js比起Node ...

  5. UVALive 4992 Jungle Outpost(半平面交)

    题意:给你n个塔(点)形成一个顺时针的凸包,敌人可以摧毁任何塔,摧毁后剩下的塔再组成凸包 在开始的凸包内选一点为主塔,保证敌人摧毁尽量多塔时主塔都还在现在的凸包内,求出最多摧毁的塔 题解:这题关键就是 ...

  6. zepto学习之路--核心函数$()的实现

    $()可以说是jquery的精华了,为dom操作带来了极大的灵活和方便.zepto号称“移动版的jquery”,那么它是怎么来实现这个核心函数呢?我们来详细探讨下. 1.首先,我们看下zepto中它是 ...

  7. java工程开发之图形化界面之(第三课)

    上面我们讲述了通过JOptionPane进行文本I/O操作,一个是通过JOptionPane来获取的 参考链接:http://blog.sina.com.cn/s/blog_993d254201013 ...

  8. Android, BaseAdapter 处理大数据量时的优化

    Android优化 最常见的就是ListView, Gallery, GridView, ViewPager 的大数据优化  图片优化  访问网络的优化优化的原则: 数据延迟加载  分批加载  本地缓 ...

  9. CentOS查询 杀死进程

    ps aux | grep XXX 查询进程 ps a 显示现行终端机下的所有程序,包括其他用户的程序. ps -A 显示所有程序. ps c 列出程序时,显示每个程序真正的指令名称,而不包含路径,参 ...

  10. CharSequence的getText()与String的getString()(转)

    CharSequence的getText()与String的getString()『Android系列七』 曾经在学习中碰见两种获取常量的方式: CharSequence chrs = getText ...