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

A 区间选点（再谈）——差分约束方法

题目

Input&&&Output

Sample

#input：
5
3 7 3
8 10 3
6 8 1
1 3 1
10 11 1
#output：
6

题解

1.差分约束：对于不等式转化为>=或<=通过最长路或最短路求解一组解
2.本题对于区间来说有三个约束：
sum[b] - sum[a-1] >= c
sum[i] - sum[i-1] >=0
sum[i] - sum[i-1] <=1
3.由于是求最小解也就是求下界因此我们选择转化为>=然后跑最长路求解下界

C++代码

#include <iostream>
#include <queue>
using namespace std;


const int maxn = 1e6+500;
int n,a,b,c;
int inf=-maxn;

struct edge
{
	int to,next,w;
}e[maxn];
int head[maxn],tot,sum[maxn],maxb=0,mina=maxn; 
bool vis[maxn];

void init(){
	for(int i=0;i<maxn;i++)
		head[i]=0,sum[i]=inf,vis[i]=0;
	tot=0;
}

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

void spfa(int s){
	queue<int> q;
	q.push(s);
	sum[s]=0,vis[s]=1;
	while(!q.empty()){
		int u=q.front();
		q.pop();
		vis[u]=0;
		for(int i=head[u];i!=0;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]=1;
				}
				
			}
		}
	}
}

int main()
{
	cin>>n;
	init();
	for(int i=0;i<n;i++){
		cin>>a>>b>>c;
		maxb=max(maxb,b);
		mina=min(mina,a);
		//a-1->b c
		add(a-1,b,c);
	}
	//0<=sum[i]-sum[i-1]<=1 i->i+1 -1, i-1->i 0 
	for(int i=mina;i<=maxb;i++)
		add(i-1,i,0),add(i,i-1,-1);
	spfa(mina-1);
	int result = sum[mina];
	for(int i=mina;i<=maxb;i++) result = max(result,sum[i]);
	cout<<result<<endl;
}

B 猫猫向前冲

题目

Input&&&Output

Sample

#input:
4 3
1 2
2 3
4 3
#output：
1 2 4 3

题解

1.本题要求字典序最小的拓扑序列
2.首先对于求拓扑序列，实际上我们就是每次找入度为0的点然后更新它的附近节点并
将入度变为零的加入队列
3.由于是最小字典序，因此我们要保证每次取点取得都是最小的因此选用小根堆，由
与小根堆定义较为复杂因此采用往大根堆中加负数模拟小根堆

C++代码

#include<iostream>
#include<queue>
#include<vector>
using namespace std;

const int maxn=600; 
int n,m,ang[maxn];
vector<int> e[maxn];//邻接表存储 
vector<int> result;

void init(int n){
	for(int i=1;i<=n;i++){
		ang[i]=0;
		e[i] = vector<int>(0);
	}
}
void solve(){
	priority_queue<int> q;//大根堆 加复数模拟小根堆 
	for(int i=1;i<=n;i++) if(ang[i]==0) q.push(-i);
	while(!q.empty()){
		int curr=-q.top();q.pop();
		result.push_back(curr);
		for(int i=0;i<e[curr].size();i++){
			int nex=e[curr][i];
			ang[nex]--;
			if(ang[nex]==0)
				q.push(-nex);
		}
	}
}
int main()
{
	int x,y;
	while(cin>>n>>m){
		init(n);
		result = vector<int>(0);
		for(int i=0;i<m;i++){
			cin>>x>>y;
			ang[y]++;
			e[x].push_back(y);
		}
		solve();
		for(int i=0;i<result.size();i++)
			i==0 ? cout<<result[i] : cout<<" "<<result[i];
		cout<<endl;
	}
	
} 

C 班长竞选

题目

Input&&&Output


Sample

#input：
2
4 3
3 2
2 0
2 1

3 3
1 0
2 1
0 2
#output：
Case 1: 2
0 1
Case 2: 2
0 1 2

题解

