bzoj3216: 话旧2

大家好，我是SB

这道题是在看3215时找不到题面点进来的（然而并没有一眼看出这两题的差别）。

一开始想了一会，感觉就是一个组合数，然后最小值==0就有点像Catalan数，于是打开geogebra画图，得到了一个：Ans=Get(Y-X,Z)-Get(Y+X+1,Z-1)，反正证（gan）出（jue）他是对的，就狂WA。。。后来找的了→ ，发现自己傻了。。

然后用lucas定理+CRT（Chinese remainder theorem），算组合数，15144 MS A了

#include<iostream>
#include<cstdio>
#include<algorithm>
#define ll long long
#define O 19940417
#define P 2848631
#define _P 7
#define N 1000005
using namespace std;
ll n,k,x[N],y[N],Ans,Mx,w[N];
ll fac[P+1],inv[P+1],_fac[_P+1],_inv[_P+1];
bool cmp(ll a,ll b)
{
	return x[a]<x[b];
}
ll C(ll x,ll y)
{
	if (x>y) return 0;
	return fac[y]*inv[fac[y-x]]%P*inv[fac[x]]%P;
}
ll Lu(ll x,ll y)
{
	if (x==y||x==0) return 1;
	return C(x%P,y%P)*Lu(x/P,y/P)%P;
}
ll _C(ll x,ll y)
{
	if (x>y) return 0;
	return _fac[y]*_inv[_fac[y-x]]%_P*_inv[_fac[x]]%_P;
}
ll _Lu(ll x,ll y)
{
	if (x==y||x==0) return 1;
	return _C(x%_P,y%_P)*_Lu(x/_P,y/_P)%_P;
}
ll Get(ll x,ll y)
{
	if (x>y) return 0;
	if ((x&1)^(y&1)) return 0;
	ll p=(y-x)/2,q=y-p;
	ll t=Lu(p,q+p),_t=_Lu(p,q+p);
	ll tmp=(t*_P%O*inv[_P]+_t*P%O*_inv[P%_P])%O;
	//cout<<(tmp+O)%O<<endl;
	return tmp;
}
int main()
{
	
	
	scanf("%lld%lld",&n,&k);
	for (ll i=1;i<=k;i++)
		scanf("%lld%lld",&x[i],&y[i]),w[i]=i;
	x[0]=y[0]=0;w[0]=0;
	x[k+1]=n;y[k+1]=0;w[k+1]=k+1;
	sort(w+1,w+k+1,cmp);
	
	inv[1]=_inv[1]=1;
	for (ll i=2;i<P;i++)inv[i]=(-P/i)*inv[P%i]%P;
	for (ll i=2;i<_P;i++)_inv[i]=(-_P/i)*_inv[_P%i]%_P;
	fac[0]=_fac[0]=1;
	for (ll i=1;i<P;i++)fac[i]=i*fac[i-1]%P;
	for (ll i=1;i<_P;i++)_fac[i]=i*_fac[i-1]%_P;
	
	Ans=1;Mx=0;
	for (ll i=0;i<=k;i++)
	{
		ll X=y[w[i]],Y=y[w[i+1]],Z=x[w[i+1]]-x[w[i]];
		if (X>Y) swap(X,Y);
		if (Z) Ans=Ans*(Get(Y-X,Z)-Get(Y+X+2,Z))%O;
		Mx=max(Mx,Y+(Z-Y+X)/2);
	}
	printf("%lld %lld\n",(Ans+O)%O,Mx);
	return 0;
}