程序设计思维与实践 Week8 作业

Problem A - 区间选点 II

Description

给定一个数轴上的 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 Onput

6

Note

Idea

如果用差分约束系统来解决题目,那就需要先找到所有的约束。 第 i 个区间 [ai, bi] 里至少有 ci 个点 那也就是说,如果设点i到1需要放的点数为Sum[i],就有Sum[bi]-Sum[ai-1]>=ci.但是也要注意到 Sum[x]<=Sum[x+1]但是Sum[x]-Sum[x-1]<=1。就是说后面的点的S值肯定>=前一顶点,但最多也就多1。所以整理所有约束条件:

Sum[bi]   = Sum[ai-1]+ci
Sum[i]    = Sum[i-1]+0
Sum[i-1]  = Sum[i]-1 

把从0---MAX(bi)视为图中的各个顶点,边界0视为源点,那么每个点到前一顶点的有向边距离认为必须>=0,每个点到后一顶点的有向边距离必须>=-1.某个点bi到点ai-1之间的距离必须>=ci,可以先认为他们最开始的权就是这些限制。然后需要求在已知限制的情况下取得图中MAX(bi)点到源点的最长路径的权值和的最小解。===>可以使用SPFA求最长路、判断条件设置>=来处理。

Codes

重测数据

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
const int MAXN=200005;
const int maxn=50005;
const int inf=1e8;
int n, MAX;
struct Edge{
	int to,w,next;
}e[MAXN];
int head[maxn],sum[maxn],tot;
bool vis[maxn];
queue<int> q;
void init(){
	tot=0;MAX=-1;
	memset(head,-1,sizeof(head));
}

void addEdge(int u,int v,int w){ 
	e[tot].to=v;
	e[tot].w=w;
	e[tot].next=head[u];
	head[u]=tot;
	tot++;
}
void spfa(int s) {
	for(int i = 0; i <= MAX; ++i) 
		sum[i] = -inf;
	memset(vis,false,sizeof(vis));
	q.push(s); sum[s] = 0; vis[s] = true;  
	while(!q.empty()) {
		int u = q.front(); q.pop(); 
		vis[u] = false; 
		 for(int i = head[u]; i!=-1; i = e[i].next) {
			int v = e[i].to;  
			if(sum[v] < sum[u] + e[i].w) {
				sum[v] = sum[u] + e[i].w;
				if(!vis[v]) {
					q.push(v);
					vis[v] = true;
				}
			}
		 } 
	}
	cout<<sum[MAX]; 
}
int main(){
    scanf("%d",&n);
    int u,v,w;
    init();
    for(int i=0;i<n;++i){
    	scanf("%d%d%d",&u,&v,&w); 
		MAX=v>MAX?v:MAX;
    	addEdge(u,v+1,w); 
    }  
    ++MAX;
    for(int i=0;i<=MAX;++i) {
    	addEdge(i,i+1,0);
    	addEdge(i+1,i,-1);
	}
    spfa(0);
	return 0;
}

Problem B - 猫猫向前冲

Description

众所周知, 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 Onput

1 2 4 3

Idea

拓扑排序问题,另外需要注意要求最小字典序。拓扑排序可以用Kahn算法。

将所有入度=0的点组成集合。
• 每次从 S 里面取出一个顶点 u (可以随便取)放入 L , 然后遍历顶点 u 的 所有边 (u, v) , 并删除之,并判断如果该边的
另一个顶点 v,如果在移除 这一条边后入度为 0 , 那么就将这个顶点放入集合 S 中。不断地重复取 出顶点然后重复这个过程
• 最后当集合为空后,就检查图中是否存在任何边。如果有,那么这个图 一定有环路,否者返回 L , L 中顺序就是拓扑排序的结果

由于要确保最小字典序,这里我选择了优先级队列(升序)。每次把最小的入度=0的顶点从队列取出,保证了在排序合理的情况下能够让编号偏小的队伍在前。

Codes

在这里插入图片描述

#include<iostream>
#include<queue>
#include<vector>
#include<cstring>
using namespace std;
const int MAXN=505;
const int MAXM=5000;
int N,M;
struct Edge{
	int to,next;
}e[MAXM];
int head[MAXN],InDeg[MAXN],tot;
bool vis[MAXM];
struct cmp{
   bool operator()(int r, int l){ return r>l; } 
};
priority_queue<int,vector<int>,cmp > q;
vector<int> ans;
void init(){
	tot=0;
	memset(head,-1,sizeof(head));
	memset(InDeg,0,sizeof(InDeg));
	ans.clear();
}
void addEdge(int p1,int p2){ 
	e[tot].to=p2; 
	e[tot].next=head[p1];
	head[p1]=tot;
	++InDeg[p2];
	tot++;
} 
void Kahn(){
	for(int i=1;i<=N;++i){
		if(!InDeg[i]) q.push(i);
	}
	memset(vis,true,sizeof(vis));
	while(q.size()){
		int p1=q.top(); q.pop();
		ans.push_back(p1);
		for(int i = head[p1]; i!=-1; i = e[i].next){
		  int p2 = e[i].to;
		  if(vis[i]){
		  	 if(--InDeg[p2]==0) q.push(p2);
		  	 vis[i]=false;//边被删除 
		  }	
		}
	}
	cout<<ans[0]; 
	for(int i=1;i<N;++i) cout<<" "<<ans[i];
	cout<<endl; 
}
int main(){ 
    int p1,p2;
	while(~scanf("%d%d",&N,&M)) {
      init();
	  for(int i=0;i<M;++i){
		scanf("%d%d",&p1,&p2);
		addEdge(p1,p2);
	  }
	  Kahn();
	} 
    return 0;
}

