Solution -「CF 1391E」Pairs of Pairs

\(\mathcal{Description}\)

  Link.

  给定一个 \(n\) 个点 \(m\) 条边的无向图，在其上找到一条包括不少于 \(\lceil\frac{n}2\rceil\) 个结点的简单路径；或者将至少 \(\lceil\frac{n}2\rceil\) 个结点划分为若干二元组，使得任意两个不同二元组内四个结点的导出子图含有至多两条边。多组数据。

  \(n,\sum n\le5\times10^5\)，\(m,\sum m\le10^6\)。

\(\mathcal{Solution}\)

  去分开刚两个 NP 问题，请。

  这种给两个问题让你解其中一个的，显然以尝试求解一个失败为条件来求解另一个，不可能分开解的。比如本题，随便选一个起点 DFS 全图，若发现深度不小于 \(\lceil\frac{n}2\rceil\) 的 DFS 树上结点，回答第一问，结束。否则，问题转化为：一个图的 DFS 满足树深 \(<\lceil\frac{n}2\rceil\)，求解第二问。

  由于 DFS 树没有横叉边，所以在树上某一结点的左子树选一个，右子树选一个来组成一对，这两对间的导出子图显然不超过两条边。所以用 std::vector 记录子树内匹配剩下的结点，在父亲处与兄弟子树剩下的结点配对即可。可以证明在树深 \(<\lceil\frac{n}2\rceil\) 的情况下一定有解。

  复杂度 \(\mathcal O(\sum n+m)\)。

\(\mathcal{Code}\)

/* Clearink */

#include <cstdio>
#include <vector>
#include <cstdlib>
#include <assert.h>

inline int rint () {
	int x = 0, f = 1; char s = getchar ();
	for ( ; s < '0' || '9' < s; s = getchar () ) f = s == '-' ? -f : f;
	for ( ; '0' <= s && s <= '9'; s = getchar () ) x = x * 10 + ( s ^ '0' );
	return x * f;
}

template<typename Tp>
inline void wint ( Tp x ) {
	if ( x < 0 ) putchar ( '-' ), x = ~ x + 1;
	if ( 9 < x ) wint ( x / 10 );
	putchar ( x % 10 ^ '0' );
}

const int MAXN = 5e5, MAXM = 1e6;
int n, m, ecnt, head[MAXN + 5], p[MAXN + 5];
bool vis[MAXN + 5], solved;
std::vector<int> path, rest[MAXN + 5];
std::vector<std::pair<int, int> > pairs;

struct Edge { int to, nxt; } graph[MAXM * 2 + 5];

inline void link ( const int s, const int t ) {
	graph[++ ecnt] = { t, head[s] };
	head[s] = ecnt;
}

inline void clear () {
	ecnt = solved = 0, path.clear (), pairs.clear ();
	for ( int i = 1; i <= n; ++ i ) head[i] = vis[i] = 0, rest[i].clear ();
}

inline void printPath () {
	solved = true;
	printf ( "PATH\n%d", ( int ) path.size () );
	for ( int i = 0; i ^ path.size (); ++ i ) {
		printf ( "%c%d", i ? ' ' : '\n', path[i] );
	}
	putchar ( '\n' );
}

inline void printPairs () {
	solved = true;
	printf ( "PAIRING\n%d\n", ( int ) pairs.size () );
	for ( auto p: pairs ) printf ( "%d %d\n", p.first, p.second );
}

inline void solve ( const int u ) {
	if ( solved ) return ;
	vis[u] = true, path.push_back ( u ), p[u] = u;
	if ( ( int ) path.size () << 1 >= n ) printPath ();
	for ( int i = head[u], v; i; i = graph[i].nxt ) {
		if ( !vis[v = graph[i].to] ) {
			solve ( v ), path.pop_back ();
			if ( solved ) return ;
			while ( !rest[p[u]].empty () && !rest[p[v]].empty () ) {
				pairs.push_back ( { rest[p[u]].back (), rest[p[v]].back () } );
				if ( ( int ) pairs.size () << 2 >= n ) return printPairs ();
				rest[p[u]].pop_back (), rest[p[v]].pop_back ();
			}
			p[u] = rest[p[u]].empty () ? p[v] : p[u];
		}
	}
	rest[p[u]].push_back ( u );
}

int main () {
	for ( int T = rint (); T --; ) {
		n = rint (), m = rint (), clear ();
		for ( int i = 1, u, v; i <= m; ++ i ) {
			u = rint (), v = rint ();
			link ( u, v ), link ( v, u );
		}
		solve ( 1 );
		assert ( solved );
	}
	return 0;
}

\(\mathcal{Details}\)

  std::vector 的操作能少压入就少压入啊……不知道怎么内存就挂掉了。