1.首先我们发现投票结果具有传递性，而且可能会有环路产生，因此我们可以认为同
一个团体里的人投票效果是等效的
2.团体：此处我们将可以互相到达的一组点称为一个团体（也就是SCC
3.然后我们将团体视为一个巨型点（缩点）
4.如何求解投票结果：（此处对于团体而言）：
若a->b 则b的支持人数为 a人数+b人数-1
按照这种思路，此处有一个巧妙的做法，也就是求反图然后dfs就可以得到它对应的
人数
5.最后的优化：我们发现很多点是没有必要dfs的，因为既然是最大那么对于原图而
言，它的出度一定为零，因此对于反图我们只需要遍历入度为零的点然后dfs就可以
了。

C++代码

#include<iostream>
#include<vector>
#include<algorithm>
#include<cstring>
#include<string>
using namespace std;
const int maxn = 5005;
int g,n,m,cou[maxn],ang[maxn];
int dfn[maxn],dcnt,scnt,c[maxn],vis[maxn];//c[i] 属于哪个分量 scnt记录scc dfn记录后序的第i个点  dcnt fs时间 
vector<int> ed[maxn];//原图 
vector<int> red[maxn];//反图 
//vector<int> tiny[maxn];//缩点图 
vector<int> tiny2[maxn];//缩点图 反图 
vector<int> result;

void init(int n){
	result.clear();
	for(int i=0;i<=n;i++)
		ed[i].clear(),red[i].clear(),tiny2[i].clear();
}

void dfs1(int x){
	vis[x] = 1;
	for(int i=0;i<ed[x].size();i++)
		if(!vis[ed[x][i]]) dfs1(ed[x][i]);
	dfn[++dcnt] = x; 
}

void dfs2(int x){
	c[x] = scnt;
	for(int i=0;i<red[x].size();i++)
		if(!c[red[x][i]]) dfs2(red[x][i]);
}

void kosaraju(){
	dcnt = scnt = 0;
	memset(c,0,sizeof c);
	memset(vis,0,sizeof vis);
	for(int i=0;i<n;i++) if(!vis[i]) dfs1(i);
	for(int i=n;i>=1;i--) if(!c[dfn[i]]) ++scnt,dfs2(dfn[i]); 
}

int sum=-1,msum=0;
void dfs(int x){
	if(vis[x] == 1) return;
	vis[x] = 1;
	msum += cou[x];
	for(int i=0;i<tiny2[x].size();i++)
		if(!vis[tiny2[x][i]]) dfs(tiny2[x][i]);
}

void solve(){
	kosaraju();
	for(int i=1;i<=scnt;i++) cou[i]=0,ang[i]=0;
	for(int i=0;i<n;i++) cou[c[i]]++;//连通分量含点数
	for(int i=0;i<n;i++){
		for(int j=0;j<ed[i].size();j++){
			if(c[i] != c[ed[i][j]]){
				//tiny[c[i]].push_back(c[ed[i][j]]);
				//cout<<"*"<<c[i]<<" : "<<c[ed[i][j]]<<endl;
				tiny2[c[ed[i][j]]].push_back(c[i]);
				ang[c[i]]++;
			}
		}
	}
	//for(int i=1;i<=scnt;i++) cout<<i<<" : "<<ang[i]<<endl;
	sum = -1;
	for(int i=1;i<=scnt;i++){//找到入度为零的点 ->原图出度为零 
		if(ang[i] == 0){
			//dfs缩点后的反图
			msum = 0;
			memset(vis, 0, sizeof vis);
			dfs(i);
			if(msum > sum){//更新 sum
				sum = msum; 
				result.clear();
				result.push_back(i);
			}else if(msum == sum){
				result.push_back(i);
			}
		}
	}
	//for(int i=0;i<result.size();i++) cout<<"****"<<i<<" : "<<result[i]<<endl; 
}

int main(){
	int A,B;
	cin>>g;
	int gr=1;
	while(g--){
		cin>>n>>m;
		init(n);
		for(int i=0;i<m;i++){
			cin>>A>>B;
			ed[A].push_back(B);
			red[B].push_back(A);
		}
		solve();
		int flag=0;
		cout<<"Case "<<gr++<<": "<<sum-1<<endl;
		for(int i = 0;i<n;i++){
			for(int j =0;j<result.size();j++){
				if(c[i] == result[j]){
					if(flag == 0) flag = 1;
					else cout<<" ";
					cout<<i;
				}
			}	
		}
		cout<<endl;
	}
	return 0;
}