G - Maximum Palindromes

G - Maximum Palindromes
题目链接

题意：给你一个字符串，然后让你找出最长的回文串有多少个。一道全排列题目。

思路：首先把[a - z]有多少个先求出来，然后对[a - z]的个数进行分类讨论，处理出来有多少个2a, 2b,```,2z;然后根据这个和剩下的a, b,````,z算数就可以了。

求除法用乘法逆元:
( a / b ) % mod =a * inv ( b , mod ) %mod =( a%mod * inv ( b , mod )%mod ) %mod

#include<bits/stdc++.h>
#define int unsigned long long int
using namespace std;

const int N = 1e5 + 10;
int mod = 1e9 + 7;//模
char st[N];//存字符串
int q, l, r;
int sum[N][30];//求字母个数，前缀和
int num[N];//处理出来1 * 2 * ··· * N - 1

void startt(){
    num[0] = 1;
    for(int i = 1; i < N; i++){
        num[i] = (num[i - 1] * i)%mod;
    }
}

int exgcd(int a, int b, int &x, int &y){
    if(!b){
        x = 1, y = 0;return a;
    }
    int res = exgcd(b, a%b, x, y);
    int tmp = x;
    x = y, y = tmp - (a / b) * y;
    return res;
}

int exgcd_inv(int a, int mod){
    int d, x, y;
    d = exgcd(a, mod, x, y);
    return (x + mod) % mod;
}//求逆元

signed main(){
    ios::sync_with_stdio(false);
    cin.tie(0); cout.tie(0);
    startt();
    cin >> st + 1;
    int len = strlen(st + 1);
    for(int i = 1; i <= len; i++){
        for(int j = 0; j < 26; j++){
            if(st[i] == (j + 'a')){
                sum[i][j] = sum[i - 1][j] + 1;
            }
            else sum[i][j] = sum[i - 1][j];
        }
    }
    cin >> q;
    while(q -- ){
        cin >> l >> r;
        int s[30];
        memset(s, 0, sizeof s);
        //求出来a ~ z每个字母的个数
        for(int i = 0; i < 26; i++){
            s[i] = sum[r][i] - sum[l - 1][i];
        }
        int tot = 0, numodd = 0, numeven = 0;
        for(int i = 0; i < 26; i++){
            if(s[i] % 2 == 0){
                numeven += s[i];
                tot += s[i];
                s[i]/= 2;
            }
            else{
                numeven += s[i] - 1;
                numodd += 1;
                tot += s[i];
                s[i]--;
                s[i]/=2;
            }
        }
        if(numodd == 0){
            int ans = 1, ans1 = 1;
            int cnt = 0;
            for(int i = 0; i < 26; i++){
                if(s[i] != 0){
                    ans1 = (ans1 * num[s[i]]) % mod;
                }
            }
            ans = num[numeven/2];
            ans = (ans % mod * exgcd_inv(ans1, mod)%mod)%mod;
            cout << ans << endl;
        }
        else{
            if(numeven == 0){
                cout << tot << endl;
            }
            else{
                int ans = 1, ans1 = 1;
                int cnt = 0;
                for(int i = 0; i < 26; i++){
                if(s[i] != 0){
                    ans1 = (ans1 * num[s[i]]) % mod;
                }
            }
            ans = num[numeven/2];
            ans = (ans % mod * exgcd_inv(ans1, mod)%mod)%mod;
                ans = (ans * numodd) %mod;
                cout << ans << endl;
            }
        }
    }
    return 0;
}