2387: [Ceoi2011]Traffic

题目链接

题目大意：平面图，保证边不在除端点外的地方相交，求每个左侧点(x=0)能到达的右侧点(x=A)的数量

题解：删掉不可达的右侧点，容易看出每个左侧点可达的右侧点区间连续，只需要求min和max
然后缩点+dp（或者在反图上按Y递增和递减从右侧染色两次，记录时间差）

我的收获：……

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <vector>
#include <stack>
#include <algorithm>
using namespace std;

const int INF=1e9;
const int N=300005;

int n,m,A,B,remain;
int id[N];

int tim,scnt;
int dfn[N],low[N];
int col[N];
int mi[N],mx[N];
bool ins[N],vis[N];

stack<int> s;
vector<int> E[N],G[N];

struct point{
    int x,y,id;
}a[N],b[N];

bool cmp(point a,point b){return a.y>b.y;}

void trajan(int x)
{
    int now=0;
    low[x]=dfn[x]=++tim;
    s.push(x);ins[x]=1;
    for(int i=0;i<E[x].size();i++){
        int v=E[x][i];
        if(!dfn[v]) trajan(v),low[x]=min(low[x],low[v]);
        else if(ins[v]) low[x]=min(low[x],dfn[v]);
    }
    if(low[x]==dfn[x]){
        scnt++;mx[scnt]=-INF;mi[scnt]=INF;
        while(now!=x){
            now=s.top();s.pop();
            col[now]=scnt,ins[now]=0;
            if(a[now].x==A) mx[scnt]=max(mx[scnt],id[now]),mi[scnt]=min(mi[scnt],id[now]);
        }
    }
} 

void rebuild()
{
    for(int x=1;x<=n;x++)
        for(int i=0;i<E[x].size();i++)
            if(col[x]!=col[E[x][i]]) G[col[x]].push_back(col[E[x][i]]);
    for(int x=1;x<=scnt;x++) sort(G[x].begin(),G[x].end());
    for(int x=1;x<=scnt;x++) G[x].erase(unique(G[x].begin(),G[x].end()),G[x].end());
}

void dp(int x)
{
    if(vis[x]) return ;
    vis[x]=1;
    for(int i=0;i<G[x].size();i++){
        int v=G[x][i];
        dp(v);
        mx[x]=max(mx[x],mx[v]);
        mi[x]=min(mi[x],mi[v]);
    }
} 

void dfs(int x)
{
    if(vis[x]) return ;
    vis[x]=1;
    for(int i=0;i<E[x].size();i++) dfs(E[x][i]);
}

void delit()
{
    for(int i=1;i<=n;i++) if(!a[i].x) dfs(a[i].id);
    for(int i=1;i<=n;i++) if(a[i].x==A&&vis[i]) b[++remain]=a[i];
    sort(b+1,b+1+remain,cmp);
    for(int i=1;i<=remain;i++) id[b[i].id]=i;
}

void output()
{
    sort(a+1,a+1+n,cmp);
    for(int i=1;i<=n;i++) if(!a[i].x) printf("%d\n",max(0,mx[col[a[i].id]]-mi[col[a[i].id]]+1));
}

void work(){
    delit();
    memset(vis,0,sizeof(vis));
    for(int i=1;i<=n;i++) if(!dfn[i]) trajan(i);
    rebuild();
    for(int i=1;i<=scnt;i++) dp(i);
    output();
}

void init()
{
    cin>>n>>m>>A>>B;
    for(int i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y),a[i].id=i;
    while(m--){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        E[x].push_back(y);
        if(z==2) E[y].push_back(x);
    }
}

int main()
{
    init();
    work();
    return 0;
}