hdu 5651

问题描述
xiaoxin巨从小就喜欢字符串，六年级的时候他就知道了什么是回文串。这时，xiaoxin巨说到：如果一个字符串 SS 是回文串，那么该字符串从前往后看和从后往前看是一样一样的。


六年级的暑假，xiaoxin很快就做完了暑假作业，然后到腾讯做起了实习生。这日，leader给了xiaoxin一个字符串，请xiaoxin帮忙写一个函数来生成所有可能的回文串，可以任意改变字符串的顺序但是不可以扔掉某个字符。并且leader告诉xiaoxin，每生成一个不一样的回文串就可以得到一颗西瓜糖。


请你帮忙计算xiaoxin的leader最多需要买多少颗西瓜糖呢？
输入描述
多组测试数据。第一行包含一个整数 T(T\leq 20)T(T≤20) 表示组数。每组测试数据给出一个只包含小写字母的字符串 S(1\leq length(S)\leq 1,000)S(1≤length(S)≤1,000)
输出描述
对于每组测试数据，输出一个数, 表示leader需要买的西瓜糖的个数，结果对 1,000,000,0071,000,000,007 取模。
输入样例
3
aa
aabb
a
输出样例
1
2
1

首先，如果不止一个字符出现的次数为奇数，则结果为0。 否则，我们把每个字符出现次数除2，也就是考虑一半的情况。 那么结果就是这个可重复集合的排列数了。 设该集合有n个数，第i个数出现次数为ai，那么结果就是 fact(n)/fact(a_1)/fact(a_2)/..../fact(a_n)fact(n)/fact(a1)/fact(a2)/..../fact(an)。

#include <cstdio>
#include <cstring>
using namespace std;
char str[1010];
int fig[200];
long long A[1010];
const long long mod = 1e9 + 7;

long long exGcd(long long a, long long b, long long &x, long long &y)
{
    if (b == 0)
    {
        x = 1;
        y = 0;
        return a;
    }
    long long res = exGcd(b, a % b, y, x);
    y -= a / b * x;
    return res;
}

long long Inv(long long a, long long m)
{
    long long d, x, y;
    d = exGcd (a, m, x, y);
    if(d == 1) return (x % m + m) % m;
    return -1;
}

int main()
{
    int T;
    scanf("%d", &T);
    A[1] = 1;
    for (int i = 2; i <= 1005; i++) A[i] = A[i - 1] * i % mod;
    while (T--)
    {
        scanf("%s", str);
        int len = strlen(str);
        if (len == 1)
        {
            printf("1\n");
            continue;
        }
        memset(fig, 0, sizeof(fig));
        for (int i = 0; i < len; i++)
            fig[str[i]] ++;
        int num = 0;
        for (int i = 0; i < 150; i++)
            if (fig[i] % 2 == 1)
                num++;
        if ((num > 1) || (num == 1 && len % 2 == 0))
        {
            printf("0\n");
            continue;
        }
        long long ans = 1;
        int left;
        if (num == 1) left = (len - 1) / 2;
        else left = len / 2;
        for (int i = 0; i < 150; i++)
        {
            if (fig[i] && fig[i] > 1)
            {
                ans = ans * A[fig[i] / 2] %  mod;
                //printf("(%d)\n", A[fig[i] / 2] );
            }
        }
        ans = A[left] * Inv(ans, mod) % mod;
        printf("%I64d\n", ans);
    }
    return 0;
}