大意:求出一个串不通过增删组成多少个回文串。
思路:有公式S (ni) !/((n1)!*(n2)!….(ni)! )=x;x就是可重复序列的最终组成不同的序列数。那么既然还有取余的问题,所以用到逆元(就是化除为乘),那么就可以用快速幂处理(参数二为mod-2).
#include<map>
#include<queue>
#include<cmath>
#include<cstdio>
#include<stack>
#include<iostream>
#include<cstring>
#include<algorithm>
#define inf 0x3f3f3f3f
#define eps 1e-8
#define mod 1000000007
#define ls l,mid,rt<<1
#define rs mid+1,rt,rt<<1|1
#define LL __int64
using namespace std;
char s[1010];
LL ha[30],f[1010]={1,1};
LL so(LL n ,LL m){
LL tmp = 1;
while(m){
if(m&1)
tmp =( tmp * n )%mod;
n = (n*n) % mod;
m >>= 1;
}
return tmp;
}
int main(){
LL n,m,i,j,k;
for(i = 2 ;i <= 1010 ;++ i){
f[i] = (f[i-1]*i) % mod;
}
scanf("%I64d",&n);
while(n--){
memset(ha,0,sizeof(ha));
scanf("%s",s);
int le = strlen(s);
if(le == 1){
puts("1");continue;
}
LL x = 0;
for(i = 0;i < le;++ i){
++ha[s[i]-97];
}
for(i = 0;i < 26; ++ i){
if(ha[i]){
x += (ha[i]/2);//cout<<x<<endl;
}
}
int cnt = 0;
bool vis = false;
for(i = 0;i <= 25 ;++ i){
if(ha[i]&1){
cnt++;
}
if(cnt >= 2){
vis = true;
break;
}
}
if(vis){
puts("0");continue;
}
LL tmp = f[x];
//cout<<tmp<<endl;
for(i = 0 ;i < 26 ;++ i)
if( ha[i] )
tmp = (tmp * so(f[ha[i]/2],mod-2) )% mod;
printf("%I64d\n",tmp);
}
return 0;
}