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的边加入,得到缩点之后的图
- 分析可得:对于属于第i个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;
}