NOI模拟 时间幻象

按照四个字符的个数分出点，连出一个图，缩点后用拓扑序求一个最长链。
忽略这个奇怪的组合数。

#include<bits/stdc++.h>
using namespace std;
const int maxn=32,maxm=55,N=maxn*maxn*maxn*maxn;
typedef long long ll;
queue<int> Q;
int P[10]={2,3,5,7,11,13,17,19,23,29};
int n,m,len,trans[maxm][4],cost[maxm][4];char s[maxn];
struct NUM{
	int num[10];
	NUM(){memset(num,0,sizeof num);}
	friend inline NUM operator*(const NUM &a,const NUM &b){
		NUM ret;
		for(int i=0;i<10;++i) ret.num[i]=a.num[i]+b.num[i];
		return ret;
	}
	friend inline NUM operator/(const NUM &a,const NUM &b){
		NUM ret;
		for(int i=0;i<10;++i) ret.num[i]=a.num[i]-b.num[i];
		return ret;
	}
	inline ll getval(ll ret=1){
		for(int i=0;i<10;++i)
			for(int j=1;j<=num[i];++j)
				ret*=P[i];
		return ret;
	}
}fac[maxn];
int Head[N],Next[N],V[N],cnt=0;ll val[N],ans=0;
int dfn[N],low[N],st[N],ins[N],tot=0,top=0;
int belong[N],in[N],num;ll siz[N],dis[N];
inline void add(int u,int v){Next[++cnt]=Head[u],V[cnt]=v,Head[u]=cnt;}
inline int id(int A,int B,int C,int D){return A*(n+1)*(n+1)*(n+1)+B*(n+1)*(n+1)+C*(n+1)+D;}
inline void tarjan(int u){
	dfn[u]=low[u]=++tot,st[++top]=u,ins[u]=1;
	for(int i=Head[u];i;i=Next[i]){
		if(!dfn[V[i]]) tarjan(V[i]),low[u]=min(low[u],low[V[i]]);
		else if(ins[V[i]]) low[u]=min(low[u],dfn[V[i]]);
	}
	if(dfn[u]==low[u]){
		++num;int k;do{
			k=st[top--],belong[k]=num,siz[num]+=val[k],ins[k]=0;
		}while(k!=u);
	}
}
int main(){
	//freopen("return.in","r",stdin);
	scanf("%d%d",&n,&m),num=(n+1)*(n+1)*(n+1)*n;
	for(int i=1;i<=n;++i){
		int tmp=i;NUM a;
		for(int j=0;j<10;++j)
			while(tmp%P[j]==0)
				a.num[j]+=1,tmp/=P[j];
		fac[i]=fac[i-1]*a;
	}
	for(int i=1;i<=m;++i){
		scanf("%s",s),len=strlen(s);
		for(int j=0;j<len;++j)
			--trans[i][s[j]-'A'],++cost[i][s[j]-'A'];
		scanf("%s",s),len=strlen(s);
		for(int j=0;j<len;++j)
			++trans[i][s[j]-'A'];
	}
	for(int i=0;i<=n;++i)
		for(int j=0;(i+j)<=n;++j)
			for(int k=0;(i+j+k)<=n;++k){
				int now=id(i,j,k,n-i-j-k);
				val[now]=(fac[n]/fac[i]/fac[j]/fac[k]/fac[n-i-j-k]).getval();
				for(int l=1;l<=m;++l){
					if((!trans[l][0])&&(!trans[l][1])&&(!trans[l][2])&&(!trans[l][3])) continue;
					if((i>=cost[l][0])&&(j>=cost[l][1])&&(k>=cost[l][2])&&(n-i-k-j>=cost[l][3]))
						add(id(i+trans[l][0],j+trans[l][1],k+trans[l][2],n-i-j-k+trans[l][3]),now);
				}
			}
	for(int u=n;u<=(n+1)*(n+1)*(n+1)*n;++u) if(val[u]&&(!dfn[u])) tarjan(u);
	for(int u=n;u<=(n+1)*(n+1)*(n+1)*n;++u) if(val[u])
		for(int i=Head[u];i;i=Next[i]) if(val[V[i]])
			if(belong[V[i]]!=belong[u])
				add(belong[u],belong[V[i]]),++in[belong[V[i]]];
	for(int i=(n+1)*(n+1)*(n+1)*n;i<=num;++i) if(!in[i]) dis[i]=siz[i],Q.push(i);
	while(!Q.empty()){
		int u=Q.front();Q.pop();
		for(int i=Head[u];i;i=Next[i]){
			dis[V[i]]=max(dis[V[i]],dis[u]+siz[V[i]]),--in[V[i]];
			if(!in[V[i]]) Q.push(V[i]);
		}
	}for(int i=(n+1)*(n+1)*(n+1)*n;i<=num;++i) ans=max(ans,dis[i]);
	printf("%lld\n",ans);
}