NOIP2020 退役记

NOIP2020 退役记

把原来写的东西全删了。毕竟留下了不少遗憾，也让我自己失望了吧。

12:30 想到 T3 怎么做，写完后发现第二个样例过不去，也就爆零了。

晚上回家理顺了一遍思路，写了半个小时不到就写完并在 oitiku 上 AC 了。我刷了不少 AGC 和 CF 的构造题，原本只是准备省选和 NOI 写的，结果 NOIP 考到了。我以为我能暴切 T3 考到省队线，结果被 T3 送退役了。

估分 \(80\sim 100+84\sim 100+0+60\sim 80=224\sim 280\)。不过应该省选还是会去的，但是翻盘的可能性不大了。

不过我也努力过了，这次的失误也是我想题较慢的缘故。若是把这道 T3 放到省选或 NOI 上，我感觉自己能现场写出来的。太遗憾了，实在是太遗憾了。

要是 T3 一点都不会做该多好，这样我就可以安心退役了。我真的心有不甘啊啊啊啊啊啊。

唉，在绝望中给人希望是一件多么残酷的事情。

T3 做法：

考虑分治，将颜色编号 \(\le mid\) 和 \(>mid\) 分为两个集合。这样问题转化为只用若干个长度为 \(m\) 的 \(01\) 串，将要求每个串全部 \(0/1\)。

接着考虑两个串 \(S,T\)。可以知道，\(cnt_{S,0/1}+cnt_{T,0/1}\) 总有一个 \(\ge m\)。那么，我们可以通过如下构造，得到一个串全部 \(0/1\)。

设空串为 \(now\)。取 \(S,T\) 的后 \(\lfloor \frac m2\rfloor\) 位，放到 \(now\) 中。如果是 \(0\) 放到 \(S\)，\(1\) 放到 \(T\)，那么可以得到一个长为 \(\lfloor \frac m2\rfloor\) 的串全部 \(0/1\)。若是 \(S\)，再将 \(S\) 全部放到 \(now\)。此时 \(S\) 为空串，\(T\) 没有性质，\(now\) 前 \(\lfloor \frac m2\rfloor\) 位相等。此处用了 \(3m\) 次操作。

对于串 \(now\) 和 \(T\)，我们取 \(now\) 的后 \(\lceil \frac m2\rceil\) 位，\(T\) 的后 \(\lfloor \frac m2\rfloor\) 位，执行类似刚才的操作。此时 \(now\) 的后 \(\lceil \frac m2\rceil\) 位相等。此处用了 \(2m\) 次操作。

现在讨论 \(now\) 前 \(\lfloor \frac m2\rfloor\) 位和后 \(\lceil \frac m2\rceil\) 位是否相等。若相等，那么我们已经构造了出来。否则，将 \(now\) 的后 \(\lceil \frac m2\rceil\) 位放到 \(S\)，然后对 \(T\) 执行刚才的过程。如果是 \(0\) 放到 \(S\)，\(1\) 放到 \(now\)。这样 \(S/now\) 一定有一个串符合条件，用了 \(\frac 32m\) 次操作。

上界大概是 \(O(F(n)\times \frac {13}{2}m)\)，\(F(n)=\sum (len_{rt}-1)\)，\(rt\in\) 长度为 \(n\) 的线段树。\(F(50)=237\)，可以放心写，不怎么卡常数。

