bzoj 1770/Luogu 2962 [Usaco2009 Nov]lights 燈

难得看到bzoj上有繁体字这么带感的东西。。。（语文老师说过 繁体字才有文化啊）

反正我看着也没什么 毕竟时不时在TVB或者在HK 什么的也能看到

-------------------------------------------可爱的分割线-----------------------------------------------------

这道题 应该是一道高斯消元，然而我不会

就是解异或方程组，用dfs跑一下自由元就好。。。

其实有另一个做法：折半搜索 orz HZWER 真的666啊

分开两半来搜 因为每个开关 你只有开或者不开两种情况 时间复杂度O(2^N)

把N/2之后来搜 然后用map记录状态 再合并就好啦

为什么他们说折半搜索很快。。412ms不快啊。难道是我偷懒用了map？....

#include<map>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long LL;
const int inf=1e9;
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
map<LL,int>f;
LL p[40],b[40],ed;
bool half;
int n,m,cnt,mn;
void dfs(int x,LL now,int u)
{
	if(x==cnt+1)
	{
		if(now==ed)mn=min(mn,u);
		if(!half)
		{
			int c=f[now];
			if(!c || c>u)f[now]=u;
		}
		else
		{
			int c=f[ed-now];
			if(c)mn=min(mn,c+u);
		}
		return;
	}
	dfs(x+1,now,u);
	dfs(x+1,now^p[x],u+1);
}
int main()
{
	int i,x,y; n=read(),m=read();
	p[1]=b[1]=1; for(i=2;i<=n+1;i++) p[i]=b[i]=b[i-1]<<1;
	ed=b[n+1]-1;
	while(m--)
	{
		x=read(),y=read();
		p[x]+=b[y],p[y]+=b[x];
	}
	mn=inf,cnt=n/2;
	dfs(1,0,0); half=1;
	cnt=n; dfs(n/2+1,0,0);
	printf("%d\n",mn);
	return 0;
}

高斯消元：

#include<cstdio>
#include<bitset>
#include<cstring>
#include<algorithm>
using namespace std;
const int inf=1e9,N=37;
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
int n,m,id[N]; bitset<N>mp[N];
void gauss()
{
	int k=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=k;j<=n;j++)
		  if(mp[j][i]) {if(j!=k)swap(mp[j],mp[k]); break;}
		if(mp[k][i])
		{
			for(int j=k+1;j<=n;j++)
			  if(mp[j][i])mp[j]^=mp[k];
			id[i]=k++;
		}
	}
}
int mn,p[N];
void dfs(int x,int c)
{
	if(c>=mn)return;
	if(x<1){mn=c; return;}
	if(!id[x]) dfs(x-1,c+(p[x]=1)),dfs(x-1,c+(p[x]=0));
	else
	{
		int k=id[x]; p[x]=mp[k][n+1];
		for(int i=x+1;i<=N;i++)if(mp[k][i])p[x]^=p[i];
		dfs(x-1,c+p[x]);
	}
}
int main()
{
	int i,x,y; n=read(),m=read();
	for(i=1;i<=n;i++)mp[i][i]=mp[i][n+1]=1;
	while(m--)
	{
		x=read(),y=read();
		mp[x][y]=mp[y][x]=1;
	}
	gauss();
	mn=inf; dfs(n,0);
	printf("%d\n",mn);
	return 0;
}