Weekly8作业——差分约束/拓扑排序/Kosaraju

最新推荐文章于 2021-11-27 16:34:42 发布

明天可以有猫吗

最新推荐文章于 2021-11-27 16:34:42 发布

阅读量156

点赞数

分类专栏：程序设计思维与实践

本文链接：https://blog.csdn.net/weixin_43793842/article/details/105571861

版权

程序设计思维与实践专栏收录该内容

20 篇文章 0 订阅

订阅专栏

A - 区间选点 II

题目描述

给定一个数轴上的 n 个区间，要求在数轴上选取最少的点使得第 i 个区间 [ai, bi] 里至少有 ci 个点

Input

输入第一行一个整数 n 表示区间的个数，接下来的 n 行，每一行两个用空格隔开的整数 a，b 表示区间的左右端点。1 <= n <= 50000， 0 <= ai <= bi <= 50000 并且 1 <= ci <= bi - ai+1。

Output

输出一个整数表示最少选取的点的个数

Sample Input
5
3 7 3
8 10 3
6 8 1
1 3 1
10 11 1
Sample Output
6

思路

构造不等式组
- 记dis[i]表示数轴上[0,i]之间选点的个数
- 对于第i个区间[ai,bi]需要满足dis[bi]-dis[ai-1]≥ci
保证dis有意义：0≤dis[i]-dis[i-1]≤1
求该差分约束系统的最小解转化为≥不等式组用SPFA跑最长路，答案为dis[max{bi}]

总结

由于ai可能等于0，所以将ai++，bi++，避免出现ai-1<0的情况
求解差分约束系统，都可以通过不等式组的转换→图论中的单源最短路问题

实现代码

#include<iostream>
#include<cmath>
#include<cstdio>
#include<string.h>
#include<algorithm>
#include<queue>
const int N = 5e5+100;
const int M = 5e5+100;
const int inf =1e5;
using namespace std;
 
struct edge{
	int to,w,next;
}e[M];
int head[N],tot,dis[N],cnt[N];
queue<int>q;
int a,b,c,n;
int Max;

void add(int x,int y,int w)
{
	e[++tot].to=y,e[tot].next=head[x];
	e[tot].w=w,head[x]=tot;
}

void spfa(int s)
{
	memset(dis,inf,sizeof(dis));
	q.push(s);dis[s]=0;
	while(!q.empty())
	{
		int u=q.front();q.pop();cnt[u]=0;
		for(int i=head[u];i!=-1;i=e[i].next)
		{
			int v=e[i].to;
			if(dis[v]<dis[u]+e[i].w)
			{
				dis[v]=dis[u]+e[i].w;
				if(!cnt[v]) cnt[v]=1,q.push(v);
			}
		}
	}
}
 
int main()
{
	memset(head,-1,sizeof(head));
	cin>>n;
	for(int i=0;i<n;i++)
	{
		cin>>a>>b>>c;
		a++;b++;//避免出现负数
		add(a-1,b,c);
		Max=max(Max,b);
	}
	for(int i=1;i<=Max;i++)
	{
		add(i-1,i,0);
		add(i,i-1,-1);
	} 
	spfa(0);
	cout<<dis[Max];
	return 0;
}

B - 猫猫向前冲

题目描述

众所周知， TT 是一位重度爱猫人士，他有一只神奇的魔法猫。
有一天，TT 在 B 站上观看猫猫的比赛。一共有 N 只猫猫，编号依次为1，2，3，…，N进行比赛。比赛结束后，Up 主会为所有的猫猫从前到后依次排名并发放爱吃的小鱼干。不幸的是，此时 TT 的电子设备遭到了宇宙射线的降智打击，一下子都连不上网了，自然也看不到最后的颁奖典礼。
不幸中的万幸，TT 的魔法猫将每场比赛的结果都记录了下来，现在他想编程序确定字典序最小的名次序列，请你帮帮他。

Input

输入有若干组，每组中的第一行为二个数N（1<=N<=500），M；其中N表示猫猫的个数，M表示接着有M行的输入数据。接下来的M行数据中，每行也有两个整数P1，P2表示即编号为 P1 的猫猫赢了编号为 P2 的猫猫。

Output

