hihoCoder 最大权闭合子图

在一个图中,我们选取一些点构成集合,记为V,且集合中的出边(即集合中的点的向外连出的弧),所指向的终点(弧头)也在V中,则我们称V为闭合图。最大权闭合图即在所有闭合图中,集合中的点的权值之和最大的V,我们称V为最大权闭合图。

(可抽象为最大权闭合图求法的题一般是含有一些点具有一定的利益收入,另一些点具有一定的利益支出,即每个点都有一个权值。拥有某个具有收入效益的点前提是必须得拥有一些具有一定利益支出的点。求拥有哪些点能够获得最大的利益。我们可以对收入的点指向支出的点进行建边,则构建的图的每个闭合子图就是一种方式,求解其中最大权闭合子图)

解决该类问题的步骤是:
一。建立源点s和汇点t,然后将源点s与所有权值为正的点相连,容量为权值;将所有权值为负的点与汇点t相连,容量为权值的绝对值;权值为0的点不做处理;同时将原来的边容量设置为无穷大。
二。求解该图的最大流即最小割的值。
有一个结论:最大权闭合子图的权值等于所有正权点之和减去最小割。
解释:
最小割一定是简单割。
简单割一定和一个一个闭合子图对应。(不明白可以尝试在图上画一个割)
每个简单割的容量C(S, T) = T中正权点之和 + S中的负权点绝对值之和
一个闭合子图的权值W = S中的正权点之和 - S中负权点绝对值之和。
所以推出 W + C(S, T) = 所有正权点之和。
进而:闭合子图的权值 = 所有正权点之和 - C(S, T)。
由于所有正权点的权值之和是一个定值,那么割的容量越小,W也就越大。因此当C(S,T)取最小割时,W也就达到了最大权。

(不懂的看)

下面贴上一份Dinic代码和一份SAP代码。


Code1:

#include <algorithm>
#include <string.h>
#include <cstdio>
#include <queue>
using namespace std;
const int inf = 0x3f3f3f3f;
const int maxn = 405;
const int maxm = maxn*maxn;
const int BAS = 200;
struct node
{
	int v, w, next;
} edge[maxm];
int no, head[maxn];
int n, m, all, S, T;
int val[maxn];
queue<int> q;
int rec[maxn], pre[maxn], block[maxn], dis[maxn];
inline void init()
{
	no = 0; all = 0;
	memset(head, -1, sizeof head);
}
inline void add(int u, int v, int w)
{
	edge[no].v = v; edge[no].w = w;
	edge[no].next = head[u]; head[u] = no++;
	
	edge[no].v = u, edge[no].w = 0;
	edge[no].next = head[v]; head[v] = no++;
}
void mapping()
{
	int a, b, x;
	S = BAS*2+1, T = BAS*2+2;
	scanf("%d %d", &n, &m);
	for(int i = 1; i <= m; ++i)
	{
		scanf("%d", &x);
		add(BAS+i, T, x);
	}
	for(int i = 1; i <= n; ++i)
	{
		scanf("%d %d", &a, &b);
		add(S, i, a); all += a;
		for(int j = 1; j <= b; ++j)
		{
			scanf("%d", &x);
			add(i, BAS+x, inf);
		}
	}
}
void reset(int S, int T)
{ 
	memset(dis, 0x3f, sizeof dis);
	memset(block, 0, sizeof block);
	q.push(S); dis[S] = 0;
	while(!q.empty())
	{
		int top = q.front(); q.pop();
		for(int k = head[top]; k != -1; k = edge[k].next)
		if(dis[edge[k].v] == inf && edge[k].w)
			dis[edge[k].v] = dis[top]+1, q.push(edge[k].v);
	}
}
int dinic(int S, int T)
{
	int ans = 0, flow = inf, top = S;
	pre[S] = S;
	reset(S, T);
	while(dis[T] != inf)
	{
		int k, tmp;
		for(k = head[top]; k != -1; k = edge[k].next)
		{
			if(edge[k].w && dis[edge[k].v]==dis[top]+1 && 
			!block[edge[k].v]) break;
		}
		if(k != -1)
		{
			tmp = edge[k].v;
			flow = min(flow, edge[k].w);
			pre[tmp] = top, rec[tmp] = k;
			top = tmp;
			if(top == T)
			{
				ans += flow;
				for(; top != S; top = pre[top])
				{
					edge[rec[top]].w -= flow;
					edge[rec[top]^1].w += flow;
					if(!edge[rec[top]].w) tmp = top;
				}
				top = pre[tmp], flow = inf;
				for(; top != S; top = pre[top])
				flow = min(flow, edge[rec[top]].w);
				top = pre[tmp];
			}
		}
		else
		{
			block[top] = 1;
			top = pre[top];
			if(block[S]) reset(S, T);
		}
	}
	return ans;
}
int main()
{
	init();
	mapping();
	printf("%d\n", all-dinic(S, T));
	return 0;
}

