【TC SRM 499】ImpossibleGame（SCC）（拓扑序DP）

传送门

---

又想起了当初在TC上刷网络流的时候被submit failed支配一下午的恐惧。。。

题解：

发现由于我们允许交换任意两个位置的字符，对我们唯一有用的就只剩下字符个数，30的不超过4个的有序非负整数拆分有$5456$种，权值显然就是$\frac{n!}{a!b!c!d!}$，直接组合数算一下即可，连边可以直接暴力连，那么缩一个SCC之后就是拓扑序DP求DAG最长链了。

由于是考场代码改的，并不是我一贯的在TopCoder上做题的代码风格，好吧估计也没人关心我的TopCoder代码风格是什么样子

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

using std::cerr;
using std::cout;

cs int N=31,M=6e3+7;//30的有序正整数拆分只有5456种 

int n;
int id[N][N][N][N],tot;
ll w[M],C[N][N];
std::vector<int> G[M]; 

inline void init(){
	for(int re i=0;i<=30;++i){
		C[i][0]=1;
		for(int re j=1;j<=i;++j)C[i][j]=C[i-1][j-1]+C[i-1][j];
	}
}
inline ll calc(int a,int b,int c,int d){
	ll ans=1;
	ans*=C[n][a];
	ans*=C[n-a][b];
	ans*=C[n-a-b][c];
	ans*=C[n-a-b-c][d];
	return ans;
}

int dfn[M],low[M],dfc;
int scc[M],inst[M],st[M],tp,sct;

void tarjan(int u){
	inst[u]=true;st[++tp]=u;
	dfn[u]=low[u]=++dfc;
	for(int re e=G[u].size()-1;~e;--e){
		int v=G[u][e];
		if(!dfn[v]){
			tarjan(v);
			low[u]=std::min(low[u],low[v]);
		}
		else if(inst[v])low[u]=std::min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u]){
		++sct;
		while(true){
			int v=st[tp--];
			scc[v]=sct;inst[v]=false;
			if(v==u)break;
		}
	}
}

namespace DAG{
	int el[M],nxt[M*100],to[M*100],ect;
	ll w[M],f[M];int d[M],q[M],qn;
	
	inline void adde(int u,int v){
		nxt[++ect]=el[u],el[u]=ect,to[ect]=v,++d[v];
	}
	
	inline void top_sort(){
		for(int re i=1;i<=sct;++i)if(!d[i])q[++qn]=i;
		for(int re i=1;i<=qn;++i){
			int u=q[i];
			for(int re e=el[u],v=to[e];e;v=to[e=nxt[e]])if(!--d[v])q[++qn]=v;
		}
		for(int re i=qn;i;--i){
			int u=q[i];ll t=0;
			for(int re e=el[u],v=to[e];e;v=to[e=nxt[e]])t=std::max(t,f[v]);
			f[u]=w[u]+t;
		}
	}
	
	inline ll solve(){
		for(int re u=1;u<=tot;++u){
			w[scc[u]]+=::w[u];
			for(int re e=G[u].size()-1;~e;--e){
				int v=G[u][e];if(scc[u]!=scc[v])adde(scc[u],scc[v]);
			}
		}
		top_sort();ll ans=0;
		for(int re i=1;i<=sct;++i)ans=std::max(ans,f[i]);
		return ans;
	}
}

using std::string;

class ImpossibleGame{
	public:
		ll getMinimum(int k,std::vector<string> before,std::vector<string> after){
			n=k;init();
			for(int a=0;a<=n;++a)
			for(int b=0;a+b<=n;++b)
			for(int c=0;a+b+c<=n;++c){
				int d=n-a-b-c;
				id[a][b][c][d]=++tot;
				w[tot]=calc(a,b,c,d);
			}
			for(int re i=0;i<before.size();++i){
				auto s=before[i],t=after[i];
				int ca=0,cb=0,cc=0,cd=0;
				int cA=0,cB=0,cC=0,cD=0;
				for(int re i=0;i<s.size();++i){
					switch(s[i]){
						case 'A':++ca;break;
						case 'B':++cb;break;
						case 'C':++cc;break;
						case 'D':++cd;break;
					}
					switch(t[i]){
						case 'A':++cA;break;
						case 'B':++cB;break;
						case 'C':++cC;break;
						case 'D':++cD;break;
					}
				}
				if(ca==cA&&cb==cB&&cc==cC&&cd==cD)continue;
				cA-=ca,cB-=cb,cC-=cc,cD-=cd;
				for(int a=ca;a<=n;++a)
				for(int b=cb;a+b<=n;++b)
				for(int c=cc;a+b+c<=n;++c){
					int d=n-a-b-c;if(d<cd)continue;
					int A=a+cA,B=b+cB,C=c+cC,D=d+cD;
					G[id[a][b][c][d]].push_back(id[A][B][C][D]);
				}
			}
			for(int re i=1;i<=tot;++i)if(!dfn[i])tarjan(i);
			return DAG::solve();
		}
};

#ifdef zxyoi

ImpossibleGame Solver;

signed main(){
	std::cout<<Solver.getMinimum(	
6,
{ "AABBC", "AAAADA", "AAACA", "CABAA", "AAAAAA", "BAAAA" },
{ "AACCB", "DAAABC", "AAAAD", "ABCBA", "AABAAA", "AACAA" }
	)<<"\n";
	return 0;
}

#endif