2021牛客暑期多校训练营6 C题: Delete Edges

C题: Delete Edges

原题链接:https://ac.nowcoder.com/acm/contest/11257/C

题目大意

有一张 $n(n\le 2000)$ 个点的完全图，你可以进行如下操作:

1. 每次选取 $3$ 个不同的点 $x,y,z$ 构成的环( $3$ 边原本存在);
2. 删除这个环中的边 $(x,y)(x,z)(y,z)$ ;

求一种操作方案，使得若干次操作结束后剩余边数少于 $n$ 。

题解

大胆假设，大胆猜测

一张 $n$ 个点的完全图有 $\frac{n(n-1)}{2}$ 条边，若我们已有一种贪心方案，其对于 $n$ 个点的完全图的操作次数为 $f(n)$ ，则 $f(n)$ 应满足:
$$\begin{array}{rl}\frac{n(n-1)}{2}(总边数)-3f(n)(删去总边数)& <n(目的剩余边数)\\ -3f(n)& <\frac{2n-n(n-1)}{2}\\ f(n)& >\frac{n^2-3n}{6}\end{array}$$

我们猜测结论并证明(怎么才能在比赛时猜出来？真巧，我也不知道)。
我们假设 $f(n)$ 每次所选的 $x,y,z$ 满足 $(1\le )x<y<z(\le n)$ , $x+y+z\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}n)$ 然后我们构造其递推式。
显然， $f(n)$ 中每一种方案都对应了 $f(n+3)$ 中一种方案，因为 $x+y+z=n$ 可以映射到 $(x+1)+(y+1)+(z+1)\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}n+3)$ ， $x+y+z=2n$ 可以映射到 $(x+2)+(y+2)+(z+2)\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}n+3)$ 。
(原 $f(n)$ 中不存在 $x+y+z=3n$ ，因为 $x=y=z=n$ ，不合法)
注意，此时 $(x+1)>1,(z+2)<n+3$，即 $f(n+3)$ 中最小值 $x=1$ 和最大数字$z=n+3$的情况未纳入计算，我们枚举这些情况:

下表为 $x=1$ 时，部分满足条件的 $y,z$ 对应值(条件: $x+y+z=n+3$ ):

x	1	1	1	…	1
y	2	3	4	…	$\lfloor \frac{n+2}{2}\rfloor$(奇$n$) $\frac{n}{2}$(偶$n$)
z	n	n-1	n-2	…	$\lceil \frac{n+2}{2}\rceil$(奇$n$) $\frac{n+4}{2}$(偶$n$)

共 $\lfloor \frac{n-1}{2}\rfloor$ 组(不考虑奇偶性)
(注意极限，若继续枚举则 $y\ge z$ 不合法)

下表为 $z=n+3$ 时，部分满足条件的 $x,y$ 对应值(条件: $x+y+z=2(n+3)$ )

z	n+3	n+3	n+3	…	n+3
x	1	2	3	…	$\frac{n+1}{2}$(奇$n$) $\lfloor \frac{n+3}{2}\rfloor$(偶$n$)
y	n+2	n+1	n	…	$\frac{n+5}{2}$(奇$n$) $\lceil \frac{n+3}{2}\rceil$(偶$n$)

共 $\lfloor \frac{n}{2}\rfloor +1$ 组(不考虑奇偶性)
(注意极限，若继续枚举则 $x\ge y$ 不合法)

