HDU5651(逆元)

最新推荐文章于 2024-06-17 00:07:21 发布

西风丶

最新推荐文章于 2024-06-17 00:07:21 发布

阅读量286

点赞数

分类专栏： STL&数据结构数论（学）

本文链接：https://blog.csdn.net/West___wind/article/details/51356316

版权

数论（学）同时被 2 个专栏收录

10 篇文章 0 订阅

订阅专栏

STL&数据结构

8 篇文章 0 订阅

订阅专栏

点击打开链接

(a/b)%c=a*b^(c-2)%c c为质数

#include<cstdio>
#include<map>
#include<cstring>
#define N 1000000007
using namespace std;
typedef long long ll;
ll quickpow(ll a,ll b)
{
    ll r=1,c=a;
    while(b)
    {
        if(b%2)
        r*=c;
        r%=N;
        c*=c;
        c%=N;
        b/=2;
    }
    return r;
}
ll f(int n)
{
    ll result=1;
    for(int i=2;i<=n;i++)
    result=(result*i)%1000000007;
    return result;
}
int main()
{
    ll T,n,i,len,num;
    ll ans,flag;
    char a[1100];
    scanf("%lld",&T);
    while(T--)
    {
        flag=1;
        scanf("%s",a);
        len=strlen(a);
        map<char,ll> b;
        for(i=0;i<len;i++)
        {
            b[a[i]]++;
        }
        map<char,ll>::iterator p;
        num=0;
        n=0;
        for(p=b.begin();p!=b.end();p++)
        {
            n++;
            if((p->second)%2==1)
            num++;
        }
        if(len%2==1&&num!=1)
        ans=0;
        else if(len%2==1&&num==1)
        {
            ans=f(len/2);
            for(p=b.begin();p!=b.end();p++)
            {
                if((p->second)>1)
                {
                    flag*=f((p->second)/2);
                    flag%=N;
                }
            }
        }
        else if(len%2==0&&num!=0)
        ans=0;
        else if(len%2==0&&num==0)
        {
            ans=f(len/2);
            for(p=b.begin();p!=b.end();p++)
            {
                if((p->second)%2==0)
                {
                    flag*=f((p->second)/2);
                    flag%=N;
                }
            }
        }
        printf("%lld\n",(ans*quickpow(flag,N-2))%N);
    }
    return 0;
}