[POI2012]Squarks

[POI2012]Squarks

题目大意：

设有\(n\)个互不相同的正整数\(\{X_1,X_2,...,X_n\}\)，任取两个\(X_i,X_j(i\ne j)\)，能算出\(X_i+X_j\)。
现在所有取法共\(\frac{n(n-1)}2\)个和，要你求出\(\{X_i\}\)的取值方案数，并求出所有方案的\(X_1,X_2,\ldots,X_n\)。

思路：

设\(\{X_i\}\)两两之和构成集合\(A\)，将\(A\)和\(X\)从小到大排序，则\(A_1=X_1+X_2,A_2=X_1+X_3\)。

而\(X_2+X_3\)的取值可能在\(X_3\sim X_n\)中，我们可以枚举\(X_2+X_3\)的取值，这样我们就可以求出\(X_1,X_2,X_3\)的值。

用multiset维护\(A\)，将\(X_1+X_2,X_1+X_3,X_2+X_3\)从集合中删去，剩下最小的一定是\(X_1+X_4\)。这样可以求出\(X_4\)。将\(X_4\)与\(X_1\sim X_3\)的和从集合中删去，剩下最小就是\(X_1+X_5\)。以此类推。

这样我们\(\mathcal O(n)\)枚举\(X_2+X_3\)的值，再\(\mathcal O(n^2\log n)\)不断推出\(X_4\sim X_n\)的值。时间复杂度\(\mathcal O(n^3\log n)\)。

源代码：

#include<set>
#include<cstdio>
#include<cctype>
#include<vector>
#include<algorithm>
inline int getint() {
	register char ch;
	while(!isdigit(ch=getchar()));
	register int x=ch^'0';
	while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
	return x;
}
const int N=301,M=(N-1)*(N-2)/2+1;
int a[M],cnt,ans[N][N],n;
std::multiset<int> set;
void solve(const int &k) {
	set.clear();
	for(register int i=3;i<=n*(n-1)/2;i++) {
		if(i!=k) set.insert(a[i]);
	}
	const int tmp=(a[1]+a[2]+a[k])>>1;
	ans[cnt][1]=tmp-a[k];
	if(ans[cnt][1]<=0) return;
	ans[cnt][2]=tmp-a[2];
	if(ans[cnt][2]<=0) return;
	ans[cnt][3]=tmp-a[1];
	if(ans[cnt][3]<=0) return;
	for(register int i=4;i<=n;i++) {
		const int val=*set.begin();
		ans[cnt][i]=val-ans[cnt][1];
		if(ans[cnt][i]<=0) return;
		for(register int j=1;j<i;j++) {
			const std::set<int>::iterator pos=set.find(ans[cnt][i]+ans[cnt][j]);
			if(pos==set.end()) return;
			set.erase(pos);
		}
	}
	cnt++;
}
int main() {
	n=getint();
	for(register int i=1;i<=n*(n-1)/2;i++) {
		a[i]=getint();
	}
	std::sort(&a[1],&a[n*(n-1)/2]+1);
	for(register int i=3;i<=n;i++) {
		if(i!=3&&a[i]==a[i-1]) continue;
		if((a[1]+a[2]+a[i])&1) continue;
		solve(i);
	}
	printf("%d\n",cnt);
	for(register int i=0;i<cnt;i++) {
		for(register int j=1;j<=n;j++) {
			printf("%d%c",ans[i][j]," \n"[j==n]);
		}
	}
	return 0;
}