C
模拟
D
只有2个数是相同的,假设他们的位置是u,v,对于长度k的子序列,重复的情况有
Ck−1u−1+n−v
种,减掉
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;
}