SPOJ_839_OPTM

原创 2012年03月30日 13:03:04

最小割,该题因为要求XOR后的最小值,我们可以通过二进制的知识知道,XOR后的结果各个位置上是互不影响的,XOR只有在两数不同的时候才会得到1.

于是,我们在建图的时候可以将原图G中每个点属于V拆成31位,每位成为一个点,属于同一位的所有点构成一副新图G‘.

假如新图中任意两位所属的原图中的点相连,那么将新图中这两点连一条容量为1的边.

另外如果该位为已知位,假如它是1,那么将它与原点连一条容量为INF的边,如果它是0,那么将它与汇点连一条容量为INF的边。

对于未知位,将它与原点和汇点分别连一条容量为1的边.

最后做一遍最大流,然后从原点出发做一遍DFS,找到与原点相连的点,找到这些点中属于未知位的点,将他们赋值为1,并更新到原图节点。

之所以满足最小割模型是因为,我们要保证每一位要不是1(与原点相连),要不就是0(与汇点相连),又要保证他们XOR出来的值最小。因为他们XOR出来的值=(Mincut - |v|)*2^x(|v|代表未知位的个数,x代表第几位),从这个公式可以知道,只要保证了最小割,那么它XOR出来的值必然最小,所以可以用最小割求解。

/*
 *NAME: OPTM
 *LANG: C++
 */
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#define INF 0x7FFFFFFF
using namespace std;

int tt,n,m,s,t,ans,edgeNum;
struct node{
    int x,y;
}link_fuck[100001];//表示大点边的相连情况
int p[5001],pp[5001],ppp[5001];//p表示该节点数字,pp表示该节点剩余数字
struct edge_node{
    int v,c,next;
}e[100000];//边集数组
int v[5000],d[5000],c[5000];//v储存边集,d储存深度,c储存该层节点数
bool cover[5001];

void insert(int a,int b,int c){
    edgeNum++;
    e[edgeNum].v=b;e[edgeNum].c=c;e[edgeNum].next=v[a];v[a]=edgeNum;
    edgeNum++;
    e[edgeNum].v=a;e[edgeNum].c=0;e[edgeNum].next=v[b];v[b]=edgeNum;
}
void graphInit(){
    s=n+1;t=n+2;edgeNum=0;
    memset(cover,false,sizeof(cover));
    memset(v,0,sizeof(v));
    memset(d,0,sizeof(d));
    memset(c,0,sizeof(c));
    memset(e,0,sizeof(e));
}
void graphMake(){
    for (int i=1;i<=n;++i)
	if (p[i]==-1)
	    insert(s,i,1),insert(i,t,1);
	else{
	    int tmp = pp[i] % 2;pp[i]/=2;
	    tmp==1?insert(s,i,INF):insert(i,t,INF);
	}
    for (int i=1;i<=m;++i)
	insert(link_fuck[i].x,link_fuck[i].y,1),insert(link_fuck[i].y,link_fuck[i].x,1);
}
void update(int len,int pos){
    cover[pos]=true;
    int tmp=v[pos];
    while (tmp>0){
	if (!cover[e[tmp].v] && e[tmp].c>0){
	    if (p[e[tmp].v]==-1)
		ppp[e[tmp].v]+=pow(2,len);
	    update(len,e[tmp].v);
	}
	tmp=e[tmp].next;
    }
}
int opp(int x){
    return (x%2==1)?x+1:x-1;
}
int maxflow(int x,int lim){
    int ret,min_num,i,j,max_flow;
    if (x==t) return lim;
    max_flow=0;min_num=t;
    j=v[x];
    while (j!=0){
	if (e[j].c>0){
	    if (d[e[j].v]+1==d[x]){
		ret=min(lim,e[j].c);
		ret=maxflow(e[j].v,ret);
		e[j].c-=ret;e[opp(j)].c+=ret;
		lim-=ret;max_flow+=ret;
		if (d[s]>=t) return max_flow;
	    }
	    min_num=min(min_num,d[e[j].v]+1);
	}
	j=e[j].next;
    }
    if (d[x]==min_num) return max_flow;
    --c[d[x]];
    if (c[d[x]]==0) d[s]=t;
    d[x]=min_num;++c[d[x]];
    return max_flow;
}
void init(){
    scanf("%d %d",&n,&m);memset(link_fuck,0,sizeof(link_fuck));
    for (int i=1;i<=m;++i) scanf("%d %d",&link_fuck[i].x,&link_fuck[i].y);
    int know;scanf("%d",&know);memset(p,0,sizeof(p));
    for (int i=0,a,b;i<know;++i){scanf("%d %d",&a,&b);p[a]=b;}
    for (int i=1;i<=n;++i) {ppp[i]=pp[i]=p[i];if (!p[i]) p[i]-=1;}
}
void print(){
    for (int i=1;i<=n;++i) 
	printf("%d\n",ppp[i]);
}
int main(){
    scanf("%d",&tt);
    for (int k=1;k<=tt;++k){
	init();
	//solve
	for (int i=0;i<=31;++i){
	    //makeGraph
	    graphInit();
	    graphMake();
	    //get minCut
	    c[0]=n+2;ans=0;
	    while (d[s]<t) ans+=maxflow(s,INF);
	    //update
	    update(i,s);
	}
	print();
    }
    return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

SPOJ 839 Optimal Marks 最小割 经典 按位建图

胡伯涛论文中的一题,经典建模,由于二进制每一位异或不会相互影响,所以我们把问题转换模型,按位处理。 即已知一些点的标号0/1(还有些可以自己任意改),和一些边,边权定义为两端点标号的异或,要求边权和...

查找字符串中第一个出现一次的字符

#include using namespace std; // look up the first character in the string. // using hash table t...

OPTM-Optimal Marks-SPOJ839最小割

You are given an undirected graph G(V, E). Each vertex has a mark which is an integer from the range...

spoj839 Optimal Marks (最小割经典模型 位拆分)

转载自:http://blog.csdn.net/baidu_23081367/article/details/52403163 题目大意: 给定一张图,每个点有一个数字,有的点没有数字。...

zoj 2587 Unique Attack(判断最小割的唯一性)

Unique Attack Time Limit: 5 Seconds      Memory Limit: 32768 KB N supercomputers in the Unit...

SPOJ 839 Optimal Marks 最小割模型的转化(按位求最大流)

此题Amber的论文上还是有讲,建图的方法就不再赘述 题意描述:一个无向图,一些顶点权值已知而一些顶点权值未知,其中图中边的权值为其关联的两个顶点的异或值,现在让你在未知权值的顶点上填上权值后使得要...

spoj839(Optimal Marks)

省赛比完了,进了工作室学做项目,闲来无聊看看网络流建模,这题想半天没想出硬是想不出和网络流有关系,看了07年胡伯涛的论文才有的思路,觉得这题很经典,所以就写了一下。   第一次用VS写的题,小小的纪...

【bzoj2400】Spoj 839 Optimal Marks 二进制+最小割

Description定义无向图中的一条边的值为:这条边连接的两个点的值的异或值。 定义一个无向图的值为:这个无向图所有边的值的和。 给你一个有n个结点m条边的无向图。其中的一些点的值是给定的,而...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)