#include <bits/stdc++.h>
using namespace std;
const int maxn=55;
const int maxa=820005;
int n,m,x_[maxa],y_[maxa],sz,now;
vector<int> g[maxn];
inline void add(int x,int y)
{
	x_[++sz]=x,y_[sz]=y;
	g[y].push_back(g[x].back());
	g[x].pop_back();
}
void solve(int l,int r,vector<int> &id)
{
	if(id.size()==1) return;
	int mid=(l+r)>>1;
	vector<int> A,B; A.clear(),B.clear();
	while(!id.empty())
	{
		int s=id.back(); id.pop_back();
		int t=id.back(); id.pop_back();
		// s, t -> now
		for(int i=0;i<m/2;i++) add(s,now),add(t,now);
		int pa=0,pb=0;
		for(int i=0;i<g[now].size();i++)
			if(g[now][i]<=mid) pa++;
			else pb++;
		int op1=(pa>=m/2)?0:1;
		while(!g[now].empty())
		{
			if(((op1==0 && g[now].back()<=mid) || (op1==1 && g[now].back()>mid)) && g[s].size()<m) add(now,s);
			else add(now,t);
		}
		for(int i=0;i<m;i++) add(s,now);
		for(int i=0;i<m/2;i++) add(t,s);
		for(int i=0;i<(m+1)/2;i++) add(now,s);
		int qa=0,qb=0;
		for(int i=0;i<g[s].size();i++)
			if(g[s][i]<=mid) qa++;
			else qb++;
		int op2=(qa>=(m+1)/2)?0:1;
		while(!g[s].empty())
		{
			if(((op2==0 && g[s].back()<=mid) || (op2==1 && g[s].back()>mid)) && g[now].size()<m) add(s,now);
			else add(s,t);
		}
		if((g[now][0]<=mid)==(g[now].back()<=mid))
		{
			if(g[now].back()<=mid) A.push_back(now);
			else B.push_back(now);
			if(id.empty())
			{
				if(g[t].back()<=mid) A.push_back(t);
				else B.push_back(t);
			}
			else id.push_back(t);
			now=s;
		}
		else
		{
			// op1 -> now, op2 -> s
			for(int i=0;i<(m+1)/2;i++) add(now,s);
			int x=0,y=0;
			for(int i=0;i<g[t].size();i++)
				if(g[t][i]<=mid) x++;
				else y++;
			if(op2==0)
			{
				if(x+g[s].size()>=m)
				{
					while(!g[t].empty())
					{
						if(g[t].back()<=mid && g[s].size()<m) add(t,s);
						else add(t,now);
					}
					A.push_back(s);
					if(id.empty())
					{
						if(g[now].back()<=mid) A.push_back(now);
						else B.push_back(now);
					}
					else id.push_back(now);
					now=t;
				}
				else
				{
					while(!g[t].empty())
					{
						if(g[t].back()>mid && g[now].size()<m) add(t,now);
						else add(t,s);
					}
					B.push_back(now);
					if(id.empty())
					{
						if(g[s].back()<=mid) A.push_back(s);
						else B.push_back(s);
					}
					else id.push_back(s);
					now=t;
				}
			}
			else
			{
				if(y+g[s].size()>=m)
				{
					while(!g[t].empty())
					{
						if(g[t].back()>mid && g[s].size()<m) add(t,s);
						else add(t,now);
					}
					B.push_back(s);
					if(id.empty())
					{
						if(g[now].back()<=mid) A.push_back(now);
						else B.push_back(now);
					}
					else id.push_back(now);
					now=t;
				}
				else
				{
					while(!g[t].empty())
					{
						if(g[t].back()<=mid && g[now].size()<m) add(t,now);
						else add(t,s);
					}
					A.push_back(now);
					if(id.empty())
					{
						if(g[s].back()<=mid) A.push_back(s);
						else B.push_back(s);
					}
					else id.push_back(s);
					now=t;
				}
			}
		}
	}
	solve(l,mid,A),solve(mid+1,r,B);
}
int main()
{
	freopen("ball.in","r",stdin);
	freopen("ball.out","w",stdout);
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			int x; scanf("%d",&x);
			g[i].push_back(x);
		}
	vector<int> id; id.clear();
	for(int i=1;i<=n;i++) id.push_back(i);
	now=n+1,solve(1,n,id);
	printf("%d\n",sz);
	for(int i=1;i<=sz;i++)
		printf("%d %d\n",x_[i],y_[i]);
	return 0;
}