【COCI 2021-2022#5】Usmjeravanje

Description

Alice生活在梦幻岛,岛中有两条河流,均从西向东流动。第一条河的岸边有a座城市,从西到东依次编号为1到a。同样,第二条河的岸边有b个城市,从西到东依次编号为1到b。同一条河岸边的两个城市i和j,如果满足i

梦幻岛的居民计划建立m条单向航线。第i条航线连接第一条河流的城市xi和第二条河流的城市yj,但方向尚未确定。梦幻岛的居民希望把城市尽可能地连接起来。

一个城市对(c,d),如果从城市c出发能到达城市d,且从城市d出发能到达城市c,则称城市对(c,d)是互连的,途中航线和河流都允许被使用。岛主找到Alice请她帮忙确定每条航线的方向,使得最大的任意两个城市均不互连的城市群最小。

Solution

看到任意两点均不互连,想到与之类似的强连通分量,预示着这题将使用 Tarjan \text{Tarjan} Tarjan 来完成。

先考虑把边的方向确定下来。

注意到有贡献的连法只有以下三种:

【COCI 2021-2022#5】Usmjeravanje_第1张图片 【COCI 2021-2022#5】Usmjeravanje_第2张图片 【COCI 2021-2022#5】Usmjeravanje_第3张图片

观察发现,这三种连法中右下的点一定是向左上连出的。

根据这个发现,我们有了以下思路:

  1. 将边排序,第一关键字 x 从小到大(上河岸),第二关键字 y 从大到小(下河岸)。

  2. 枚举每条边,如果当前边下河岸的点在目前所有边的最右端,那么说明这个点在某种连法中一定处在右下的位置,所以此边的类型为 1。

  3. 其余边类型为 0。

现在已经确定了边的方向,剩下的就是模板题:求强连通的个数。

Code

#include
#include
#include
#define N 400005
using namespace std;
int na,nb,n,m,tot,num,mx,cnt,ans,sum[N],sta[N*10],dfn[N<<1],low[N<<1];
bool bj[N<<1];
struct node {int x,y,id;} edg[N];
struct edge {int to,next,head;} a[N*3];
bool cmp(node x,node y)
{
    if (x.x<y.x) return true;
    if (x.x>y.x) return false;
    return x.y>y.y;
}
void add(int x,int y) {a[++tot].to=y;a[tot].next=a[x].head;a[x].head=tot;}
void dfs(int x)
{
    dfn[x]=low[x]=++cnt;
    sta[++num]=x;
    bj[x]=true;
    for (int i=a[x].head;i;i=a[i].next)
    {
        int y=a[i].to;
        if (!dfn[y]) dfs(y),low[x]=min(low[x],low[y]);
        else if (bj[y]) low[x]=min(low[x],dfn[y]);
    }
    if (dfn[x]==low[x])
    {
        ++ans;
        while (sta[num]!=x) bj[sta[num--]]=false;
        bj[sta[num--]]=false;
    }
}
int main()
{
    freopen("neverland.in","r",stdin);
    freopen("neverland.out","w",stdout);
    scanf("%d%d",&na,&nb);
    for (int i=1;i<na;++i)
        add(i,i+1);
    for (int i=1;i<nb;++i)
        add(i+na,i+na+1);
    scanf("%d",&m);
    for (int i=1;i<=m;++i)
        scanf("%d%d",&edg[i].x,&edg[i].y),edg[i].id=i;
    sort(edg+1,edg+m+1,cmp);
    for (int i=1;i<=m;++i)
    {
        if (edg[i].y>mx) mx=edg[i].y,add(edg[i].y+na,edg[i].x),sum[edg[i].id]=1;
        else add(edg[i].x,edg[i].y+na),sum[edg[i].id]=0;
    }
    n=na+nb;
    for (int i=1;i<=n;++i)
        if (!dfn[i]) dfs(i);
    printf("%d\n",ans);
    for (int i=1;i<=m;++i)
        printf("%d ",sum[i]);
    return 0;
}

你可能感兴趣的:(算法,深度优先,c++)