Code2:

#include <algorithm>
#include <string.h>
#include <cstdio>
#include <queue>
using namespace std;
const int inf = 0x3f3f3f3f;
const int maxn = 410;
const int maxm = maxn*maxn;
const int BAS = 200;
struct node
{
	int v, w, next;
} edge[maxm];
int no, head[maxn];
int N, n, m, all, S, T;
int val[maxn];
queue<int> q;
int rec[maxn], pre[maxn], now[maxn], dis[maxn], gap[maxn];
inline void init()
{
	no = 0; all = 0;
	memset(head, -1, sizeof head);
}
inline void add(int u, int v, int w)
{
	edge[no].v = v; edge[no].w = w;
	edge[no].next = head[u]; head[u] = no++;
	
	edge[no].v = u, edge[no].w = 0;
	edge[no].next = head[v]; head[v] = no++;
}
void mapping()
{
	int a, b, x;
	S = BAS*2+1, T = BAS*2+2;
	scanf("%d %d", &n, &m);
	N = n+m+2;
	for(int i = 1; i <= m; ++i)
	{
		scanf("%d", &x);
		add(BAS+i, T, x);
	}
	for(int i = 1; i <= n; ++i)
	{
		scanf("%d %d", &a, &b);
		add(S, i, a); all += a;
		for(int j = 1; j <= b; ++j)
		{
			scanf("%d", &x);
			add(i, BAS+x, inf);
		}
	}
}
void PRE_SAP(int S, int T)
{
	memset(gap, 0, sizeof gap);
	memset(dis, 0, sizeof dis);
	//注意这里now数组要把所有用到的标号都存过来 
	//初始化一定要注意把所涉及的都覆盖到 
	for(int i = 0; i <= T; ++i) now[i] = head[i];
	q.push(T); dis[T] = 0;
	while(!q.empty())
	{
		int top = q.front(); q.pop();
		++gap[dis[top]];
		for(int k = head[top]; k != -1; k = edge[k].next)
		if(dis[edge[k].v] == inf && edge[k^1].w)
			dis[edge[k].v] = dis[top]+1, q.push(edge[k].v);
	}
}
int SAP(int S, int T)
{
	int k, ans = 0, flow = inf, top = S;
	pre[S] = S;
	PRE_SAP(S, T);
	while(dis[S] < N)
	{
		if(top == T)
		{
			ans += flow;
			for(; top != S; top = pre[top])
			{
				edge[rec[top]].w -= flow;
				edge[rec[top]^1].w += flow;
			}
			flow = inf;
		}
		for(k = now[top]; k != -1; k = edge[k].next)
		{
			if(edge[k].w && dis[top] == dis[edge[k].v]+1)
			{
				flow = min(flow, edge[k].w);
				pre[edge[k].v] = top, rec[edge[k].v] = k;
				now[top] = k; top = edge[k].v;
				break;
			}
		}
		if(k == -1)
		{
			int mins = N;
			if(--gap[dis[top]] == 0) break;
			k = now[top] = head[top];
			for(; k != -1; k = edge[k].next)
			if(edge[k].w && mins > dis[edge[k].v])
				mins = dis[edge[k].v];
			++gap[dis[top] = mins+1];
			top = pre[top];
		}
	}
	return ans;
}
int main()
{
	init();
	mapping();
	printf("%d\n", all-SAP(S, T));
	return 0;
}

继续加油~
发布了288 篇原创文章 · 获赞 413 · 访问量 18万+
展开阅读全文

用c++ 求智能手机的滑动解锁方式一共多少种?3*3格

03-17

