BZOJ 5299: [Cqoi2018]解锁屏幕

状压DP

#include<cstdio>
using namespace std;
const int mod=1e8+7;
int F[1000005][25],dis[25][25],lim[1000005];
struct node{
	int x,y;
}E[25];
int main(){
	int n;
	scanf("%d",&n);
	for (int i=0; i<n; i++) scanf("%d%d",&E[i].x,&E[i].y);
	for (int S=0; S<n; S++)
		for (int T=S+1; T<n; T++)
			for (int k=0; k<n; k++){
				if (k==S || k==T) continue;
				int x=E[S].x-E[k].x,y=E[S].y-E[k].y;
				int X=E[T].x-E[k].x,Y=E[T].y-E[k].y;
				if (X*y==x*Y) {
					if (E[k].x>E[S].x && E[k].x>E[T].x) continue;
					if (E[k].x<E[S].x && E[k].x<E[T].x) continue;
					if (E[k].y>E[S].y && E[k].y>E[T].y) continue;
					if (E[k].y<E[S].y && E[k].y<E[T].y) continue;
					dis[S][T]|=(1<<k);
					dis[T][S]|=(1<<k);
				}
			}
	for (int i=0; i<n; i++) F[1<<i][i]=1;
	for (int i=1; i<(1<<n); i++) lim[i]=lim[i>>1]+(i&1);
	for (int i=0; i<(1<<n); i++)
		for (int S=0; S<n; S++)
			if (F[i][S]){
				for (int T=0; T<n; T++)
					if (!(i&(1<<T)) && (i&dis[S][T])==dis[S][T]) (F[i|(1<<T)][T]+=F[i][S])%=mod;
			}
	int ans=0;
	for (int i=0; i<(1<<n); i++)
		for (int S=0; S<n; S++)
			if (lim[i]>=4) (ans+=F[i][S])%=mod;
	printf("%d\n",ans);
	return 0;
}