HDU 5651 xiaoxin juju needs help 逆元&组合数学

大意:求出一个串不通过增删组成多少个回文串。

思路:有公式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;
}

你可能感兴趣的:(组合,逆元)