Tarjan(强连通分量缩点) - SCOI 2011 - 糖果 - 洛谷 P3275

Tarjan(强连通分量缩点) - SCOI 2011 - 糖果 - 洛谷 P3275

幼儿园里有 N 个小朋友,老师现在想要给这些小朋友们分配糖果,要求每个小朋友都要分到糖果。

但是小朋友们也有嫉妒心,总是会提出一些要求,比如小明不希望小红分到的糖果比他的多,于是在分配糖果的时候, 老师需要满足小朋友们的 K 个要求。

幼儿园的糖果总是有限的,老师想知道他至少需要准备多少个糖果,才能使得每个小朋友都能够分到糖果,并且满足小朋友们所有的要求。

输入格式

输入的第一行是两个整数 N,K。

接下来 K 行,表示分配糖果时需要满足的关系,每行 3 个数字 X,A,B。

如果 X=1.表示第 A 个小朋友分到的糖果必须和第 B 个小朋友分到的糖果一样多。

如果 X=2,表示第 A 个小朋友分到的糖果必须少于第 B 个小朋友分到的糖果。

如果 X=3,表示第 A 个小朋友分到的糖果必须不少于第 B 个小朋友分到的糖果。

如果 X=4,表示第 A 个小朋友分到的糖果必须多于第 B 个小朋友分到的糖果。

如果 X=5,表示第 A 个小朋友分到的糖果必须不多于第 B 个小朋友分到的糖果。

小朋友编号从 1 到 N。

输出格式

输出一行,表示老师至少需要准备的糖果数,如果不能满足小朋友们的所有要求,就输出 −1。

数据范围

1 ≤ N < 1 0 5 , 1 ≤ K ≤ 1 0 5 , 1 ≤ X ≤ 5 , 1 ≤ A , B ≤ N 1≤N<10^5, 1≤K≤10^5, 1≤X≤5, 1≤A,B≤N 1N<105,1K105,1X5,1A,BN

输入样例:

5 7
1 1 2
2 3 2
4 4 1
3 4 5
5 4 5
2 3 5
4 5 1

输出样例:

11

分析:

本 题 有 差 分 约 束 解 法 : 本题有差分约束解法: 《差分约束 - SCOI 2011 - 糖果 - 洛谷 P3275》

这 里 通 过 t a r j a n 缩 点 来 解 决 。 这里通过tarjan缩点来解决。 tarjan

本 题 要 求 最 小 值 , 因 此 是 跑 最 长 路 。 本题要求最小值,因此是跑最长路。

跑 完 最 长 路 , 有 不 等 关 系 d i s [ j ] ≥ d i s [ u ] + w [ i ] 。 跑完最长路,有不等关系dis[j]≥dis[u]+w[i]。 dis[j]dis[u]+w[i]

故 将 所 有 不 等 关 系 转 化 为   ′ ≥ ′   号 。 故将所有不等关系转化为\ '≥'\ 号。   

① 、 X = 1 : A = B 即 , A ≥ B 且 B ≥ A < = > 添 加 B − > A 和 A − > B 两 条 边 , 边 权 均 为 0 。 ①、X=1:A=B即,A≥B且B≥A\quad<=>\quad添加B->A和A->B两条边,边权均为0。 X=1A=BABBA<=>B>AA>B0

② 、 X = 2 : A < B 即 , B ≥ A + 1 < = > 添 加 A − > B 一 条 边 , 边 权 为 1 。 ②、X=2:A<B即,B≥A+1\quad<=>\quad添加A->B一条边,边权为1。 X=2A<BBA+1<=>A>B1

③ 、 X = 3 : A ≥ B 即 , A ≥ B < = > 添 加 B − > A 一 条 边 , 边 权 为 0 。 ③、X=3:A≥B即,A≥B\quad<=>\quad添加B->A一条边,边权为0。 X=3ABAB<=>B>A0

④ 、 X = 4 : A > B 即 , A ≥ B + 1 < = > 添 加 B − > A 一 条 边 , 边 权 为 1 。 ④、X=4:A>B即,A≥B+1\quad<=>\quad添加B->A一条边,边权为1。 X=4A>BAB+1<=>B>A1

⑤ 、 X = 5 : A ≤ B 即 , B ≥ A < = > 添 加 A − > B 一 条 边 , 边 权 为 0 。 ⑤、X=5:A≤B即,B≥A\quad<=>\quad添加A->B一条边,边权为0。 X=5ABBA<=>A>B0

⑥ 、 要 求 每 个 小 朋 友 都 要 分 到 糖 果 , 即 x i ≥ 1 , x i 表 示 第 i 个 小 朋 友 分 得 的 糖 果 数 量 。 ⑥、要求每个小朋友都要分到糖果,即x_i≥1,x_i表示第i个小朋友分得的糖果数量。 xi1xii

下面解释为何可以使用tarjan缩点来解决:

本 题 所 有 的 边 权 均 为 0 或 1 。 本题所有的边权均为0或1。 01

我 们 的 第 一 个 目 标 是 判 断 是 否 有 解 , 即 判 断 是 否 有 正 环 。 我们的第一个目标是判断是否有解,即判断是否有正环。