给出一个符合要求的排名。输出时猫猫的编号之间有空格，最后一名后面没有空格!
其他说明：符合条件的排名可能不是唯一的，此时要求输出时编号小的队伍在前；输入数据保证是正确的，即输入数据确保一定能有一个符合要求的排名。

Sample Input
4 3
1 2
2 3
4 3
Sample Output
1 2 4 3

思路

使用优先队列替换队列，保证字典序最小
将所有入度为0的点插入队列中
不断地从队列中取出队首元素，将其储存在数组中，将其能够到达的点的入度-1，若删除边后点的入度为0，就把它插入队列中
当队列为空时，检查图中是否存在边，若存在，则该图一定有环路，否则，数组中的元素即为结果

总结

注意输出，最后一名后面没有空格
字典序最小→优先队列

实现代码

#include<bits/stdc++.h>
using namespace std;
 
struct edge{
	int to,w,next;
}e[506*506];
int head[506],tot,num[506];
priority_queue<int,vector<int>,greater<int>> q;
vector<int>res;
int m,n,p1,p2;

void add(int x,int y)
{
	e[++tot].to=y;
	e[tot].next=head[x];
	head[x]=tot;
}

void toposort()
{
	for(int i=1;i<=n;i++)
	{
		if(num[i]==0) q.push(i);
	}
	while(!q.empty())
	{
		int u=q.top();
		q.pop();
		res.push_back(u);
		for(int i=head[u];i!=-1;i=e[i].next)
		{
			int v=e[i].to;
			if(--num[v]==0)
			{
				q.push(v);
			}
		}
	}
}
 
int main()
{
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		for(int i=0;i<506;i++)
		{
			head[i]=-1;
			num[i]=0;
			tot=0;
			res.clear();
			while(!q.empty())
				q.pop();
		}
		for(int i=0;i<m;i++)
		{
			cin>>p1>>p2;
			add(p1,p2);
			num[p2]++;
		}
		toposort();
		for(int i=0;i<n;i++)
		{
			if(i!=n-1) cout<<res[i]<<" ";
			else cout<<res[i];
		}
		cout<<endl;
	}
	return 0;
}

C - 班长竞选

题目

大学班级选班长，N 个同学均可以发表意见若意见为 A B 则表示 A 认为 B 合适，意见具有传递性，即 A 认为 B 合适，B 认为 C 合适，则 A 也认为 C 合适勤劳的 TT 收集了M条意见，想要知道最高票数，并给出一份候选人名单，即所有得票最多的同学，你能帮帮他吗？

Input

本题有多组数据。第一行 T 表示数据组数。每组数据开始有两个整数 N 和 M (2 <= n <= 5000, 0 <m <= 30000)，接下来有 M 行包含两个整数 A 和 B(A != B) 表示 A 认为 B 合适。

Output

对于每组数据，第一行输出 “Case x: ”，x 表示数据的编号，从1开始，紧跟着是最高的票数。接下来一行输出得票最多的同学的编号，用空格隔开，不忽略行末空格！

Sample Input
2
4 3
3 2
2 0
2 1

3 3
1 0
2 1
0 2
Sample Output
Case 1: 2
0 1
Case 2: 2
0 1 2

思路

用Kosaraju求SCC
缩点：将每个SCC看成一个点，遍历反图中的边，将不再同一SCC的边加入，得到缩点之后的图
分析可得：对于属于第ｉ个SCC的点来说，答案分为两部分，令SCC[i]表示第i个SCC中点的个数
　 - 当前SCC中的点：ans=ans+SCC[i]-1
　 - 其他SCC中的点：sum（SCC[j]），其中j可以到达i
可得答案一定在出度为0的SCC中，因此只要将边反向，对每个出度为0的点进行dfs，计算能够到达的点的sum（SCC[j]），即为答案

总结

~~要用到的数组好多，记得一定要初始化~~
需要用到三张图：原图、反图、缩点后的图；三次DFS
票数最多的点一定出现在出度为0的SCC中，可以通过反证法证明
注意输出格式：行末无空格

实现代码