综上所述， $f(n)$ 转移式如下:
$$\begin{array}{rl}f(n+3)& =f(n)+\lfloor \frac{n-1}{2}\rfloor +\lfloor \frac{n}{2}\rfloor +1\\ f(n+3)& =f(n)+\{\begin{array}{rc}\frac{n-1}{2}+\frac{n-1}{2}+1(奇n)& \\ & \\ \frac{n-2}{2}+\frac{n}{2}+1(偶n)& \end{array}\\ f(n+3)& =f(n)+n(合并)\end{array}$$

因为转移式不方便直接使用，我们有以下两种方案:

1. 采用矩阵快速幂对线代进行加速
2. 转化为通项式

在文末简单介绍了矩阵快速幂的构造，供参考

因为题目数据范围 $n\le 2000$ ，矩阵快速幂效果并不显著，且不必要，在此讲述通项式转化法:

若 $n$ 为 $3$ 的倍数，观察:$f(3)=1\stackrel{+3}{⟶}f(6)=4\stackrel{+6}{⟶}f(9)=10\stackrel{+9}{⟶}...$

易得 $f(n)$ 为 $f(3)=1$ 加等差数列 $3,6,9,...,n-3$ ，通项式:
$$f(n)=1+\frac{\frac{n-3}{3}n}{2}=\frac{n^2-3n+6}{6}$$

若 $n\%3=1$ ，观察:$f(4)=1\stackrel{+4}{⟶}f(7)=5\stackrel{+7}{⟶}f(10)=12\stackrel{+10}{⟶}...$

易得 $f(n)$ 为 $f(4)=1$ 加等差数列 $4,7,10,...,n-3$ ，通项式:
$$f(n)=1+\frac{\frac{n-4}{3}(n+1)}{2}=\frac{n^2-3n+2}{6}$$

若 $n\%3=2$ ，观察:$f(5)=2\stackrel{+5}{⟶}f(8)=7\stackrel{+8}{⟶}f(11)=15\stackrel{+11}{⟶}...$

易得 $f(n)$ 为 $f(5)=2$ 加等差数列 $5,8,11,...,n-3$ ，通项式:
$$f(n)=2+\frac{\frac{n-5}{3}(n+2)}{2}=\frac{n^2-3n+2}{6}$$

合并后如下:
$$f(n)=\{\begin{array}{rc}\frac{n^2-3n+6}{6}& >\frac{n^2-3n}{6}(n\%3=0)\\ & \\ \frac{n^2-3n+2}{6}& >\frac{n^2-3n}{6}(n\%3!=0)\end{array}$$

显然，该方案满足第一部分推算出的要求，合理。
输出方案时，枚举 $x,y$ ，然后计算输出对应的 $z$ 即可。

参考代码

#include<bits/stdc++.h>
#define For(i,n,m) for(int i=n;i<=m;i++)
using namespace std;

int n;

int main()
{
	cin>>n;
	if(n==3){
		puts("1");
		puts("1 2 3");
	}
	else{
		printf("%d\n",n%3?(n*n-3*n+2)/6:(n*n-3*n+6)/6)//方案总操作数
		int k;
		For(i,1,n){
			For(j,i+1,n){
				k=(2*n-i-j-1)%n+1;//和为n或2n(3n需i=j=k=n,不合法),最大2n,减法计算后取模即可。注意先-1再+1,使得1<=k<=n
				if(k>j)printf("%d %d %d\n",i,j,k);//若k<j为非法方案则不输出
			}
		}
	}
	return 0;
}

附:矩阵快速幂计算方案操作数

我们首先设原矩阵为 $\left[\begin{array}{c}f(n)\\ f(n-1)\\ f(n-2)\end{array}\right]$ ，目标矩阵为 $\left[\begin{array}{c}f(n+1)\\ f(n)\\ f(n-1)\end{array}\right]$ 初步操作式为:
$$\left[\begin{array}{c}f(n+1)\\ f(n)\\ f(n-1)\end{array}\right]=\left[\begin{array}{c}???\\ ???\\ ???\end{array}\right]\left[\begin{array}{c}f(n)\\ f(n-1)\\ f(n-2)\end{array}\right]$$

倒推出 $?$ 部分，发现计算 $f(n+1)=f(n-2)+n-2$ 需要常数 $n-2$ ，添加一项:
$$\left[\begin{array}{c}f(n+1)\\ f(n)\\ f(n-1)\\ n-1\end{array}\right]=\left[\begin{array}{c}0011\\ 1000\\ 0100\\ ????\end{array}\right]\left[\begin{array}{c}f(n)\\ f(n-1)\\ f(n-2)\\ n-2\end{array}\right]$$

发现此时需更新 $n-2$ 为 $n-1(n-2+1)$ ，再添加一个常数项 $1$ ，得最终转移式:
$$\left[\begin{array}{c}f(n+1)\\ f(n)\\ f(n-1)\\ n-1\\ 1\end{array}\right]=\left[\begin{array}{c}00110\\ 10000\\ 01000\\ 00011\\ 00001\end{array}\right]\left[\begin{array}{c}f(n)\\ f(n-1)\\ f(n-2)\\ n-2\\ 1\end{array}\right]$$