CF1137C:Museums Tour（缩点、分层图）

解析

比起CF，更像是主题库的一道题
也是很好的一道题

由于可以在环里转圈，容易想到缩点
但是出现了一个问题：缩点是无法统计一个强连通分量里有哪些环及其大小的
而本题中显然环的大小和d的关系至关重要

考虑转换思路
把每个城市x割成d个点(x,1)到(x,d)
每个点表示第i天到达该城市
对新图跑tarjan
每个分量记录其不同的博物馆的数目
可以利用一个带标记的桶来解决

最后跑一遍简单的dp即可

为什么这样不会重复统计？
假设(x,i)走到了不在同一个分量里的(x,i+w)，使城市x被统计了两遍
那么x就在一个大小为w的环上
那么从(x,i+w)再在环上跑(d-1)遍，必然能回到(x,i)
与假设中“不在同一个分量矛盾”
得证

代码

#include<bits/stdc++.h>
const int N=1e7+100;
const int mod=1e9+7;
#define ll long long
using namespace std;
inline ll read() {
	ll x(0),f(1);char c=getchar();
	while(!isdigit(c)) {if(c=='-')f=-1;c=getchar();}
	while(isdigit(c)) {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
	return x*f;
}

int n,m,d;
#define id(i,j) ((i-1)*d+j)
struct node{
	int to,nxt;
}p[N];
int fi[N],cnt;
inline void addline(int x,int y){
	p[++cnt]=(node){y,fi[x]};fi[x]=cnt;
	return;
}
bool jd[N];
int bac[N];
int dfn[N],low[N],top,tim,tot,zhan[N],col[N];
inline void add(int x){
	if(jd[x]&&bac[(x+d-1)/d]!=tot){
		bac[(x+d-1)/d]=tot;low[tot]++;
	}
	return;
}
inline void tarjan(int x){
	dfn[x]=low[x]=++tim;zhan[++top]=x;
	for(int i=fi[x];~i;i=p[i].nxt){
		int to=p[i].to;
		if(!dfn[to]){
			tarjan(to);
			low[x]=min(low[x],low[to]);
		}
		else if(!col[to]) low[x]=min(low[x],dfn[to]);
	}
	if(low[x]==dfn[x]){
		col[x]=++tot;
		add(x);
		while(zhan[top]!=x){
			int now=zhan[top--];
			col[now]=tot;
			add(now);
		}
		top--;
	}
	return;
}
//int dp[N];
inline int find(int x){
	if(dfn[x]!=-1) return dfn[x];
	dfn[x]=low[x];
	for(int i=fi[x];~i;i=p[i].nxt){
		int to=p[i].to;
		dfn[x]=max(dfn[x],low[x]+find(to));
	}
	return dfn[x];
}
int main(){
	//printf("%d\n",sizeof(jd)/1024/1024);
	//memset(dp,-1,sizeof(dp));
	memset(fi,-1,sizeof(fi));cnt=-1;
	n=read();m=read();d=read();
	for(int i=1;i<=m;i++){
		int x=read(),y=read();
		for(int j=1;j<=d;j++) addline(id(x,j),id(y,j%d+1));
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=d;j++){
			int o=id(i,j);
			scanf("%1d",&jd[o]);
		}
	}
	tot=id(n,d);
	for(int i=1;i<=id(n,d);i++){
		if(!dfn[i]) tarjan(i);
	}
	for(int i=1;i<=id(n,d);i++){
		for(int j=fi[i];~j;j=p[j].nxt){
			int to=p[j].to;
			if(col[i]!=col[to]) addline(col[i],col[to]);
		}
	}
	memset(dfn,-1,sizeof(dfn));
	printf("%d\n",find(col[id(1,1)]));
}
/*
1
281239
*/