洛谷P3376 网络最大流

不是讲网络流,就是存个板子
另外我的Dinic跑得比EK慢一倍可还行(

附两份比较好的教程,均来自洛谷日报
E K \sf \color{blue}EK EK
D i n i c \sf \color{blue}Dinic Dinic


下面是代码

EK

#include<cstdio>
#include<cstring> 
const int N=10010; 
int n,m,l,r,t,p1,p2,x,y,c,mn,ans,v[N],last[N],Q[N];
struct node { int v,e; }pre[N];
struct Edge { int y,pre,v; }E[200010];
inline int min(int x,int y) { return (x<y?x:y); }

int bfs() {
	memset(v,0,sizeof v); memset(pre,0,sizeof pre);
	v[p1]=1; Q[l=r=1]=p1;
	while (l<=r) {
		t=Q[l]; ++l;
		for (int i=last[t]; i; i=E[i].pre)
			if (!v[y=E[i].y]&&E[i].v) {
				pre[y]={t,i}; if (y==p2) return 1;
				v[y]=1; Q[++r]=y;
			}
	}
	return 0;
}
int main() {
	scanf("%d%d%d%d",&n,&m,&p1,&p2);
	for (int i=1; i<=m; ++i) {
		scanf("%d%d%d",&x,&y,&c); 
		E[i<<1]={y,last[x],c}; last[x]=(i<<1);
		E[i<<1|1]={x,last[y],0}; last[y]=(i<<1|1);
	}
	while (bfs()) {
		mn=19260817;
		for (int i=p2; i!=p1; i=pre[i].v)
			mn=min(mn,E[pre[i].e].v);
		for (int i=p2; i!=p1; i=pre[i].v)
			E[pre[i].e].v-=mn,E[pre[i].e^1].v+=mn;
		ans+=mn;
	}
	printf("%d\n",ans);
	return 0;
}

Dinic

#include<cstdio>
#include<cstring>
const int N=10010,inf=0x3f3f3f3f;
int n,m,t1,t2,s,ans,l,r,x,y,c,last[N],v[N],dep[N],Q[N];
struct node { int y,pre,v; }E[200010];
inline int min(int x,int y) { return (x<y?x:y); } 

int bfs() {
	memset(dep,0x3f,sizeof dep); int t,y;
	for (Q[l=r=1]=t1,dep[t1]=0; l<=r; ++l) {
		v[t=Q[l]]=0;
		for (int i=last[t]; i; i=E[i].pre)
			if (E[i].v&&dep[y=E[i].y]>dep[t]+1)
				dep[y]=dep[t]+1,v[y]||(Q[++r]=y,v[y]=1);
	}
	return (dep[t2]!=inf); 
}
int dfs(int x,int s) {
	if (x==t2) return s; int s2,y;
	for (int i=last[x]; i; i=E[i].pre)
		if (E[i].v&&dep[y=E[i].y]==dep[x]+1)
			if (s2=dfs(y,min(E[i].v,s)))
				return E[i].v-=s2,E[i^1].v+=s2,s2;
	return 0;
}
int main() {
	scanf("%d%d%d%d",&n,&m,&t1,&t2);
	for (int i=1; i<=m; ++i) {
		scanf("%d%d%d",&x,&y,&c);
		E[i<<1]={y,last[x],c}; last[x]=(i<<1);
		E[i<<1|1]={x,last[y],0}; last[y]=(i<<1|1);
	}
	while (bfs()) while (s=dfs(t1,inf)) ans+=s;
	printf("%d\n",ans);
	return 0;
}
洛谷 P1681 最大正方形II 是一个动态规划问题,要求给定一个由 '0' 和 '1' 组成的矩阵,找出其中最大的正方形,并输出其边长。 以下是一个 C++ 编写的解答示例: ```cpp #include <iostream> #include <vector> #include <algorithm> using namespace std; int maximalSquare(vector<vector<char>>& matrix) { int rows = matrix.size(); if (rows == 0) return 0; int cols = matrix[0].size(); vector<vector<int>> dp(rows + 1, vector<int>(cols + 1, 0)); int maxSide = 0; for (int i = 1; i <= rows; i++) { for (int j = 1; j <= cols; j++) { if (matrix[i-1][j-1] == '1') { dp[i][j] = min(min(dp[i-1][j], dp[i][j-1]), dp[i-1][j-1]) + 1; maxSide = max(maxSide, dp[i][j]); } } } return maxSide * maxSide; } int main() { int n, m; cin >> n >> m; vector<vector<char>> matrix(n, vector<char>(m)); for (int i = 0; i < n; i++) { for (int j = 0; j < m; j++) { cin >> matrix[i][j]; } } cout << maximalSquare(matrix) << endl; return 0; } ``` 在上述代码中,我们首先定义了一个名为 `maximalSquare` 的函数,该函数接受一个二维字符矩阵 `matrix` 作为参数,返回最大正方形的边长。 在 `main` 函数中,我们首先从标准输入读取矩阵的行数和列数,并创建一个大小为 `n x m` 的二维字符矩阵。然后,我们按行读取矩阵的数据,并调用 `maximalSquare` 函数进行求解。最后,输出最大正方形的边长。 在动态规划的解法中,我们使用一个二维数组 `dp` 来记录以当前位置为右下角的最大正方形的边长。遍历矩阵中的每个元素,如果当前元素为 '1',则根据其左方、上方和左上方的最大正方形边长计算出当前位置的最大正方形边长,并更新 `dp` 数组和最大边长变量。 请注意,以上代码仅为示例,可能需要根据具体题目要求进行适当修改。同时,为了简化示例,未进行输入验证,请确保输入的矩阵符合题目要求。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值