描述：
给出一个单词，在单词中插入若干字符使其为回文串，求回文串的个数（|s|<=200,n<=10^9）
这道题超神奇，不可多得的一道好题
首先可以搞出一个dp[l][r][i]表示回文串左边i位匹配到第l位，右边i位匹配到第r位的状态数，可以发现可以用矩阵乘法优化（某人说看到n这么大就一定是矩阵乘法了= =）
但这样一共有|s|^2个节点，时间复杂度无法承受
我们先把状态树画出来：例如add
可以发现是个DAG
我们考虑把单独的每条链拿出来求解，那么最多会有|s|条不同的链，链长最多为|s|，时间复杂度为O(|s|^4log n)还是得跪
好像没什么思路了对吧= =（我第一步转化就没想到了= =）
我们考虑记有24个自环的为n24，25个自环的为n25，可以发现n24+n25*2=|s|或|s|+1也就是说对于一个确定的n24，一定有一个确定的n25
那么这样构图：
可以发现所有状况都被包括进来了！！！
那么一共有2|s|个节点，时间复杂度降了一个|s|，看上去好像还是不行
压常数= =
可以发现这个是棵树，也就是说如果按拓扑序编号的话，到时的矩阵左下角将是什么都没有的
那么就直接for i = 1 to n j = i to n k=i to j 就行了 = =
总结下吧
这道题为何神奇呢
首先它把一个DAG的图拆成了若干条相似的链
然后它又把这些链和成了一个更和谐的图
最后再观察题目性质得到一个比较神奇的优化方法
这给了我们什么启迪呢= =
首先遇到某些DAG我们可以考虑拆成若干条相似的链
遇到某些链我们可以考虑把他们合成一个图
最重要的是，还是得参透题目的性质
这道题基本都是依靠题目的性质到达下一步的，只有真正读懂读透这道题，我们才能想出更好的解法
CODE:
#include<cstdio> #include<iostream> #include<cstring> #include<algorithm> using namespace std; #define maxn 410 #define mod 10007 typedef int ll; struct marix{ int r,c;ll a[maxn][maxn]; inline void init(int x){r=c=x;for (int i=1;i<=x;i++) a[i][i]=1;} }x,y; inline void muti(marix &ans,const marix x,const marix y){ ans.r=ans.c=x.r; for (int i=1;i<=x.r;i++) for (int j=i;j<=y.c;j++) { int tmp=0; for (int k=i;k<=j;k++) (tmp+=x.a[i][k]*y.a[k][j])%=mod; ans.a[i][j]=tmp; } } inline void power(marix &ans,marix x,int y) { ans.init(x.r); for (;y;y>>=1) { if (y&1) muti(ans,ans,x); muti(x,x,x); } } ll f[210][210][210]; char s[maxn]; inline ll calc(int l,int r,int x) { ll &u=f[x][l][r]; if (u!=-1) return u; u=0; if (l==r) return u=x==0; if (s[l]==s[r]) { if (l+1==r) return u=x==0; return u=calc(l+1,r-1,x); } if (x>0) return u=(calc(l+1,r,x-1)+calc(l,r-1,x-1))%mod; return u; } int main(){ int n,m; memset(f,-1,sizeof(f)); scanf("%s",s+1); scanf("%d",&n); m=strlen(s+1); n+=m; int l=(n+1)/2,n24=m-1,n25=(m+1)/2,n26=n25; x.r=x.c=n24+n25+n26; for (int i=1;i<=n24;i++) x.a[i][i]=24,x.a[i][i+1]=1; for (int i=n24+1;i<=n25+n24;i++) x.a[i][i]=25,x.a[i][i+n25]=1; for (int i=n24+n25+1;i<=n25+n24+n26;i++) x.a[i][i]=26; for (int i=n24+1;i<n25+n24;i++) x.a[i][i+1]=1; marix y; power(y,x,l-1); muti(x,y,x); ll ans; for (int i=0;i<=n24;i++) { int j=(m-i+1)/2,k=l-i-j; if (k<0) continue; ll sum=calc(1,m,i); (ans+=sum*x.a[n24-i+1][n24+j+n25]%mod)%=mod; if ((n&1)&&(m-i&1^1)) (ans=ans-sum*y.a[n24-i+1][n24+j]%mod+mod)%=mod; } printf("%d ",ans); return 0; }