cf949C 建模,SCC缩点

/*
给定n个数据中心,m份资料,每份资料在其中的两个中心备份,一天可供下载的时间是h小时
中心i在第hi小时需要维护,无法下载
现在要将一些中心的维护时间往后推1小时,使得任意时刻每份资料都可以被下载,请问最少选择多少个数据中心, 某个中心维护时,在其中资料无法下载,必须到其他点下载,
如果该点对应的点也在维护,那么这个对应点的维护必须往后推
对应点往后推时继续和其余点矛盾,那么其余点也要往后 所以,如果同一份数据的两个点相隔一小时,那么这两个点要么不推迟,要推迟就要一起推迟, 建立模型:点x和点y有相同的资料,并且hx+1==hy,那么加有向边(x,y)
最后得到的图中,一个强连通分量内的所有点必须一起选择
缩点构成DAG,选择出度为0的点即可
*/
#include<bits/stdc++.h>
using namespace std;
#define maxn 100005
struct Edge{int to,nxt;}edge[maxn<<];
int n,m,t,h[maxn],head[maxn],tot,d[maxn][]; int c[maxn],low[maxn],dfn[maxn],stk[maxn],ins[maxn],ind,top;
int cnt;
vector<int>scc[maxn];
void tarjan(int x){
dfn[x]=low[x]=++ind;
stk[++top]=x;ins[x]=;
for(int i=head[x];i!=-;i=edge[i].nxt){
int y=edge[i].to;
if(!dfn[y]){//树枝边
tarjan(y);
low[x]=min(low[x],low[y]);
}
else if(ins[y])//后向边
low[x]=min(low[x],dfn[y]);
}
if(dfn[x]==low[x]){
cnt++;
int y;
do{
y=stk[top--];
ins[y]=;
c[y]=cnt;scc[cnt].push_back(y);
}while(x!=y);
}
} void init(){
tot=;
memset(head,-,sizeof head);
}
void addedge(int u,int v){
edge[tot].nxt=head[u];
edge[tot].to=v;
head[u]=tot++;
}
int in[maxn],out[maxn];
int main(){
init();
cin>>n>>m>>t;
for(int i=;i<=n;i++)scanf("%d",&h[i]);
for(int i=;i<=m;i++)scanf("%d%d",&d[i][],&d[i][]);
for(int i=;i<=m;i++){//建图
if((h[d[i][]]+)%t==h[d[i][]])
addedge(d[i][],d[i][]);
if((h[d[i][]]+)%t==h[d[i][]])
addedge(d[i][],d[i][]);
}
for(int i=;i<=n;i++)
if(!dfn[i])
tarjan(i);
for(int x=;x<=n;x++)
for(int i=head[x];i!=-;i=edge[i].nxt){
int y=edge[i].to;
if(c[x]==c[y])continue;
in[c[y]]++;out[c[x]]++;
}
int ans=0x3f3f3f3f,k;
for(int i=;i<=cnt;i++)
if(out[i]== && scc[i].size()<ans)
k=i,ans=scc[i].size();
cout<<ans<<endl;
for(int i=;i<scc[k].size();i++)
printf("%d ",scc[k][i]);
}
上一篇:[Reinforcement Learning] Policy Gradient Methods


下一篇:MySQL 数据库慢查询日志分析脚本