通 过 t a r j a n 算 法 缩 点 后 , 在 拓 扑 图 中 , 正 环 必 定 是 存 在 于 某 个 强 连 通 分 量 中 。 通过tarjan算法缩点后,在拓扑图中,正环必定是存在于某个强连通分量中。 tarjan

由 于 所 有 边 权 均 为 0 或 1 , 因 此 在 一 个 强 连 通 分 量 中 , 只 要 存 在 一 条 边 权 大 于 0 的 边 , 就 会 存 在 正 环 。 由于所有边权均为0或1,因此在一个强连通分量中,只要存在一条边权大于0的边,就会存在正环。 010

若 不 存 在 正 环 , 说 明 每 个 强 连 通 分 量 中 的 边 权 均 为 0 , 若不存在正环,说明每个强连通分量中的边权均为0, 0

那 么 对 于 同 一 个 强 连 通 分 量 而 言 , 连 通 分 量 内 部 所 有 的 点 到 拓 扑 图 起 点 的 距 离 均 相 同 。 那么对于同一个强连通分量而言,连通分量内部所有的点到拓扑图起点的距离均相同。

这 样 , 我 们 就 能 够 把 一 个 强 连 通 分 量 看 作 一 个 整 体 , 当 作 拓 扑 图 中 的 一 个 点 , 来 跑 最 短 路 。 这样,我们就能够把一个强连通分量看作一个整体,当作拓扑图中的一个点,来跑最短路。

第 二 个 目 标 是 求 所 有 点 到 源 点 的 距 离 之 和 。 第二个目标是求所有点到源点的距离之和。

注 意 , 求 距 离 之 和 时 , 要 拆 点 。 每 个 连 通 分 量 中 的 点 的 距 离 都 要 累 加 上 , 故 我 们 需 要 统 计 每 个 连 通 分 量 中 点 的 数 量 。 注意,求距离之和时,要拆点。每个连通分量中的点的距离都要累加上,故我们需要统计每个连通分量中点的数量。

总结步骤:

① 、 建 图 。 ①、建图。

② 、 缩 点 。 ②、缩点。

③ 、 建 图 ( 拓 扑 图 ) , 并 判 断 是 否 有 环 。 ③、建图(拓扑图),并判断是否有环。 ()

④ 、 拆 点 , 求 和 。 ④、拆点,求和。

代码:

#include<cstdio>
#include<cstring>
#include<algorithm>

#define ll long long

using namespace std;

const int N=100010, M=600010;

int n,m;
int h[N],hs[N],e[M],ne[M],w[M],idx;
int stk[N],top;
bool in_stk[N];
int dfn[N],low[N],Size[N],timestamp;
int dis[N];
int id[N],scc_cnt;

void add(int h[],int a,int b,int c)
{
    e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}

void tarjan(int u)
{
    dfn[u]=low[u]=++timestamp;
    stk[++top]=u,in_stk[u]=true;
    for(int i=h[u];~i;i=ne[i])
    {
        int j=e[i];
        if(!dfn[j])
        {
            tarjan(j);
            low[u]=min(low[u],low[j]);
        }
        else if(in_stk[j]) low[u]=min(low[u],dfn[j]);
    }

    if(dfn[u]==low[u])
    {
        ++scc_cnt;
        int y;
        do
        {
            y=stk[top--];
            in_stk[y]=false;
            id[y]=scc_cnt;
            Size[scc_cnt]++;
        }while(y!=u);
    }
}

bool build()
{
    for(int i=0;i<=n;i++)   //建图
        for(int j=h[i];~j;j=ne[j])
        {
            int k=e[j];
            int a=id[i],b=id[k];
            if(a==b && w[j]>0) return false;    //同意连通块中权值大于0
            else if(a!=b) add(hs,a,b,w[j]);
        }
        
    return true;
}

ll solve()
{
    for(int i=scc_cnt;i;i--)
        for(int j=hs[i];~j;j=ne[j])
        {
            int k=e[j];
            dis[k]=max(dis[k],dis[i]+w[j]);
        }
        
    ll res=0;
    for(int i=1;i<=scc_cnt;i++) res+=dis[i]*Size[i];
    return res;
}


int main()
{
    scanf("%d%d",&n,&m);
    
    memset(h,-1,sizeof h);
    memset(hs,-1,sizeof hs);
    for(int i=1;i<=n;i++) add(h,0,i,1);  //xi>=x0+1
    
    int t,a,b;
    for(int i=0;i<m;i++)
    {
        scanf("%d%d%d",&t,&a,&b);
        if(t==1) add(h,a,b,0),add(h,b,a,0);  //A=B
        else if(t==2) add(h,a,b,1);   //A<B  =>  B>=A+1
        else if(t==3) add(h,b,a,0);   //A>=B
        else if(t==4) add(h,b,a,1);   //A>B  =>  A>=B+1
        else if(t==5) add(h,a,b,0);   //B>=A
    }
    
    tarjan(0);  //缩点
    
    if(!build()) puts("-1");
    
    else printf("%lld\n",solve());
    
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值