bzoj4558

真是一道非常excited的题目啊…JLOI有毒

题目大意:给一个(N+1)*(M+1)的网格图,格点坐标为(0~N,0~M),现在挖去了K个点,求剩下多少个正方形(需要注意的是正方形可以是斜着的,多斜都可以)

N,M<=10^6,K<=2*10^3。

首先我们发现有一个非常感人的K=0部分分…

我们考虑K=0怎么做。

对于一个形如这样的正方形,我们叫它(a,b)正方形好了。

我们可以很容易地发现一个(a,b)正方形实际上要占下(a+b)*(a+b)这么大一块网格。

然后我们考虑a+b的大小,这样a就是[0,a+b]这么大,这样就可以得到一个答案。

代码如下:

现在我们发现,有了这些障碍物,我们只要能求出总共的正方形个数、经过一个障碍点的正方形个数、经过两个障碍点的正方形个数、经过三个点的、经过四个点的即可。

经过三个和经过四个直接二分查找一下显然是trivial的,经过两个点的要考虑是作为边往两侧延伸和作为对角线的情况,也比较trivial。

总共的正方形个数我们已经求出来了,现在我们就要考虑经过某一个障碍点的正方形个数。

对于一个点和它相关的只有四个属性,u,d,l,r对吧。

首先我们考虑直的正方形,即(0,x)或(x,0)正方形,因为这类正方形容易被重复统计。

容易发现这类正方形个数为min(u,l)+min(u,r)+min(l,d)+min(d,r)。

其它的正方形显然都是在四个象限中某两个相邻象限的。

为了简化起见,我们先考虑l,r,d这一象限的。

还是一样,设正方形为(a,b)正方形,我们枚举a+b,假设a+b=c。

因为正方形不是直的,所以a,b≠0。

现在我们考虑求出a的取值范围。

容易发现a<=r,a<=c-1,a>=1,a>=c-l(由于b<=l)。

那么我们可以列出一个形如这样的式子来计算:

这样显然不够优秀,我们可以人工分类讨论一下…

额其实注意到当r=c-1时c=r+1,当c-l=1时c=l+1,那么min和max的两个“分界点”是l+1和r+1,在分界点中间显然都是一些一次函数,那么就都是等差数列,于是我们就可以避免人工分类讨论了。

有了这个函数calc(l,r,d),那么calc(u,d,l,r)显然就等于

一些奇怪的细节详见代码…

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <algorithm>
#include <math.h>
using namespace std;
typedef long long ll;
#define MOD 100000007
ll calc(ll l,ll r,ll d)
{
if(!l||!r||!d) return 0;
ll ans=0;
ll upp=min(l+r,d);
ll ps[3]={l+1,r+1,upp};
sort(ps,ps+3);
ll cl=1;
for(int i=0;i<3;i++)
{
ll cr=ps[i];
if(cr>upp) break;
if(cr<2||cl==cr) continue;
++cl;
ll vl=min(r,cl-1)-max(cl-l,1LL)+1;
ll vr=min(r,cr-1)-max(cr-l,1LL)+1;
ans+=(vl+vr)*(cr-cl+1)/2;
ans%=MOD;
cl=cr;
}
return ans;
}
ll calc(ll u,ll d,ll l,ll r)
{
return calc(u,d,l)+calc(u,d,r)+calc(l,r,u)+calc(l,r,d)
+min(u,r)+min(u,l)+min(d,l)+min(d,r);
}
typedef pair<ll,ll> pll;
pll ps[233333];
#define X first
#define Y second
ll n,m,k,ans=0;
bool ok(pll a)
{
return a.X>=0&&a.X<n&&a.Y>=0&&a.Y<m;
}
ll tointt(double x)
{
if(fabs(x-ll(x+0.5))<1e-5) return x+0.5;
return -1;
}
double chk(double x,double y)
{
ll xx=tointt(x),yy=tointt(y);
if(xx>=0&&xx<n&&yy>=0&&yy<m) return 1;
return 0;
}
int main()
{
cin>>n>>m>>k; ++n; ++m;
ll cnt3=0,cnt4=0;
for(ll g=1;g<=n&&g<=m;g++) ans+=(n-g)%MOD*(m-g)%MOD*g%MOD, ans%=MOD;
for(int i=1;i<=k;i++)
{
ll x,y;
scanf("%lld%lld",&x,&y);
ans-=calc(x,n-1-x,y,m-1-y);
ans%=MOD;
ps[i]=pll(x,y);
}
sort(ps+1,ps+1+k);
for(int i=1;i<=k;i++)
{
for(int j=i+1;j<=k;j++)
{
do{
double mx=(ps[i].X+ps[j].X)/2.0,my=(ps[i].Y+ps[j].Y)/2.0;
double dx=ps[i].X-mx,dy=ps[i].Y-my;
if(chk(mx-dy,my+dx)&&chk(mx+dy,my-dx)) ans++;
}while(0);
for(int p=-1;p<=1;p+=2)
{
ll dx=ps[j].X-ps[i].X,dy=ps[j].Y-ps[i].Y;
pll n1=pll(ps[j].X-dy*p,ps[j].Y+dx*p);
pll n2=pll(ps[i].X-dy*p,ps[i].Y+dx*p);
if(ok(n1)&&ok(n2));else continue;
ans++;
int cp=0;
if(binary_search(ps+1,ps+1+k,n1)) ++cp;
if(binary_search(ps+1,ps+1+k,n2)) ++cp;
if(cp==1) cnt3++;
else if(cp==2) cnt3++, cnt4++;
}
}
}
ans-=cnt3/2; ans-=cnt4/4;
printf("%d\n",int(((ans%MOD)+MOD)%MOD));
}

