ZROI-Day2比赛解题报告

ZROIDay2-比赛解题报告

版权原因不提供题面信息

---

序

这几天作息有点鬼畜,虽然昨晚很晚睡但是早上精神还不错,看到题发现T1很友好？T2woc这暴力都好难打?T3多项式？！这样下去比赛会不会出现更多高端操作，恐怕凉凉

A

感谢出题人,暴力好打分又多,正解也不难想,这题基本上部分分都打了一遍

50pts

对于每一个炮将其所在的所在的交叉行（暂且这么说）\(O(1)\) 标记，然后\(O(N^2)\)遍历一遍统计就好了

70pts

核心思想是计算出放一个炮新增的贡献,即它能覆盖的点数减去已经覆盖的点数,最后\(N^2\)减去总和既是答案

拿这部分分还花了不少功夫,终于运用人类智慧找出一些规律,也就是对于一个炮\((x,y)\),它左斜行和右斜行能覆盖的点的个数,然后又发现对于一个确定的左斜行（即确定的\(x+y\)），可以通过\(x,y\)计算出与它有公共点的左斜行的\(x-y\)的相关信息，于是新增一个炮,在他所在的左斜行加上+1标记,加上左斜行覆盖点数,在枚举与左斜行有公共点的右斜行,假设这个右斜行已有标记,则贡献-1。对于右斜行也类似

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <cctype>
#include <queue>
#include <cmath>
#define ll long long
#define ri register int
#define ull unsigned loong long
const int maxn=100005;
const int inf=0x7ffffffff;
template <class T>inline void read(T &x){
	x=0;int ne=0;char c;
	while(!isdigit(c=getchar()))ne=c=='-';
	x=c-48;
	while(isdigit(c=getchar()))x=(x<<3)+(x<<1)+c-48;
	x=ne?-x:x;return;
}
int n,m;
int p[maxn<<1],q[maxn<<1];
int cntp=0,cntq=0;
inline void solve_2(){
	int x,y;ll ans=0;
	while(m--){
      read(x),read(y);
      if(!p[x+y]){
      	p[x+y]=1;
      	cntp++;
      	if(x+y<=n){
      	  int k=x+y-2;//
      	  ans+=x+y-1;//ok
      	  for(ri i=-k;i<=k;i+=2){
      	  	ans-=q[i+n];
		  }
		}
      	else{
      	  int k=n*2-(x+y);
      	  ans+=k+1;
      	  for(ri i=(x+y)-n*2;i<=n*2-(x+y);i+=2){
      	  	ans-=q[i+n];
		  }
		}
	  }
	  if(!q[x-y+n]){
	  	q[x-y+n]=1;
	  	cntq++;
	  	if(x-y<=0){
	  		int k=x-y+n-1;
	  		ans+=n+(x-y);
	  		for(ri i=n+1-k;i<=n+1+k;i+=2){
	  			ans-=p[i];
			}
		}
		else{
			int k=x-y-n+1;
			ans+=n-(x-y);
			for(ri i=n+1+k;i<=n+1-k;i+=2){
				ans-=p[i];
			}
		}
	  }
	}
	printf("%lld\n",1ll*n*n-ans);
}
int main(){
	int x,y;
	ll ans=0;
	freopen("dat.in","r",stdin);
	freopen("bf.out","w",stdout);
	read(n),read(m);
    solve_2();
	return 0;
}

100pts

发现每个左斜行能确定的右斜行的x-y范围是连续的的奇数或偶数,于是用线段树维护标记和区间和,复杂度$O( m $ $log $ \(N)\)

