zoukankan      html  css  js  c++  java
  • bzoj 1770: [Usaco2009 Nov]lights 燈【高斯消元+dfs】

    参考:https://blog.csdn.net/qq_34564984/article/details/53843777
    可能背了假的板子……
    对于每个灯建立方程:与它相邻的灯的开关次数的异或和为1
    异或高斯消元,然后dfs,遇到自由元就分两种情况走,答案取max,加上最优性剪枝

    #include<iostream>
    #include<cstdio>
    #include<cmath>
    using namespace std;
    const int N=45;
    int n,m,a[N][N],ans[N],mn=1e9,tot;
    int read()
    {
    	int r=0,f=1;
    	char p=getchar();
    	while(p>'9'||p<'0')
    	{
    		if(p=='-')
    			f=-1;
    		p=getchar();
    	}
    	while(p>='0'&&p<='9')
    	{
    		r=r*10+p-48;
    		p=getchar();
    	}
    	return r*f;
    }
    void gaosi()
    {
    	// for(int i=1;i<=n;i++)
    	// {
    		// int now=i;
    		// for(int j=i+1;j<=n;j++)
    			// if(abs(a[now][i])<abs(a[j][i]))
    				// now=j;
    		// for(int j=i;j<=n+1;j++)
    			// swap(a[now][j],a[i][j]);
    		// for(int j=i;j<=n+1;j++)
    			// a[i][j]/=a[i][i];
    		// for(int j=i+1;j<=n;j++)
    		// {
    			// for(int k=i+1;k<=n+1;k++)
    				// a[j][k]^=(a[j][i]^a[i][k]);
    			// a[j][i]=0;
    		// }
    	// }
    	for(int i=1;i<=n;i++)  
        {  
            int j=i;  
            while(j<=n&&!a[j][i])
    			j++;  
            if(j>n)
    			continue;  
            if(i!=j)
    			for(int k=1;k<=n+1;k++)
    				swap(a[i][k],a[j][k]);  
            for(int j=1;j<=n;j++)  
                if(i!=j&&a[j][i])  
                    for(int k=1;k<=n+1;k++)  
                        a[j][k]^=a[i][k];  
        }  
    }
    void dfs(int u)
    {
    	if(tot>=mn)
    		return;
    	if(!u)
    	{
    		mn=tot;
    		return;
    	}
    	if(a[u][u])
    	{
    		int t=a[u][n+1];
    		for(int i=u+1;i<=n;i++)
    			if(a[u][i])
    				t^=ans[i];
    		ans[u]=t;
    		if(t)
    			tot++;
    		dfs(u-1);
    		if(t)
    			tot--;
    	}
    	else
    	{
    		tot++,ans[u]=1;
    		dfs(u-1);
    		tot--,ans[u]=0;
    		dfs(u-1);
    	}
    }
    int main()
    {
    	n=read(),m=read();
    	for(int i=1;i<=n;i++)
    		a[i][i]=a[i][n+1]=1;
    	for(int i=1;i<=m;i++)
    	{
    		int x=read(),y=read();
    		a[x][y]=a[y][x]=1;
    	}
    	gaosi();
    	dfs(n);
    	printf("%d
    ",mn);
    	return 0;
    }
    
  • 相关阅读:
    安装和配置SQL
    自定义中间件封装为模块
    留言板
    Spring Factories机制
    JVM参数配置
    全链路监控的方案概述与比较
    Lombok注解-@SneakyThrows
    @FeignClient注解 中属性 contextId使用
    @validate或@valid注解进行数据校验的解决方案
    maven的settings.xml标签详解
  • 原文地址:https://www.cnblogs.com/lokiii/p/8989722.html
Copyright © 2011-2022 走看看