下面是我写的C++,思想是枚举方向,如果下一个方向被标记,则按原来方向再搜下一个点,否则跳到下一个点并标记 按说 枚举每个点为起点的时候 1,3,7,9跑出来的ans是相同的; 2,4,6,8跑出来的ans也是相同的,为什么 我这代码跑出来2,4,6,8,不是相同的,并且最后加起来的答案应该是30多万,我的答案相差很大,求大神解!!! ``` #include <iostream> #include <string.h> using namespace std; const int maxn=1e5+10; int vis[10],vis1[10],in[10],out[10],head[200],q,ans,a[5][5]; int dir[12][2]={1,0,0,1,-1,0,0,-1,-1,-1,-1,1,1,1,1,-1,-2,-1,2,1,-2,1,2,-1}; struct point { int to,next; }pt[200]; struct pp{ int x; int y; }sto[10]; void add(int u,int v) { pt[q].next=head[u]; pt[q].to=v; head[u]=q++; } void dfs(int po,int x,int y,int len) { //cout<<po<<"~~~~~~~~"<<endl; if (len>=4) ans++; int flag=0; /*for (int i=head[po];i!=-1;i=pt[i].next) { cout<<"%%%%%"<<endl; int to=pt[i].to; if (!vis[to]) { flag++; vis[to]=1; dfs(to,sto[to].x,sto[to].y,len+1); vis[to]=0; } }*/ if (!flag) { for (int i=0;i<12;i++) //枚举12个方向 { int xx=x+dir[i][0]; int yy=y+dir[i][1]; if (xx>=1&&xx<=3&&yy>=1&&yy<=3&&!vis[a[xx][yy]]) //下一步如果没被标记; { vis[a[xx][yy]]=1; dfs(a[xx][yy],xx,yy,len+1); vis[a[xx][yy]]=0; } else if (xx>=1&&xx<=3&&yy>=1&&yy<=3&&i<=7)// 下一步被标记,按照方向再走一步,看是否被标记 { int xxx=xx+dir[i][0]; int yyy=yy+dir[i][1]; if (xxx>=1&&xxx<=3&&yyy>=1&&yyy<=3&&!vis[a[xxx][yyy]]) { vis[a[xxx][yyy]]=1; dfs(a[xxx][yyy],xxx,yyy,len+1); vis[a[xxx][yyy]]=0; } } } } } int main() { int t,len,s=1,n,u,v; for (int i=1;i<=3;i++) for (int j=1;j<=3;j++) a[i][j]=s++; cin>>t; for (int i=1;i<=t;i++) { int k=1; ans=0; cin>>n; memset(vis,0,sizeof(vis)); memset(head,-1,sizeof(head)); q=0; for (int i=1;i<=n;i++) { cin>>u>>v; add(u,v); add(v,u); } for (int i=1;i<=3;i++) //枚举每个点当起点 for (int j=1;j<=3;j++) { ans=0; len=1; //sto[k].x=i; //sto[k].y=j; vis[k]=1; dfs(k++,i,j,len); cout<<ans<<"######"<<endl; vis[k-1]=0; for (int s=1;s<=9;s++) cout<<vis[s]<<" "; cout<<endl; } //cout<<ans<<endl; } return 0; } ``` 问答

hihocoder 1014的相关问题

09-11

