SGU201 Non Absorbing DFA

题目大意:给定一个自动机，给出所有转移。转移分两种，一种是直接转移，字符c被接受，另一种是转移但字符不会被接受，即转移到的状态仍以字符c为转移。求给定长度的有多少个不同的字符串转被自动机接受。

若是普通的自动机转移，则dp求解，即以长度为阶段，自动机的所有状态作为状态转移即可。

但由于多了一种转移，只需要初始化即可，即确定每个转移最终会转移到的状态，注意会出现无限循环，判定为不可转移。

结果很大需要高精度。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int mod=100000000;
const int length=101;
char ch[101];
int n,m;
int len;
int S;
bool vst[1001];
int a[1001][26];
bool b[1001][26];
int d[1001][61][length];
bool bj[1001];
void dfs(int u,int v) {
	int j;
	bj[u]=1;
	j=a[u][v];
	if (b[j][v]&&!bj[j]) dfs(j,v);
	if (!b[j][v]) a[u][v]=a[j][v],b[u][v]=0;
}
void init() {
	int i,j;
	for (j=0; j<m; j++) {
		memset(bj,0,sizeof bj);
		for (i=1; i<=n; i++)
			if (b[i][j]&&!bj[i])
				dfs(i,j);
	}
}
void add(int a[length],int b[length]) {
	int i;
	if (a[0]>b[0]) b[0]=a[0];
	for (i=1; i<=b[0]; i++) b[i]+=a[i];
	for (i=1; i<=b[0]; i++) {
		b[i+1]+=b[i]/mod;
		b[i]%=mod;
	}
	if (b[b[0]+1]) b[0]++;
}
int ans[length];
void print(int a[length]) {
	int i;
	if (a[0]==0) {
		printf("0\n");
		return;
	}
	printf("%d",a[a[0]]);
	for (i=a[0]-1; i; i--) {
		printf("%08d",a[i]);
	}
	printf("\n");
}
void dp() {
	int i,j,k;
	d[S][0][0]=d[S][0][1]=1;
	for (k=0; k<=len; k++) {
		for (i=1; i<=n; i++) {
			for (j=0; j<m; j++)
				if (!b[i][j]) add(d[i][k],d[a[i][j]][k+1]);
		}
	}
	for (i=1; i<=n; i++)
		if (vst[i]) add(d[i][len],ans);
	print(ans);
}
int main() {
	int i,j,t,x;
	scanf("%s",ch);
	m=strlen(ch);
	scanf("%d",&n);
	scanf("%d",&S);
	scanf("%d",&t);
	while (t--) {
		scanf("%d",&x);
		vst[x]=1;
	}
	for (i=1; i<=n; i++)
		for (j=0; j<m; j++) {
			scanf("%d",&x);
			a[i][j]=x;
		}
	for (i=1; i<=n; i++)
		for (j=0; j<m; j++) {
			scanf("%d",&x);
			b[i][j]=x;
		}
	scanf("%d",&len);
	init();
	dp();
	return 0;
}