hdu 5651 xiaoxin juju needs help 逆元 两种求解方式

xiaoxin juju needs help

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 1159    Accepted Submission(s): 335

Problem Description

As we all known, xiaoxin is a brilliant coder. He knew **palindromic** strings when he was only a six grade student at elementry school.

This summer he was working at Tencent as an intern. One day his leader came to ask xiaoxin for help. His leader gave him a string and he wanted xiaoxin to generate palindromic strings for him. Once xiaoxin generates a different palindromic string, his leader will give him a watermelon candy. The problem is how many candies xiaoxin's leader needs to buy?

 

Input

This problem has multi test cases. First line contains a single integer T(T≤20) which represents the number of test cases.
For each test case, there is a single line containing a string S(1≤length(S)≤1,000).

 

Output

For each test case, print an integer which is the number of watermelon candies xiaoxin's leader needs to buy after mod 1,000,000,007.

 

Sample Input

3
aa
aabb
a

 

Sample Output

1
2
1

 

Source

BestCoder Round #77 (div.2)

题意：给你一个字符串，判断有多少种方式使得这个字符串回文；

思路：先判断0的情况，其次标记26个字母，ans=c(len/2,a/2)*c(len-a/2,b/2)*c(len/2-a/2-b/2,c/2)......;

　　   因为在求组合数需要取模，所以需要利用逆元的方式求解；

　　   逆元详解：http://blog.csdn.net/acdreamers/article/details/8220787

　　　扩张欧几里德求解逆元

　　　第二个是用费马小定理求逆元

#include<iostream>
#include<cstdio>
#include<cmath>
#include<string>
#include<queue>
#include<algorithm>
#include<stack>
#include<cstring>
#include<vector>
#include<list>
#include<set>
#include<map>
using namespace std;
#define ll __int64
#define mod 1000000007
int scan()
{
    int res =  , ch ;
    while( !( ( ch = getchar() ) >= '' && ch <= '' ) )
    {
        if( ch == EOF )  return  <<  ;
    }
    res = ch - '' ;
    while( ( ch = getchar() ) >= '' && ch <= '' )
        res = res *  + ( ch - '' ) ;
    return res ;
}
char a[];
ll flag[];
void extend_Euclid(ll a, ll b, ll &x, ll &y)
{
    if(b == )
    {
        x = ;
        y = ;
        return;
    }
    extend_Euclid(b, a % b, x, y);
    ll tmp = x;
    x = y;
    y = tmp - (a / b) * y;
}
ll combine1(ll n,ll m) //计算组合数C(n,m)
{
    ll sum=; //线性计算
    for(ll i=,j=n;i<=m;i++,j--)
    {
        sum*=j;
        sum%=mod;
        ll x,y;
        extend_Euclid(i,mod,x,y);
        sum*=(x%mod+mod)%mod;
        sum%=mod;
    }
    return sum;
}
int main()
{
    ll x,y,z,i,t;
    scanf("%I64d",&z);
    while(z--)
    {
        memset(flag,,sizeof(flag));
        scanf("%s",a);
        x=strlen(a);
        for(i=;i<x;i++)
        flag[a[i]-'a']++;
        ll sum=;
        for(i=;i<;i++)
        {
            if(flag[i]%)
            sum++;
        }
        if(x%==&&sum)
        printf("0\n");
        else if(x%==&&sum>)
        printf("0\n");
        else
        {
            ll f=x/;
            ll ans=;
            for(i=;i<;i++)
            {
                if(flag[i]/)
                {
                    ans*=combine1(f,flag[i]/);
                    f-=flag[i]/;
                    ans%=;
                }
            }
            printf("%I64d\n",ans);
        }
    }
    return ;
}

#include<iostream>
#include<cstdio>
#include<cmath>
#include<string>
#include<queue>
#include<algorithm>
#include<stack>
#include<cstring>
#include<vector>
#include<list>
#include<set>
#include<map>
using namespace std;
#define ll __int64
#define mod 1000000007
int scan()
{
    int res =  , ch ;
    while( !( ( ch = getchar() ) >= '' && ch <= '' ) )
    {
        if( ch == EOF )  return  <<  ;
    }
    res = ch - '' ;
    while( ( ch = getchar() ) >= '' && ch <= '' )
        res = res *  + ( ch - '' ) ;
    return res ;
}
char a[];
ll flag[];
ll poww(ll a,ll n)//快速幂
{
   ll r=,p=a;
   while(n)
   {
       if(n&) r=(r*p)%mod;
       n>>=;
       p=(p*p)%mod;
   }
   return r;
}
ll combine1(ll n,ll m) //计算组合数C(n,m)
{
    ll sum=; //线性计算
    for(ll i=,j=n;i<=m;i++,j--)
    {
        sum*=j;
        sum%=mod;
        sum*=poww(i,);
        sum%=mod;
    }
    return sum;
}
int main()
{
    ll x,y,z,i,t;
    scanf("%I64d",&z);
    while(z--)
    {
        memset(flag,,sizeof(flag));
        scanf("%s",a);
        x=strlen(a);
        for(i=;i<x;i++)
        flag[a[i]-'a']++;
        ll sum=;
        for(i=;i<;i++)
        {
            if(flag[i]%)
            sum++;
        }
        if(x%==&&sum)
        printf("0\n");
        else if(x%==&&sum>)
        printf("0\n");
        else
        {
            ll f=x/;
            ll ans=;
            for(i=;i<;i++)
            {
                if(flag[i]/)
                {
                    ans*=combine1(f,flag[i]/);
                    f-=flag[i]/;
                    ans%=;
                }
            }
            printf("%I64d\n",ans);
        }
    }
    return ;
}