大家好,想请教大家关于hihocoder 1014的相关问题 问题描述: /*时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi和小Ho是一对好朋友,出生在信息化社会的他们对编程产生了莫大的兴趣,他们约定好互相帮助,在编程的学习道路上一同前进。 这一天,他们遇到了一本词典,于是小Hi就向小Ho提出了那个经典的问题:“小Ho,你能不能对于每一个我给出的字符串,都在这个词典里面找到以这个字符串开头的所有单词呢?” 身经百战的小Ho答道:“怎么会不能呢!你每给我一个字符串,我就依次遍历词典里的所有单词,检查你给我的字符串是不是这个单词的前缀不就是了?” 小Hi笑道:“你啊,还是太年轻了!~假设这本词典里有10万个单词,我询问你一万次,你得要算到哪年哪月去?” 小Ho低头算了一算,看着那一堆堆的0,顿时感觉自己这辈子都要花在上面了... 小Hi看着小Ho的囧样,也是继续笑道:“让我来提高一下你的知识水平吧~你知道树这样一种数据结构么?” 小Ho想了想,说道:“知道~它是一种基础的数据结构,就像这里说的一样!” 小Hi满意的点了点头,说道:“那你知道我怎么样用一棵树来表示整个词典么?” 小Ho摇摇头表示自己不清楚。 提示一:Trie树的建立 “你看,我们现在得到了这样一棵树,那么你看,如果我给你一个字符串ap,你要怎么找到所有以ap开头的单词呢?”小Hi又开始考校小Ho。 “唔...一个个遍历所有的单词?”小Ho还是不忘自己最开始提出来的算法。 “笨!这棵树难道就白构建了!”小Hi教训完小Ho,继续道:“看好了!” 提示二:如何使用Trie树 提示三:在建立Trie树时同时进行统计! “那么现在!赶紧去用代码实现吧!”小Hi如是说道 输入 输入的第一行为一个正整数n,表示词典的大小,其后n行,每一行一个单词(不保证是英文单词,也有可能是火星文单词哦),单词由不超过10个的小写英文字母组成,可能存在相同的单词,此时应将其视作不同的单词。接下来的一行为一个正整数m,表示小Hi询问的次数,其后m行,每一行一个字符串,该字符串由不超过10个的小写英文字母组成,表示小Hi的一个询问。 在20%的数据中n, m<=10,词典的字母表大小<=2. 在60%的数据中n, m<=1000,词典的字母表大小<=5. 在100%的数据中n, m<=100000,词典的字母表大小<=26. 本题按通过的数据量排名哦~ 输出 对于小Hi的每一个询问,输出一个整数Ans,表示词典中以小Hi给出的字符串为前缀的单词的个数。 样例输入 5 babaab babbbaaaa abba aaaaabaa babaababb 5 babb baabaaa bab bb bbabbaab 样例输出 1 0 3 0 0*/ 我的代码: #include <stdio.h> #include <stdlib.h> struct trie { int num; struct trie *data[26]; }*p; void trie1(struct trie *x) { int i; x->num = 0; for(i = 0; i < 26; i++) x->data[i] = NULL; } void insertTrie(char s[]) { int i; trie1(p); p->num++; for(i=0;s[i];i++) { int j = s[i] - 'a'; if(p->data[j] == NULL) { p->data[j] = (struct trie *)malloc(sizeof(struct trie)); trie1(p->data[j]); } p = p->data[j]; p->num++; } } int searchTrie(char s[]) { int i; //struct trie *q = root ; //trie1(p); for(i = 0; p&&s[i]; i++) { int j = s[i] - 'a'; p = p->data[j]; } if(!p) return 0; else return p->num; } int main() { char s[12]; int n, m; while(~scanf("%d", &n)) { //root = (struct Trie *)malloc(sizeof(struct trie)); while(n--) { scanf("%s", s); insertTrie(s); } scanf("%d", &m); while(m--) { scanf("%s", s); printf("%d\n", searchTrie(s)); } } return 0; } 不知道有什么问题啊,但是运行出错,希望各位大神帮我看看,谢谢! 问答

我的hihocoder这个最长回文串为什么报wrong answer?

01-29

import java.util.Scanner; public class Main { public static void main(String[] args) { Scanner scanner = new Scanner(System.in); int n = scanner.nextInt(); for(int i = 0;i<n;i++){ String str = scanner.next(); System.out.println(findMax(buildStr(str))/2); } } public static int findMax(String str){ int pArr[] = new int[str.length()]; int pr = 0; int index = 0; for(int i = 0;i<str.length();i++){ if(i==0){ pArr[i] = expand(i,str); pr=pArr[i]+index; }else{ if(i<pr){ if(2*index-i-pArr[2*index-i]>index-pArr[index]){ pArr[i] = 2*index-i-(index-pArr[index]); index = i; pr = index+pArr[i]; }else if(2*index-i-pArr[2*index-i]<index-pArr[index]){ pArr[i] = pArr[2*index-i]; }else{ pArr[i] = expandFromPr(i,pr,str); index = i; pr = index+pArr[i]; } }else{ pArr[i] = expand(i,str); index = i; pr = index+pArr[i]; } } } int max = 0; for(int k:pArr){ max = Math.max(k*2-1, max); } return max; } public static int expandFromPr(int i,int pr,String str){ int count = pr-i; int j = pr-i; while(i-j>=0&&i+j<str.length()&&str.charAt(i-j)==str.charAt(i+j)){ j++; count++; } return count; } public static int expand(int i,String str){ int count = 1; int j = 1; while(i-j>=0&&i+j<str.length()&&str.charAt(i-j)==str.charAt(i+j)){ j++; count++; } return count; } public static String buildStr(String str){ StringBuilder strb = new StringBuilder(); int index = 0; int length = str.length(); length=length*2+1; for(int i = 0;i<length;i++){ if(i%2==0) strb.append('#'); else strb.append(str.charAt(index++)); } return strb.toString(); } } 问答

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 编程工作室 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览