JLOI2016 方的更多相关文章

  1. bzoj4558[JLoi2016]方 容斥+count

    4558: [JLoi2016]方 Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 452  Solved: 205[Submit][Status][D ...

  2. 4558: [JLoi2016]方

    4558: [JLoi2016]方 https://lydsy.com/JudgeOnline/problem.php?id=4558 分析: 容斥原理+各种神奇的计数. 如果没有被删除了的点的话,直 ...

  3. 【BZOJ 4558】 4558: [JLoi2016]方 (计数、容斥原理)

    未经博主同意不能转载 4558: [JLoi2016]方 Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 362  Solved: 162 Descri ...

  4. bzoj4558: [JLoi2016]方

    Description 上帝说,不要圆,要方,于是便有了这道题.由于我们应该方,而且最好能够尽量方,所以上帝派我们来找正方形 上帝把我们派到了一个有N行M列的方格图上,图上一共有(N+1)×(M+1) ...

  5. [BZOJ4558]:[JLoi2016]方(容斥+模拟)

    题目传送门 题目描述 上帝说,不要圆,要方,于是便有了这道题.由于我们应该方,而且最好能够尽量方,所以上帝派我们来找正方形上帝把我们派到了一个有N行M列的方格图上,图上一共有$(N+1)\times ...

  6. BZOJ.4558.[JLOI2016]方(计数 容斥)

    BZOJ 洛谷 图基本来自这儿. 看到这种计数问题考虑容斥.\(Ans=\) 没有限制的正方形个数 - 以\(i\)为顶点的正方形个数 + 以\(i,j\)为顶点的正方形个数 - 以\(i,j,k\) ...

  7. bzoj千题计划281:bzoj4558: [JLoi2016]方

    http://www.lydsy.com/JudgeOnline/problem.php?id=4558 容斥原理 全部的正方形-至少有一个点被删掉的+至少有两个点被删掉的-至少有3个点被删掉的+至少 ...

  8. 【BZOJ】4558: [JLoi2016]方

    [题意]给定有(n+1)*(m+1)个点的网格图,其中指定k个点不合法,求合法的正方形个数(四顶点合法). [算法]计数 [题解]斜着的正方形很麻烦,所以考虑每个斜正方形其外一定有正的外接正方形. 也 ...

  9. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

随机推荐

  1. tomcat中的JSP引擎

    jsp文件解析成class文件过程需要利用tomcat的jasper组件. Jasper是tomcat中使用的JSP引擎,在Tomcat 6中使用的是Jasper 2,相对于原来的版本作了不少的改进, ...

  2. NSFileManager和NSFileHandler的作用

    NSFileManager类可以管理文件系统中的文件和目录,也可以定位.判断是否存在.创建.拷贝.删除文件和目录,还可以获得文件和目录的信息: 对于读写文件,NSFileManager只可以读写字符串 ...

  3. 看苹果官方API

    command+shift+0会出现如下图 然后输入你想找的API 记得找带Reference这种标记的文档

  4. 【原】iOS:一种直接修改frame的某个属性的方法

    在iOS中view的frame属性使用地太频繁了,尤其是调UI的时候.我们知道,正常情况下我们无法对frame的某个属性(x,y,width,height等)进行单独修改,比如: someView.f ...

  5. HTML5设计网页熔岩灯导航(navigation bar)插件 已经加上完整源代码

    导航栏(navigation bar): 1.指位于页眉区域的,在页眉横幅图片上边或下边的一排水平导航按钮,它起着链接博客的各个页面的作用. 2.网页设计中不可缺少的部分,它是指通过一定的技术手段,为 ...

  6. MVC模式与三层架构和表示层

    1.MVC模式     - Model-View-Controller     - 模型-视图-控制器     - Model(模型)         > 模型分为业务模型,和数据模型     ...

  7. [MySQL Reference Manual] 7 备份和恢复

    7. 备份和恢复 本章主要会介绍: 1.备份的类型:逻辑备份,物理备份,全备和增量4种 2.创建备份的方法 3.还原方法,包括还原到时间点 4.备份计划,压缩和加密 5.表维护,恢复损坏的表 7. 备 ...

  8. 用struts和hibernate结合MVC层实例

    1.倒包hibernate11个包+sturts2 13个包 2.创建web.xml <?xml version="1.0" encoding="UTF-8&quo ...

  9. Java对象的序列化

    1.概念 序列化:把Java对象转换为字节序列的过程. 反序列化:把字节序列恢复为Java对象的过程. 2.用途 对象的序列化主要有两种用途: 1) 把对象的字节序列永久地保存到硬盘上,通常存放在一个 ...

  10. 定义类Human,具有若干属性和功能;定义其子类Man、Woman; 在主类Test中分别创建子类、父类和上转型对象,并测试其特性。

    package com.hanqi.test; public class Humen { private String speak; private String read; private Stri ...