队爷的讲学计划 CH Round #59 - OrzCC杯NOIP模拟赛day1

题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的讲学计划

题解:刚开始理解题意理解了好半天,然后发现好像就是tarjan之后求个最长链?

然后题目规定了从1出发,那么就只将scc[1]放入初始队列。

然后只有10分。。。

后来我发现了这样的情况:因为只有scc[1]放入了队列,所以其它入度为0的并没有被放入队列,这使得一些点的inp无法减到0,以致无法更新答案

然后我就先从scc[1]dfs了一遍,先统计出每个点可通过scc[1]到达的度数,再拓扑排序,就过了。

好题,怒赞!

代码:

 #include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 1000000000
#define maxn 550000
#define maxm 500+100
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define mod 1000000007
using namespace std;
inline int read()
{
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}
return x*f;
}
int n,m,cnt,ti,tot,head[maxn],inp[maxn],dfn[maxn],low[maxn],sta[maxn],top;
int scc[maxn],s[maxn],f[maxn],a[maxn],b[maxn],c[maxn];
ll g[maxn];
struct edge{int go,next,w;}e[maxn];
queue<int> q;
inline void insert(int x,int y,int z)
{
e[++tot].go=y;e[tot].next=head[x];e[tot].w=z;head[x]=tot;
}
inline void tarjan(int x)
{
dfn[x]=low[x]=++ti;sta[++top]=x;
for(int i=head[x],y;i;i=e[i].next)
if(!dfn[y=e[i].go])
{
tarjan(y);low[x]=min(low[x],low[y]);
}
else if(!scc[y])low[x]=min(low[x],dfn[y]);
if(low[x]==dfn[x])
{
cnt++;int now=;
while(now!=x)
{
now=sta[top--];scc[now]=cnt;s[cnt]++;
}
}
}
inline void dfs(int x)
{
for(int i=head[x];i;i=e[i].next)
{
int y=e[i].go;
inp[y]++;
if(inp[y]==)dfs(y);
}
}
int main()
{
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
n=read();m=read();
for1(i,m){a[i]=read();b[i]=read();c[i]=read();insert(a[i],b[i],c[i]);}
for1(i,n)if(!dfn[i])tarjan(i);
memset(head,,sizeof(head));tot=;
for1(i,m)if(scc[a[i]]!=scc[b[i]])insert(scc[a[i]],scc[b[i]],c[i]);
for1(i,cnt)f[i]=,g[i]=inf;
f[scc[]]=g[scc[]]=s[scc[]];g[scc[]]--;
dfs(scc[]);
q.push(scc[]);
//for1(i,cnt)if(!inp[i])q.push(i),f[i]=s[i],g[i]=s[i]-1;
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=head[x];i;i=e[i].next)
{
int y=e[i].go;
inp[y]--;
if(!inp[y])q.push(y);
if(f[x]+s[y]>f[y])f[y]=f[x]+s[y],g[y]=g[x]+e[i].w+s[y]-;
else if(f[x]+s[y]==f[y])g[y]=min(g[y],g[x]+e[i].w+s[y]-);
}
}
int ans=;
for1(i,cnt)
if(f[i]>f[ans])ans=i;
else if(f[i]==f[ans]&&g[i]<g[ans])ans=i;
printf("%d %lld\n",f[ans],g[ans]);
return ;
}
上一篇:Nginx下TIME_WAIT过多的调优


下一篇:JavaScript语言精粹读书笔记 - JavaScript函数