解析
比起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
*/