然后鬼畜的是刚码完大样例过不了,以为是奇偶数搞错,魔改了半天后发现有一个判定没写到循环里。。。然后又魔改还是不对,于是开始对拍手动gdb调试,发现有一颗线段树操作的上限因为是x+y要设成\(2n\),查完这个错后据考试结束还有不到半小时,真TM刺激

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <cctype>
#include <ctime>
#include <queue>
#include <cmath>
#define ll long long
#define ri register int
#define ull unsigned loong long
const int maxn=200015;
const int inf=0x7fffffff;
template <class T>inline void read(T &x){
	x=0;int ne=0;char c;
	while(!isdigit(c=getchar()))ne=c=='-';
	x=c-48;
	while(isdigit(c=getchar()))x=(x<<3)+(x<<1)+c-48;
	x=ne?-x:x;return;
}
int n,m,N;
bool p[maxn<<1],q[maxn<<1];
int L,R,dta,t;
struct Segment_Tree_1{
	int ns[maxn<<2],odds[maxn<<2];//fff是没用的
	void update(int now,int l,int r,int fff){//printf("%d****%d\n",l,r);
		if(l==r){
			if(l&1)odds[now]++;
			else ns[now]++;
			return ;
		}
		int mid=(l+r)>>1;
		if(t<=mid)update(now<<1,l,mid,fff);
		else update(now<<1|1,mid+1,r,fff);
		odds[now]=odds[now<<1]+odds[now<<1|1];
		ns[now]=ns[now<<1]+ns[now<<1|1];
		return ;
	}
	int odd_query(int now,int l,int r){
		if(L<=l&&r<=R){
			return odds[now];
		}
		int mid=(l+r)>>1,ans=0;
		if(L<=mid)ans+=odd_query(now<<1,l,mid);
		if(mid<R)ans+=odd_query(now<<1|1,mid+1,r);
		return ans;
	}
	int n_query(int now,int l,int r){
		if(L<=l&&r<=R){
			return ns[now];
		}
		int mid=(l+r)>>1,ans=0;
		if(L<=mid)ans+=n_query(now<<1,l,mid);
		if(mid<R)ans+=n_query(now<<1|1,mid+1,r);
		return ans;
	}
}P;
struct Segment_Tree_2{
	int ns[maxn<<1],odds[maxn<<1];
	void update(int now,int l,int r,int fff){
		if(l==r){
			if(l&1)odds[now]++;
			else ns[now]++;
			return ;
		}
		int mid=(l+r)>>1;
		if(t<=mid)update(now<<1,l,mid,fff);
		else update(now<<1|1,mid+1,r,fff);
		odds[now]=odds[now<<1]+odds[now<<1|1];
		ns[now]=ns[now<<1]+ns[now<<1|1];
		return ;
	}
	int odd_query(int now,int l,int r){
		if(L<=l&&r<=R){
			return odds[now];
		}
		int mid=(l+r)>>1,ans=0;
		if(L<=mid)ans+=odd_query(now<<1,l,mid);
		if(mid<R)ans+=odd_query(now<<1|1,mid+1,r);
		return ans;
	}
	int n_query(int now,int l,int r){
		if(L<=l&&r<=R){
			return ns[now];
		}
		int mid=(l+r)>>1,ans=0;
		if(L<=mid)ans+=n_query(now<<1,l,mid);
		if(mid<R)ans+=n_query(now<<1|1,mid+1,r);
		return ans;
	}
}Q1,Q2;
inline void solve(){
	int x,y;ll ans=0;
	while(m--){
      read(x),read(y);
      if(!p[x+y]){
      	t=x+y;
      	P.update(1,1,N,0);
		p[x+y]=1;
      	if(x+y<=n){
      	  int k=x+y-2;
      	  ans+=x+y-1;
      	  int lxl=-k,rr=k;
      	  if(lxl==rr){
      	  	L=R=0;
			if(q[n])ans--;
		  }
      	  else {
      	    L=0,R=abs(lxl);
			if(R&1)ans-=Q1.odd_query(1,0,n);
			else ans-=Q1.n_query(1,0,n);
		    L=1,R=rr;
		    if(R&1)ans-=Q2.odd_query(1,1,n);
		    else ans-=Q2.n_query(1,1,n);
		  }
		}
      	else{
      	  int k=n*2-(x+y);
      	  ans+=k+1;
      	  int lxl=(x+y)-n*2,rr=n*2-(x+y);
      	  //printf("2--%d ",ans);
      	  if(lxl==rr){
      	  	L=R=0;
			if(q[n])ans--;//ans-=Q1.n_query(1,0,n);
		  }
      	  else {
      	    L=0,R=abs(lxl);
			if(R&1)ans-=Q1.odd_query(1,0,n);
			else ans-=Q1.n_query(1,0,n);
		    L=1,R=rr;
		    if(R&1)ans-=Q2.odd_query(1,1,n);
		    else ans-=Q2.n_query(1,1,n);
		  }
		}
	  }
	  if(!q[x-y+n]){
	  	t=x-y;
	  	if(t<=0){
	  		t=abs(t);
		    Q1.update(1,0,n,1);
		}
		else {
			Q2.update(1,1,n,1);
		}
	  	q[x-y+n]=1;
	  	if(x-y<=0){
	  		int k=x-y+n-1;
	  		ans+=n+(x-y);
	  		L=n+1-k,R=n+1+k;
	  		//printf("3--%d ",ans);
	  		if(L&1)ans-=P.odd_query(1,1,N);
	  		else ans-=P.n_query(1,1,N);
		}
		else{
			int k=x-y-n+1;
			ans+=n-(x-y);
			L=n+1+k,R=n+1-k;
			//printf("%d %d\n",L,R);
			//printf("4--%d ",ans);
			if(R&1)ans-=P.odd_query(1,1,N);
			else {
			ans-=P.n_query(1,1,N);
	 	    }
		 }
	   }
	}
	printf("%lld\n",1ll*n*n-ans);
}
int main(){
	int x,y;
	ll ans=0;
	double st=clock();
	freopen("dat.in","r",stdin);
	freopen("std.out","w",stdout);
	read(n),read(m);N=n*2;
    solve();
    double ed=clock();
    printf("%lf\n",ed-st);
	return 0;
}

