模板:
//模板
int scc_cnt,id[N],_size[N];
int dfn[N],low[N];
int stk[N],top;
bool in_stk[N];
int timestamp;
void tarjan(int u)
{
dfn[u]=low[u]=++timestamp;
stk[++top]=u,in_stk[u]=true;
for(int i=h[u];i!=-1;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],low[j]);
}
}
if(dfn[u]==low[u]){
int y;
scc_cnt++;
do{
y=stk[top--];
in_stk[y]=false;
id[y]=scc_cnt;
_size[scc_cnt]++;
}while(u!=y);
}
}
例题:
银河中的恒星浩如烟海,但是我们只关注那些最亮的恒星。
我们用一个正整数来表示恒星的亮度,数值越大则恒星就越亮,恒星的亮度最暗是 1。
现在对于 N 颗我们关注的恒星,有 M 对亮度之间的相对关系已经判明。
你的任务就是求出这 N𝑁 颗恒星的亮度值总和至少有多大。
输入格式
第一行给出两个整数 N 和 M。
之后 M 行,每行三个整数 T,A,B,表示一对恒星 (A,B) 之间的亮度关系。恒星的编号从 1 开始。
如果 T=1,说明 A 和 B 亮度相等。
如果 T=2,说明 A 的亮度小于 B 的亮度。
如果 T=3,说明 A 的亮度不小于 B 的亮度。
如果 T=4,说明 A 的亮度大于 B 的亮度。
如果 T=5,说明 A 的亮度不大于 B 的亮度。输出格式
输出一个整数表示结果。
若无解,则输出 −1。
数据范围
N≤100000,M≤100000
输入样例:
5 7 1 1 2 2 3 2 4 4 1 3 4 5 5 4 5 2 3 5 4 5 1
输出样例:
11
ac代码 :
这里有两种思路正解是强连通分量线性o(n)过,但acwing数据很弱可以差分约束spfa过,不过洛谷hack数据过不了,这里两种写法都给出来。
//差分约束
#define _CRT_SECURE_NO_WARNINGS
#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#include<stack>
using namespace std;
typedef long long LL;
const int N = 1e5 + 10,M=3*N;
int h[M], w[M], ne[M], e[M];
int d[N];
int cnt[N];
int idx;
int n, m;
bool st[N];
void add(int a, int b, int c)
{
e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx++;
}
bool spfa()
{
memset(d, -0x3f, sizeof(d));
stack<int> q;
q.push(0);
d[0] = 0;
st[0] = true;
while (q.size())
{
int t = q.top();
q.pop();
st[t] = false;
for (int i = h[t]; i != -1; i = ne[i])
{
int j = e[i];
if (d[j] < d[t] + w[i]) {
d[j] = d[t] + w[i];
cnt[j] = cnt[t] + 1;
if (cnt[j] > n) {
return false;
}
if (st[j])
{
continue;
}
q.push(j);
st[j] = true;
}
}
}
return true;
}
int main()
{
scanf("%d%d", &n, &m);
memset(h, -1, sizeof(h));
for (int i = 1; i <= m; i++)
{
int a, b, c;
scanf("%d%d%d", &a, &b, &c);
if (a == 1) {
add(c, b, 0);
add(b, c, 0);
}
else if (a == 2)
{
add(b, c, 1);
}
else if (a == 3)
{
add(c, b, 0);
}
else if (a == 4)
{
add(c, b, 1);
}
else {
add(b, c, 0);
}
}
for (int i = 1; i <= n; i++)
{
add(0, i, 1);
}
LL sum = 0;
if (spfa())
{
for (int i = 1; i <= n; i++)
{
sum += d[i];
}
printf("%lld", sum);
}
else {
printf("-1");
}
}
//强连通分量
#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10,M=6*N;
int h[N],hs[N],w[M],ne[M],e[M],idx;
int stk[N],scc_cnt,se[N],top,id[N];
int dfn[N],low[N];
int timestamp;
bool in_stk[N];
int d[N];
int n,m;
void add(int h[],int a,int b,int c)
{
w[idx]=c,e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void targan(int u)
{
dfn[u]=low[u]=++timestamp;
stk[++top]=u,in_stk[u]=true;
for(int i=h[u];i!=-1;i=ne[i]){
int j=e[i];
if(!dfn[j]){
targan(j);
low[u]=min(low[u],low[j]);
}else if(in_stk[j])
{
low[u]=min(low[u],low[j]);
}
}
if(dfn[u]==low[u])
{
scc_cnt++;
int y;
do{
y=stk[top--];
in_stk[y]=false;
id[y]=scc_cnt;
se[scc_cnt]++;
}while(y!=u);
}
}
int main()
{
scanf("%d%d",&n,&m);
memset(h,-1,sizeof(h));
memset(hs,-1,sizeof(hs));
for(int i=1;i<=m;i++)
{
int t,a,b;
scanf("%d%d%d",&t,&a,&b);
if(t==1)
{
add(h,a,b,0);
add(h,b,a,0);
}else if(t==2)
{
add(h,a,b,1);
}else if(t==3)
{
add(h,b,a,0);
}else if(t==4)
{
add(h,b,a,1);
}else if(t==5)
{
add(h,a,b,0);
}
}
for(int i=1;i<=n;i++){
add(h,0,i,1);
}
targan(0);
for(int i=0;i<=n;i++)
{
for(int j=h[i];j!=-1;j=ne[j]){
int k=e[j];
int a=id[i];
int b=id[k];
if(a==b){
if(w[j]){
printf("-1");
return 0;
}
}else {
add(hs,a,b,w[j]);
}
}
}
for(int i=scc_cnt;i;i--){
for(int j=hs[i];j!=-1;j=ne[j]){
int k=e[j];
d[k]=max(d[k],d[i]+w[j]);
}
}
long long ans=0;
for(int i=1;i<=scc_cnt;i++)
{
ans+=(long long)d[i]*se[i];
}
printf("%lld",ans);
return 0;
}