[2018.10.20 T2] 面包

暂无链接

面包

【问题描述】

小 X 有一块$N\ast M$的巨大面包，上面有$K$个葡萄干，第$i$个葡萄干的位置在第$A_i$行$B_i$列，任意两个葡萄干都不在同一位置。
现在小 X 的朋友小 C 来了，小 X 要从这块面包上切下呈矩形的一块给小 C。但是如果切下来的面包上一个葡萄干都没有，小 X 会很没面子，于是小 X 的一个合法的切割方案要满足切下的部分是含有至少一个葡萄干的矩形。现在小 X 想知道：所有合法的切割方案切下的部分包含葡萄干的个数的方差是多少？为了方便，请输出答案对$998,244,353$取模的结果，保证数据是在一定范围内随机生成的。

【输入格式】

第一行，三个正整数$N，M$和$K$ 。
接下来$K$行，每行两个整数$A_i$和$B_i$ 。

【输出格式】

一行，包含一个非负整数，表示答案对$998,244,353$取模的结果。
注意，你的答案应该在$[0,998,244,353)$的范围内。

【输入输出样例 1】

bread. in

1 2 2
1 1
1 2

bread.out

887328314

【输入输出样例说明 1】

有三种合法方案，其中两种方案包含一个葡萄干，一种方案包含两个葡萄干，所以方差为$\frac{9}{2}$，对$998244353$取模得到$887328314$

【输入输出样例 2】

见下发文件中的$bread2.in/bread2.out$。

【输入输出样例 3】

见下发文件中的$bread3.in/bread3.out$。

【数据范围与约定】

对于所有的数据，$1\le N,M\le 10^9$；$1\le K\le min(N\ast M,2000)$；$1\le A_i\le N$；$1\le B_i\le M$。
对于每个子任务的特殊限制：
Subtask1(20pts): $N,M\le 10$；
Subtask2(25pts): $N,M\le 100$；
Subtask3(15pts): $K\le 3$；
Subtask4(40pts): 无特殊限制。

题解

按照我不知道的一个套路，我们先化简一波方差的式子，设$x_i$为第$i$个合法矩形中的点数，$n$为合法的矩形个数：
$$\begin{array}{rl}& \frac{\sum (x_i-\overline{x}{)}^2}{n}\\ =& \frac{\sum (x_i^2-2x_i\overline{x}+{\overline{x}}^2)}{n}\\ =& \frac{\sum x_i^2}{n}-\frac{2\sum x_i\overline{x}}{n}+\frac{\sum {\overline{x}}^2}{n}\\ =& \frac{\sum x_i^2}{n}-\frac{2n{\overline{x}}^2}{n}+\frac{n{\overline{x}}^2}{n}\\ =& \frac{\sum x_i^2}{n}-{\overline{x}}^2\end{array}$$

第二个东西非常好求，我们直接枚举每个点，看它对几个矩形有贡献，就求到了$\sum x_i$；第一个东西又要用到一个套路，$x_i^2$相当于一个矩形里的点对数（包括自己跟自己），我们只需要$O(k^2)$枚举点对，看它们对几个矩形产生了贡献，就能统计出$\sum x_i^2$。

那么问题来了，$n$怎么求？

我们可以先枚举右下角的行数，在那一行上从左到右加入在右下角上面的点，维护一个左上角坐标范围的单调栈，就能统计出有多少个矩形包含了至少一个点，大致如下图：

单调栈维护之后：

愉快地统计出$n$之后当然要愉快的$AC$啦。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int M=2005,mod=998244353;
struct sd{int x,y;}ptx[M],pty[M],sta[M];
int sum[M],tot,top,sum1,sum2,n,m,k;
bool cmpy(sd a,sd b){return a.y==b.y?a.x>b.x:a.y<b.y;}
bool cmpx(sd a,sd b){return a.x<b.x;}
ll power(ll x,ll p){ll r=1;for(;p;p>>=1,x=x*x%mod)if(p&1)r=r*x%mod;return r;}
void in(){scanf("%d%d%d",&n,&m,&k);for(int i=1;i<=k;++i)scanf("%d%d",&ptx[i].x,&ptx[i].y),pty[i]=ptx[i];}
void ac()
{
	for(int i=1;i<=k;++i)(sum1+=1ll*ptx[i].x*(n-ptx[i].x+1)%mod*ptx[i].y%mod*(m-ptx[i].y+1)%mod)%=mod;
	for(int i=1;i<=k;++i)for(int j=1;j<=k;++j)
	(sum2+=1ll*min(ptx[i].x,ptx[j].x)*(n-max(ptx[i].x,ptx[j].x)+1)%mod*min(ptx[i].y,ptx[j].y)%mod*(m-max(ptx[i].y,ptx[j].y)+1)%mod)%=mod;
	sort(ptx+1,ptx+1+k,cmpx);sort(pty+1,pty+1+k,cmpy);ptx[k+1].x=n+1;
	for(int i=k,j,tmp;i>=1;--i)if(ptx[i].x!=ptx[i+1].x)
	{
		for(top=0,tmp=0,j=1;j<=k;++j)
		{
			if(pty[j].x>ptx[i].x)continue;
			(tmp+=1ll*sum[top]*(pty[j].y-sta[top].y)%mod)%=mod;
			for(;pty[j].x>sta[top].x&⊤--top);
			sta[++top]=pty[j];
			sum[top]=(sum[top-1]+1ll*pty[j].x*(pty[j].y-sta[top-1].y)%mod)%mod;
		}
		(tot+=1ll*(tmp+1ll*sum[top]*(m-sta[top].y+1)%mod)*(ptx[i+1].x-ptx[i].x)%mod)%=mod;
	}
	tot=power(tot,mod-2);
	printf("%lld",(1ll*sum2*tot%mod-1ll*sum1*tot%mod*sum1%mod*tot%mod+mod)%mod);
}
int main(){in(),ac();}