#include<iostream>
#include<vector>
#include<cstring>
using namespace std;
const int N=5006,M=30006;
int n,m,c[N],dfn[N],vis[N],dcnt,scnt,head1[N],head2[N],head3[N],tot1,tot2,tot3,scc[N],in_deg[N],num,mscc[N],res[N];
bool judge[N];
int a,b,Max;
vector<int> G1[N],G2[N];//G1-原图，G2-反图
struct edge {
	int u, v, next;
}e1[M],e2[M],e3[M];
//原图
void add(int u, int v,edge e[],int &tot,int head[]) 
{
	e[++tot].u = u;
	e[tot].v = v;
	e[tot].next = head[u];
	head[u] = tot;
}

void dfs1(int x) {
	vis[x] = 1;
	for (int i = head1[x]; i!=-1; i = e1[i].next)
	{
		int y = e1[i].v;
		if (!vis[y]) 
		{ 
			dfs1(y);
		}
	}
	dfn[++dcnt] = x;
}

void dfs2(int x) {
	c[x] = scnt;
	scc[scnt]++;
	for (int i = head2[x]; i!=-1; i = e2[i].next)
	{
		int y = e2[i].v;
		if (!c[y])
		{
			dfs2(y);
		}
	}
}

void dfs3(int x)
{
	num=num+scc[x],judge[x]=true;
	for(int i=head3[x];i!=-1;i=e3[i].next)
	{
		if(!judge[e3[i].v])
		{
			dfs3(e3[i].v);
		}
	}
}
 
void kosaraju() 
{
	dcnt = 0, scnt = 0;
	memset(c, 0, sizeof(c));
	memset(vis, 0, sizeof(vis));
	for (int i = 1; i <= n; i++)
	{
		if (!vis[i])
		{
			dfs1(i);
		}		
	}
	for (int i=n; i>=1; i--)
	{
		if (!c[dfn[i]])
		{
			++scnt;
			dfs2(dfn[i]);
		}
	}
}

int main()
{
	int t;
	cin>>t;
	for(int p=0;p<t;p++)
	{
		for(int i=0;i<N;i++)
		{
			head1[i]=-1,head2[i]=-1,head3[i]=-1,tot1=0,tot2=0,tot3=0;
			dcnt=0,scnt=0,num=0,Max=0;
			c[i]=0,dfn[i]=0,vis[i]=0,scc[i]=0,in_deg[i]=0,mscc[i]=0,res[i]=0;
			judge[i]=false;
		}
		cin>>n>>m;
		for(int i=0;i<m;i++)
		{
			cin>>a>>b;
			a++,b++;
			add(a,b,e1,tot1,head1);
			add(b,a,e2,tot2,head2);
		}
		kosaraju();
		for(int q=1;q<=n;q++)
		{
			for(int i=head2[q];i!=-1;i=e2[i].next)
			{
				if(c[q]==c[e2[i].v])
				{
					continue;
				}
				add(c[q],c[e2[i].v],e3,tot3,head3);
				in_deg[c[e2[i].v]]++;
			}
		}
		for(int i=1;i<=scnt;i++)
		{
			if(in_deg[i]==0)
			{
				num=0;
				for(int j=0;j<N;j++)
					judge[j]=false;
				dfs3(i);
				res[i]=res[i]+num-1;
			}
		}
		for(int i=1;i<=scnt;i++)
		{
			if(Max<res[i])
				Max=res[i];
		}
		for(int i=1;i<=scnt;i++)
		{
			if(res[i]==Max)
				mscc[i]=1;
		}
		cout<<"Case "<<p+1<<": "<<Max<<endl;
		bool it=false;
		for(int i=1;i<=n;i++)
		{
			if(mscc[c[i]]!=0)
			{
				if(it==false)
				{
					cout<<i-1;it=true;
				}
				else 
					cout<<" "<<i-1;
			}
		}
		cout<<endl;
	}
	return 0;
}

明天可以有猫吗

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
Weekly8作业——差分约束/拓扑排序/Kosaraju

A - 区间选点 II题目描述给定一个数轴上的 n 个区间，要求在数轴上选取最少的点使得第 i 个区间 [ai, bi] 里至少有 ci 个点Input输入第一行一个整数 n 表示区间的个数，接下来的 n 行，每一行两个用空格隔开的整数 a，b 表示区间的左右端点。1 <= n <= 50000， 0 <= ai <= bi <= 50000 并且 1 &lt...
复制链接

扫一扫