ARC 077

C
模拟
D
只有2个数是相同的，假设他们的位置是u,v，对于长度k的子序列，重复的情况有$C_{u-1+n-v}^{k-1}$种，减掉

E
考虑把值定为i对n-1次变化花费的减少量之和
对于一次变化l->r，
若l< r，当i∈(l+1,r)，能减少i-(l+1)的花费；
若l>r，当i∈(l+1,m)时，能减少i-(l+1)的花费，
当i∈(1,r)时，能减少m-l+i-1的花费，相当于一个(1-(m-l),r)的l< r的区间

那么把值定为i的时候，他的贡献就是i*包含他的区间个数-区间左端点的和，枚举i找最大贡献

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 210000;

int n,m;
int v[maxn];
struct node{int p,i;}a[maxn<<1]; int tot;
inline bool cmp(const node x,const node y){return x.p<y.p;}
int l[maxn],r[maxn],id;

ll re;

int main()
{
    scanf("%d%d",&n,&m); for(int i=1;i<=n;i++) scanf("%d",&v[i]);
    for(int i=1;i<n;i++)
    {
        if(v[i]<v[i+1])
        {
            re+=v[i+1]-v[i];
            if(v[i]+1!=v[i+1])
            {
                ++id; l[id]=v[i]+1; r[id]=v[i+1]+1;
                a[++tot]=(node){l[id],id};
                a[++tot]=(node){r[id],id};
            }
        }
        else
        {
            re+=m-v[i]+v[i+1];
            if(v[i]%m+1!=v[i+1])
            {
                int kk=m-v[i];
                ++id; l[id]=1-kk; r[id]=v[i+1]+1;
                a[++tot]=(node){l[id],id};
                a[++tot]=(node){r[id],id};
            }
            if(v[i]<m-1)
            {
                ++id; l[id]=v[i]+1,r[id]=m+1;
                a[++tot]=(node){l[id],id};
            }
        }
    }sort(a+1,a+tot+1,cmp);

    int pos=1;
    ll ans=0,sum=0,num=0;
    for(int i=1;i<=m;i++)
    {
        while(pos<=tot&&a[pos].p<=i)
        {
            id=a[pos].i;
            if(a[pos].p!=r[id]) sum+=l[id],num++;
            else sum-=l[id],num--;
            pos++;
        }
        ans=max(ans,num*(ll)i-sum);
    }
    printf("%lld\n",re-ans);

    return 0;
}

F：
翻译一下题解qwq

观察一次变化SS->STST
因为要求T的长度尽量小，所以有|T|<=|S|
因为是在原串的末尾加上若干字符，所以原串不变，SS和STST的一段前缀相同，所以T一定是S的前缀，画个图yy一下能发现T是S的最小周期(S=TTT..T或S=TTTT..TA(A是T的前缀))
因为左右两部分相同，我们只讨论左半部分的变化，用g(S)代表S的变化，分两类情况讨论：
S=TT..T即|T|整除|S|时，易知变化g(S)=ST,g(ST)=STT,g(STT)=STTT..，每次在末尾加上一个T
当S=TT..TA时（A是T的前缀），g(S)=ST=TT..TAT，发现ST的最小正周期是S,g(ST)=STS，重复这个过程得g^(i+2)=g^(i+1)+g^(i)，很像斐波那契数列，到80~90项左右长度就超过10^18，dp一下每类字符个数，找解的时候递归下去找

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
#define inf 1e18
using namespace std;

const int maxn = 210000;

ll L,R;
int n,pos;
bool flag;
int ex[maxn],id,mx;
char str[maxn];
struct node
{
    ll n,s[26];
    node(){n=0;memset(s,0,sizeof s);}
}f[maxn]; int ed;
inline node operator +(node &x,node &y)
{
    node re; re.n=x.n+y.n;
    for(int i=0;i<26;i++) re.s[i]=x.s[i]+y.s[i];
    return re;
}

void dp()
{
    f[0].n=pos;
    for(int i=0;i<pos;i++) f[0].s[str[i]-'a']++;
    f[1].n=n;
    for(int i=0;i<n;i++) f[1].s[str[i]-'a']++;

    ed=1;
    while(f[ed].n<inf) 
        f[ed+1]=f[ed]+f[ed-1],ed++;
}
node solve1(ll x)
{
    ll div=x/f[0].n,oth=x%f[0].n;
    node re; for(int i=0;i<26;i++) re.s[i]=f[0].s[i]*div;
    if(oth)
    {
        for(int i=0;i<oth;i++) re.s[str[i]-'a']++;
    }
    return re;
}
node solve2(ll x,int now)
{
    while(now&&f[now-1].n>=x) now--;
    node re;
    if(now<=1)
    {
        for(int i=0;i<x;i++) re.s[str[i]-'a']++;
        return re;
    }
    re=f[now-1];
    node temp=solve2(x-f[now-1].n,now-2);
    re=re+temp;
    return re;
}
node calc(ll x)
{
    node re;
    if(!x) return re;
    if(flag) re=solve1(x);
    else re=solve2(x,ed);
}

int main()
{
    scanf("%s",str); n=strlen(str); n>>=1; str[n]='$';
    id=1,mx=0;
    for(int i=1;i<n;i++)
    {
        if(mx>=i) ex[i]=min(ex[i-id],mx-i);
        while(str[i+ex[i]]==str[ex[i]]) ex[i]++;
        if(i+ex[i]-1>mx) id=i,mx=i+ex[i]-1;
    }

    flag=false;
    for(pos=1;pos<n&&pos+ex[pos]-1!=n-1;pos++);
    if(n%pos==0) flag=true;
    dp();

    scanf("%lld%lld",&L,&R);
    node a1=calc(R),a2=calc(L-1);
    for(int i=0;i<26;i++) printf("%lld ",a1.s[i]-a2.s[i]);

    return 0;
}