B

毒瘤期望,n=4的暴力枚举都及其毒瘤,没时间打

题解暂时没搞懂

C

多项式算了吧，考场上暴力模拟感谢出题人拿了40pts,题解推了一大波东西然后什么阶乘卷积。。。

贴一份老师的std

#include<stdio.h>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<cmath>
#include<string>
#include<assert.h>
#define rep(i,j,k) for(int i=(int)j;i<=(int)k;i++)
#define per(i,j,k) for(int i=(int)j;i>=(int)k;i--)
#define pii pair<int,int>
#define fi first
#define se second
#define pb push_back
using namespace std;
typedef long long LL;
const int P=998244353;
const int G=3;
const int N=310000;
inline int Pow(int a,int b){
	int c=1;
	for(;b;b>>=1,a=a*1ll*a%P)if(b&1)c=c*1ll*a%P;
	return c;
}
int w[2][N];
int rev[N];
inline void fft(int *a,int n,int ff){
	rep(i,0,n-1)if(i<rev[i])swap(a[i],a[rev[i]]);
	for(int i=1;i<n;i<<=1)
	for(int j=0,t=n/(i<<1);j<n;j+=(i<<1))
	for(int k=0,l=0;k<i;k++,l+=t){
		int x=a[i+j+k]*1ll*w[ff][l]%P;
		int y=a[j+k];
		a[j+k]=(x+y)%P;
		a[i+j+k]=(y+P-x)%P;
	}
	if(ff==1){
		int v=Pow(n,P-2);
		rep(i,0,n-1)a[i]=a[i]*1ll*v%P;
	}
}
inline void initfft(int n){
	rep(i,0,n-1){
		int x=i;
		int y=0;
		for(int k=1;k<n;k<<=1,x>>=1)(y<<=1)|=(x&1);
		rev[i]=y;
	}
	w[0][0]=w[1][0]=1;
	int V=Pow(G,(P-1)/n);
	int VV=Pow(V,P-2);
	rep(i,1,n-1){
		w[0][i]=w[0][i-1]*1ll*V%P;
		w[1][i]=w[1][i-1]*1ll*VV%P;
	}
}
int a[N],n;
int fac[N],inv[N];
int ans[N];
inline void init(int n){
	fac[0]=1;rep(i,1,n)fac[i]=fac[i-1]*1ll*i%P;
	inv[n]=Pow(fac[n],P-2);per(i,n-1,0)inv[i]=inv[i+1]*1ll*(i+1)%P;
}
int main(){
	init(200000);
	scanf("%d",&n);
	assert(1<=n&&n<=100000);
	rep(i,0,n-1){
		scanf("%d",&a[i]);
		assert(0<=a[i]&&a[i]<P);
	}
	rep(i,0,n-1)a[i]=a[i]*1ll*fac[i]%P;

	initfft(1<<18);
	fft(a,1<<18,0);
	rep(i,0,(1<<18)-1)a[i]=a[i]*1ll*a[i]%P;
	fft(a,1<<18,1);

	rep(d,0,n-1){
		ans[d]=a[n-1+d];
		ans[d]=ans[d]*1ll*Pow(2,d)%P;
		ans[d]=ans[d]*1ll*inv[d]%P;
	}
	rep(i,0,n-1)printf("%d ",ans[i]);
	return 0;
}