USACO wormhole

  新题啊，网上代码好像都挺长的，自己写了个70行的，加上注释应该会好读一些吧，思路应该都没啥问题，枚举所有可能，然后判断就OK了，代码如下:

/*
ID: jinusac1
PROG: wormhole
LANG: C++
*/
#include <iostream>
#include <cstdio>
using namespace std;
int n,ans=0,e[20],g[20]={0};//g数组表示每个点延x轴正方向走可以遇到的最近的点，没有为0
long long x[20],y[20];
bool is[20]={false};//枚举点对时用到的判重
bool judge(int pos)
{
	int pt=pos;
	bool yes=false;
	while(true){
		pt=e[pt];
		//if(pt==pos) {yes=true;break;}这个语句得注意一下，此时是出来的，不能视为循环
		pt=g[pt];
		if(!pt) break;
		if(pt==pos) {yes=true;break;}
	}
	return yes;
}
void dfs(int total)
{
	if(total==(n>>1)){
		bool mark=false;
		for(int i=1;i<=n;i++) if(judge(i)){
			mark=true;
			break;
		}
		if(mark) ans++;
		return ;
	}
	int a;
	for(int i=1;i<=n;i++) if(!is[i]){
		is[i]=true;a=i;
		break;
	}
	for(int i=1;i<=n;i++){
		if(!is[i]){
			is[i]=true;
			e[i]=a;e[a]=i;
			dfs(total+1);
			is[i]=false;
		}
	}
	is[a]=false;
}
int main()
{
	freopen("wormhole.in","r",stdin);
	freopen("wormhole.out","w",stdout);
	cin>>n;
	for(int i=1;i<=n;i++) cin>>x[i]>>y[i];
	for(int i=1;i<=n;i++){
		int pos=0;
		for(int j=1;j<=n;j++){
			if(y[i]==y[j]&&x[j]>x[i]){
				if(pos==0||x[j]<x[pos])
					pos=j;
			}
		}
		g[i]=pos;
	}//构造g数组
	dfs(0);
	cout<<ans<<endl;
	return 0;
}