北京集训DAY6

/*
    sort排序 漏了一个条件只有70 我真是 zz 
    对于一个合法的序列 拿走一个数插入数列中 
    造成的影响就是一些数向前移动 还有数向后移动 
    不管怎么动  向后移动或向前移动的数一定只有一个 
    后者都不动 
*/
#include <algorithm>
#include <cctype>
#include <cstdio>
#define N 10005000
int n, cnt1, cnt2, sum;
struct node
{
    int num, pos;
    bool operator < (node a)const
    {
        if(num==a.num) return pos<a.pos;
        return num<a.num;
    }
}a[N];
int main(int argc, char *argv[])
{
    //freopen("sort.in", "r", stdin);
    //freopen("sort.out", "w", stdout);
    scanf("%d", &n);
    for (int i = 1; i <= n; ++i) scanf("%d", &a[i].num), a[i].pos=i;
    std :: sort(a + 1, a + n + 1);
    for (int i = 1; i <= n; ++i)
    {
        if (a[i].pos > i) cnt1++;
        if (a[i].pos < i) cnt2++;
    }
    if (cnt1 == 1 || cnt2 == 1 || (!cnt2 && !cnt1)) puts("YES");
    else puts("NO");
    return 0;
    //fclose(stdin); fclose(stdout);
}

/*
    一个模板吧 。
    对于60%的数据ai互质 显然可以用CRT来求解
    但是后40%ai不互素 CRT求解显然不对 
    这用了另一种方法 合并多项式 
    即把4个式子利用exgcd合并成1个式子
    即可求得ans 
*/
#include<iostream>
#include<cstdio>
typedef long long LL; 
using namespace std;

LL n,m[5],a[5];
LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if (b == 0)
    {
        x = 1, y = 0;
        return a;
    }
    LL r = exgcd(b, a % b, x, y);
    LL tmp = x;
    x = y;
    y = tmp - a / b * y;
    return r;
}
inline LL crt()
{
    LL a1 = a[1], a2, m2, d, c, m1=m[1];
    for (LL i = 2; i <= 4; ++i)
    {
        a2 = a[i], m2 = m[i];
        c = a2 - a1;
        LL x, y;
        d = exgcd(m1, m2, x, y);
        x = x * c / d;
        int mod = m2 / d;
        x = (mod + x % mod) % mod;
        a1 += m1 * x;
        m1 *= mod;
    }
    return a1;
}
int main(int argc, char *argv[])
{
    freopen("mod.in","r",stdin);
    freopen("mod.out","w",stdout);
    for(int i=1;i<=4;i++) cin>>m[i]>>a[i];
    cout<<crt()<<endl;
    return 0;
    fclose(stdin); fclose(stdout);
}

/*
    我们发现回文串是可以二分的
    比如：bbbabcbaabcba 我们以第1个c为
    中心，发现回文半径是3，大于3一定不是回文串。当回文串长度为偶数的时候
    我们需要特殊处理一下。
    现在有一个结论：本质不同的回文串个数只有O(N)个
    本质不同：字符串本身是不同的。
    每一次处理完回文串，我们要把他的hash值记录下来。
*/
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <map>
using namespace std;

typedef unsigned long long ULL;
typedef long long LL;

char s[10005];
ULL h[10005],rh[10005],pw[10005];
int L;

ULL hs(int l,int r){
    return h[r]-h[l-1]*pw[r-l+1];
}

ULL rhs(int l,int r){
    return rh[l] - rh[r+1]*pw[r-l+1];
}

struct N{
    int a[26];
    bool ok(){
        int b[26];
        for(int i=0;i<26;i++) b[i]=a[i];
        sort(b,b+26);
        for(int i=0;i<25;i++){
            if(b[i]>0&& b[i] == b[i+1]) return true;
        }
        return false;
    }
    void clear(){
        memset(a,0,sizeof a);
    }
};

LL ans=0;
map<ULL,LL> num;
map<ULL,N> A;
void solve_odd(){
    for(int i=1;i<=L;i++){
        int l = 1,r = min(i,L-i+1)+1;
        while(r-l>1){
            int mid = (l+r)/2; 
            if(hs(i-mid+1,i+mid-1)== rhs(i-mid+1,i+mid-1)) l=mid;
            else r=mid;
        }
        int p=l;
        int tmp = p;
        while(tmp>=1&&num.find(hs(i-tmp+1,i+tmp-1))==num.end()) tmp--;
        LL sum = 0;
        N st;
        st.clear();
        if(tmp>=1){
            sum=num[hs(i-tmp+1,i+tmp-1)];
            st = A[hs(i-tmp+1,i+tmp-1)];
        }
        while(tmp<p){
            st.a[s[i+tmp]-'a']+= (tmp == 0?1:2);
            if(st.ok()) sum++;
            num[hs(i-tmp,i+tmp)] = sum;
            A[hs(i-tmp,i+tmp)] = st;
            tmp++;
        }
        ans+=sum;  
    }
}

void solve_even() {
    A.clear();
    num.clear();
    for(int i=1;i<L;i++){
        int l = 1,r = min(i,L-i)+1;
        while(r-l>1){
            int mid = (l+r)/2;
            if(hs(i-mid+1,i+mid)== rhs(i-mid+1,i+mid)) l=mid;
            else r=mid;
        }
        int p=l;
        int tmp = p;
        while(tmp>=1&&num.find(hs(i-tmp+1,i+tmp))==num.end()) tmp--;
        LL sum = 0;
        N st;
        st.clear();
        if(tmp>=1){
            sum=num[hs(i-tmp+1,i+tmp)];
            st = A[hs(i-tmp+1,i+tmp)];
        }
        while(tmp<p){
            st.a[s[i+tmp+1]-'a']+= 2;
            if(st.ok()) sum++;
            num[hs(i-tmp,i+tmp+1)] = sum;
            A[hs(i-tmp,i+tmp+1)] = st;
            tmp++;
        }
        ans+=sum;
    }
}

int main(){
    freopen("str.in","r",stdin);
    freopen("str.out","w",stdout);
    scanf("%s",s+1);
    L=strlen(s+1);
    s[0]='#';
    pw[0]=1;
    for(int i=1;i<=L;i++) pw[i] = pw[i-1]*13131 ;
    for(int i=1;i<=L;i++) h[i] = h[i-1]*13131 + s[i];
    for(int i=L;i>=1;i--) rh[i] = rh[i+1]*13131 + s[i];    
    solve_odd();
    solve_even();
    printf("%lld
",ans);
    fclose(stdout);
    return 0;
}