Problem C - 班长竞选

Description

大学班级选班长,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 Onput

Case 1: 2
0 1
Case 2: 2
0 1 2

Note

PE多次,这次题目说的不要忽略行末空格我以为是行末空格也要带上,然而并不是。另外刚开始写了一遍然后运行发现一个样例都没有得到正确结果,之后分析是把后序,逆后序,生成的SCC编号之间的对应关系搞错了。

Idea

整体思路:如果把每个人看作一个点,那么那些互相可以通过意见传递发现两个人意见相互连通,那他们受到的支持票数是一样的。这些点就相当于一个有向图中的强连通分量。我们就可以先把点缩小至对应的强连通分量,因为一个分量中的点的票数都是一样的。找到强连通分量后,在原图遍历所有出度=0的强连通分量SCC(也就是反图的入度=0的SCC)。答案必然在这样的SCC中,所以只需要遍历每个入度=0的SCC,然后DFS计算能够到达的所有SCC,然后记得初始-1,因为本人不能给自己投票。

具体:

1.求SCC需要DFS得到后序序列,然后从后序序列的最后一个点往前DFS,如过遍历点时没有标记就说明是新的SCC,那些在逆后序序列中可以
从一个点访问到的都是同一个SCC。
2.遍历图中所有边,当边的两个顶点不属于同一SCC,from点[反图就是to点]所属的SCC出度[反图的入度]+1。
3.遍历所有点i,当 c[i] (所属的SCC)的入度InDeg[c[i]]=0,开始遍历寻找能够到达的所有SCC j,然后把他们的点数scc[j]都加到ANS[c[i]]
中,每次遍历完成后更新最大值。
4.遍历所有点,当对应的SCC的ANS值等于最大值时,打印

Codes

在这里插入图片描述

#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>
using namespace std;
const int MAXN=5005;
const int MAXM=30005;
int N,M;
struct Edge{
	int from,to,next;
}e1[MAXM],e2[MAXM];
int head1[MAXN],head2[MAXN],dfn[MAXN];
int scnt,c[MAXN],scc[MAXN],Out[MAXN],InDeg[MAXN];//scnt=scc数量,c[i]表示第i号同学处于c[i]号SCC中 InDeg是SCC入
int dcnt,tot; bool vis[MAXN],vic[MAXN];
int ans[MAXN];
void init(){
	tot=0;
	memset(head1,-1,sizeof(head1));
	memset(head2,-1,sizeof(head2));
    memset(Out,0,sizeof(Out));
}
void addEdge(int u,int v){ 
	e1[tot].from=u;e1[tot].to=v;   e1[tot].next=head1[u];  //正图 
	head1[u]=tot; 
	e2[tot].from=v; e2[tot].to=u;   e2[tot].next=head2[v];  //反图 
	head2[v]=tot; ++Out[u];  tot++;
}
void dfs1(int x){
	vis[x]=true;
	for(int i=head1[x];i!=-1;i=e1[i].next){
		if(!vis[e1[i].to]) dfs1(e1[i].to); 
    }
	dfn[dcnt++]=x;
}
void dfs2(int x){
	c[x]=scnt; ++scc[scnt];//第scnt号scc加一 
	for(int i=head2[x];i!=-1;i=e2[i].next){
	 if(!c[e2[i].to]) dfs2(e2[i].to);	
	}
}
void dfs(int x,int j){//j表示 源 强连通分量 
   vis[x]=true; 
   if(!vic[c[x]]) ans[j]+= scc[c[x]]; vic[c[x]]=true;
   for(int i=head2[x];i!=-1;i=e2[i].next){
   	 if(!vis[e2[i].to]) dfs(e2[i].to,j);
   }
}
void getCase(){
	dcnt=scnt=0;
	memset(scc,0,sizeof(scc));
	memset(vis,false,sizeof(vis));
	memset(c,0,sizeof(c)); 
	memset(InDeg,0,sizeof(InDeg)); 
	memset(ans,-1,sizeof(ans));
	for(int i=0;i<N;++i){
		if(!vis[i]) dfs1(i);
	} 
	for(int i=N-1;i>=0;--i){ //逆后序 
	   if(!c[dfn[i]]){
	   	++scnt; dfs2(dfn[i]);  //得到scnt个scc,i点的scc编号是c[i] 
	   }
	} 
	for(int i=0;i<M;++i){
	  int fr=e1[i].from,to=e1[i].to;
	  if(c[fr]!=c[to])  ++InDeg[c[fr]];
	} 
    int j=0,maxm=-1;
	for(int i=0;i<N;++i){
    	j=c[i];
		if(!InDeg[j]){
		memset(vis,false,sizeof(vis));
		memset(vic,false,sizeof(vic));	
		  dfs(i,j); InDeg[j]=1;
		  if(maxm<ans[j]) {  maxm=ans[j]; }
		}
	}
	cout<<maxm<<endl;bool tag=false;
	for(int i=0;i<N;++i){
		if(ans[c[i]]==maxm){
		 if(tag) cout<<" "; tag=true;
		 cout<<i;	
		} 
	} cout<<endl;
}
int main(){
	int T,A,B; scanf("%d",&T); 
	for(int i=1;i<T+1;++i){
	    init();
		scanf("%d%d",&N,&M);
		for(int j=0;j<M;++j){
			scanf("%d%d",&A,&B);
			addEdge(A,B);
		}
		cout<<"Case "<<i<<": ";
		getCase();
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值