jzoj3577 [CEOI2011]Traffic 强连通分量缩点

Description

---

Gdynia 的中心坐落于Kacza 河中游的一个岛屿上。每天早晨，成千辆小车从河流的西岸的住宅区出发穿过岛屿（使用连接岛西路口的桥梁）到达东岸的工业区（使用连接岛东路口的桥梁）。

这个岛屿酷似一个边平行于坐标轴的矩形。因此，我们视他为一个笛卡尔坐标系上的A*B的矩形，对角的坐标为(0,0) 和(A,B)。

在这个岛屿上，有n 个路口从1 到n 编号。路口i 位于坐标(xi; yi)。如果一个路口坐标类似于(0,y)，那么它在岛西。类似的，坐标形于(A,y) 的路口坐落在岛东。路口通过街道连接起来。每一条街道是一条连接两个路口的线段。街道可以是双向或者单向的。没有两条街道有公共点（除了作为线段端点的路口）。岛上没有桥或者地道。其他的道路网络形状是不被认可的。特别的是，街道可以与岛屿的边缘重合，或者存在没有连接街道的路口。

因为交通密度不断增长，所以市长雇佣你去检查现在的道路网络是否足够。他要求你写一个程序，确定从每个岛西的路口出发能到达多少个岛东的路口。

Solution

---

注意到一个条件是没有两条边相交，那么我们把东岸的点按y坐标降序排列后，任意一个西岸的点能到达的东岸一定是连续的一段。tarjan缩点之后求出一个起点能到达的最大和最小标号的终点即可

本来还以为会爆栈的= =!

Code

---

#include <stdio.h>#include <string.h>#include <stack>#include <algorithm>#include <queue>#define rep(i,st,ed) for (int i=st;i<=ed;++i)#define fill(x,t) memset(x,t,sizeof(x))#define read2(x,y) read(x),read(y)#define read3(x,y,z) read(x),read(y),read(z)#define max(x,y) ((x)>(y)?(x):(y))#define min(x,y) ((x)<(y)?(x):(y))#define INF 0x3f3f3f3f#define N 600001#define E 4000001struct edge{int x,y,next;}e[E];struct pos{int x,y;}p[N];std:: queue<int>que;std:: stack<int>stack;int id[N],inStack[N],dfn[N],low[N],scc[N],data[N],vis[N];int mn[N],mx[N];int ls[N],edCnt=0;int n,m,A,B;void read(int &x) {    x=0; int v=1; char ch=getchar();    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());    x*=v;}void addEdge(int x,int y) {    e[++edCnt]=(edge){x,y,ls[x]}; ls[x]=edCnt;}void dfs(int now) {    if (vis[now]) return ;    vis[now]=1;    for (int i=ls[now];i;i=e[i].next) dfs(e[i].y);}void tarjan(int now) {    stack.push(now);    inStack[now]=1;    dfn[now]=low[now]=++dfn[0];    for (int i=ls[now];i;i=e[i].next) {        if (!dfn[e[i].y]) {            tarjan(e[i].y);            low[now]=min(low[now],low[e[i].y]);        } else if (!scc[e[i].y]) {            low[now]=min(low[now],dfn[e[i].y]);        }    }    if (dfn[now]==low[now]) {        ++scc[0];        mx[scc[0]]=-INF;        mn[scc[0]]=INF;        for (int tmp=0;tmp!=now;) {            tmp=stack.top(); stack.pop();            if (id[tmp]) {                mx[scc[0]]=max(mx[scc[0]],id[tmp]);                mn[scc[0]]=min(mn[scc[0]],id[tmp]);            }            scc[tmp]=scc[0];            inStack[tmp]=0;        }    }}bool cmp(int a,int b) {return p[a].y>p[b].y;}void dp(int now) {    if (vis[now]) return ;    vis[now]=1;    for (int i=ls[now];i;i=e[i].next) {        dp(e[i].y);        mx[now]=max(mx[now],mx[e[i].y]);        mn[now]=min(mn[now],mn[e[i].y]);    }}int main(void) {    read2(n,m); read2(A,B);    rep(i,1,n) {read2(p[i].x,p[i].y);}    rep(i,1,m) {        int x,y,opt;        read3(x,y,opt);        addEdge(x,y);        if (opt!=1) {addEdge(y,x);}    }    int cnt=0;    rep(i,1,n) if (!p[i].x) dfs(i);    rep(i,1,n) if (p[i].x==A&&vis[i]) data[++cnt]=i;    std:: sort(data+1,data+cnt+1,cmp);    rep(i,1,cnt) {id[data[i]]=i;}    rep(i,1,n) if (!dfn[i]) {tarjan(i);}    fill(ls,0); int tmp=edCnt;    rep(i,1,tmp) if (scc[e[i].x]!=scc[e[i].y]) {        addEdge(scc[e[i].x],scc[e[i].y]);    }    fill(vis,0);    rep(i,1,scc[0]) dp(i);    cnt=0;    rep(i,1,n) if (!p[i].x) {data[++cnt]=i;}    std:: sort(data+1,data+cnt+1,cmp);    rep(i,1,cnt) {        printf("%d\n", max(0,mx[scc[data[i]]]-mn[scc[data[i]]]+1));    }    return 0;}