`
xitonga
  • 浏览: 590208 次
文章分类
社区版块
存档分类
最新评论

FZU 1704 Turn off the light

阅读更多
高斯消元 套模板,
求出方程的自由变量个数,由于变量取值只有0,1;
故2^(var-k)种方法,var-k最大取100,用高精
#include <iostream>
#include <cstdlib>
#include <algorithm>
#include <cstdio>
#include <string>
#include <cstring>
using namespace std;
#define CL(a,b) memset(a,b,sizeof(a))
#define INF 100000000
const int M(200);
int a[M][M];
int free_x[M],num;
int x[M],y[M];
int gauss(int equ,int var)
{
	int i,j,k,col;
	int max_r;
//	num=0;
//	CL(free_x,0);
//	CL(x,0);
	for(k=0,col=0;k<equ&&col<var;k++,col++)
	{
		max_r=k;
		for(i=k+1;i<equ;i++)
			if(a[i][col]>a[max_r][col])
				max_r=i;
		if(a[max_r][col]==0)
		{
			k--;
			continue;
		}
		if(k!=max_r)
			for(i=0;i<var+1;i++)
				swap(a[k][i],a[max_r][i]);
		for(i=k+1;i<equ;i++)
			if(a[i][col]!=0)
				for(j=col;j<var+1;j++)
					a[i][j]=a[i][j]^a[k][j];
	}
	return var-k;
}
int gj[M];
int gaojin(int n)
{
	int i,j,k,jw;
	gj[0]=1;
	j=1;
	for(i=0;i<n;i++)
	{
		jw=0;
		for(k=0;k<j;k++)
		{
			gj[k]<<=1;
			gj[k]+=jw;
			jw=0;
			if(gj[k]>=10)
			{
				jw=gj[k]/10;
				gj[k]%=10;
			}
		}
		if(jw>0)
			gj[j++]=jw;
	}
	return j;
}
int main()
{
	int t,n,i,c,d,m,j;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%d",&n,&m);
		CL(a,0);
		for(i=0;i<n;i++)
			scanf("%d",&a[i][m]);
		for(i=0;i<m;i++)
		{
			scanf("%d",&d);
			for(j=0;j<d;j++)
			{
				scanf("%d",&c);
				a[c-1][i]=1;
			}
		}
	/*	for(i=0;i<n;i++)
		{
			for(j=0;j<n+1;j++)
				printf("%d ",a[i][j]);
			printf("\n");
		}*/
		int k=gauss(n,m);
		CL(gj,0);
		k=gaojin(k);
		for(i=k-1;i>=0;i--)
			printf("%d",gj[i]);
		printf("\n");
	}
}
分享到:
评论

相关推荐

Global site tag (gtag.js) - Google Analytics