【CSP-S 2019模拟】题解

T1：

傻逼题
显然对于一个$i$，最大的为$\left(\genfrac{}{}{0ex}{}{i}{\frac{i}{2}}\right)$
用堆维护每一行的最大值
先用对数判大小再在取模的情况下计算答案

#include<bits/stdc++.h>
using namespace std;
const int RLEN=1<<20|1;
inline char gc(){
    static char ibuf[RLEN],*ib,*ob;
    (ob==ib)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
    return (ob==ib)?EOF:*ib++;
}
#define gc getchar
inline int read(){
    char ch=gc();
    int res=0,f=1;
    while(!isdigit(ch))f^=ch=='-',ch=gc();
    while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
    return f?res:-res;
}
#define ll long long
#define re register
#define pii pair<int,int>
#define fi first
#define se second
#define pb push_back
#define cs const
#define bg begin
#define poly vector<int>
cs int mod=1e9+7;
inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
inline int mul(int a,int b){return 1ll*a*b%mod;}
inline void Add(int &a,int b){(a+=b)>=mod?a-=mod:0;}
inline void Dec(int &a,int b){(a-=b)<0?a+=mod:0;}
inline void Mul(int &a,int b){a=1ll*a*b%mod;}
inline int ksm(int a,int b,int res=1){for(;b;b>>=1,Mul(a,a))(b&1)&&(Mul(res,a),1);return res;}
inline int Inv(int x){return ksm(x,mod-2);}
inline void chemx(int &a,int b){a<b?a=b:0;}
inline void chemn(int &a,int b){a>b?a=b:0;}
cs int N=1000005;
pii ans[N];
int tot;
namespace Comb{
	int fac[N],ifac[N];
	inline void init(cs int len=N-5){
		ifac[0]=fac[0]=1;
		for(int i=1;i<=len;i++)fac[i]=mul(fac[i-1],i);
		ifac[len]=Inv(fac[len]);
		for(int i=len-1;i;i--)ifac[i]=mul(ifac[i+1],i+1);
	}
	inline int C(int n,int m){
		return n<m?0:mul(fac[n],mul(ifac[m],ifac[n-m]));
	}
	inline int calc(){
		int res=0;
		for(int i=1;i<=tot;i++)
		Add(res,C(ans[i].fi,ans[i].se));
		return res;
	}
}
#define ldb long double
struct node{
	ldb val;int pos,tim;
	node(ldb _v=0,int _p=0,int _t=0):val(_v),pos(_p),tim(_t){}
	friend inline bool operator <(cs node &a,cs node &b){
		return (a.val==b.val)?(a.pos==b.pos?a.tim<b.tim:a.pos<b.pos):a.val<b.val;
	}
};
priority_queue<node> q;
int vis[N],n,k;
ldb v[N],fac[N];
inline void init(cs int len=N-5){
	for(int i=1;i<=len;i++)v[i]=log2(i);
	for(int i=1;i<=len;i++)fac[i]=fac[i-1]+v[i];
}
inline ldb C(int n,int m){
	return n<m?0:(fac[n]-fac[m]-fac[n-m]);
}
int main(){
	Comb::init(),init();
	n=read(),k=read();
	for(int i=0;i<=n;i++){
		if(i&1)q.push(node(C(i,i/2),i,2));
		else q.push(node(C(i,i/2),i,1));
		vis[i]=i/2;
	}
	while(k>0){
		node now=q.top();q.pop();
		int t=min(k,now.tim);k-=t;
		while(t--)ans[++tot]=pii(now.pos,vis[now.pos]);
		vis[now.pos]--;
		if(vis[now.pos]>=0)q.push(node(C(now.pos,vis[now.pos]),now.pos,2));
	}
	cout<<Comb::calc();
}

T2：

显然暴力是每次建出$Trie$统计答案

考虑这个不好维护，用莫队做
先把整体的$Trie$建出来
每次把一个串的权值加到$Trie$上
由于对串分块复杂度会和串长有关
考虑把所有串拼起来分块

这样暴力删每个字符
复杂度就是对的了，每次判断$Trie$上每个节点是否被加入/删除
用一个链表维护每个空段

但是链表只能直接删除和回退
所以要写回滚莫队

复杂度$O(n\sqrt{n})$

#include<bits/stdc++.h>
using namespace std;
const int RLEN=1<<20|1;
inline char gc(){
    static char ibuf[RLEN],*ib,*ob;
    (ob==ib)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
    return (ob==ib)?EOF:*ib++;
}
#define gc getchar
inline int read(){
    char ch=gc();
    int res=0,f=1;
    while(!isdigit(ch))f^=ch=='-',ch=gc();
    while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
    return f?res:-res;
}
#define ll long long
#define re register
#define pii pair<int,int>
#define fi first
#define se second
#define pb push_back
#define cs const
#define poly vector<int>
inline void chemx(int &a,int b){a<b?a=b:0;}
cs int N=300005;
int A,B,C,n,m,a[N],num;
int len[N],pos[N],vis[N];
int v[N],g[N],bel[N],L,blo,cnt,fir[N],bg[N],ed[N];
char s[N]; 
ll ans,anc[N],V;
inline ll S(int x){
	return 1ll*x*(x+1)/2;
}
inline ll gcd(ll a,ll b){
	return b?gcd(b,a%b):a;
}
struct opt{
	int l,r,id;
	friend inline bool operator <(cs opt &a,cs opt &b){
		return bel[a.l]==bel[b.l]?a.r>b.r:a.l<b.l;
	}
}p[N];
int pre[N],suf[N],top;
struct node{
	int l,p,r;
	node(int _l=0,int _p=0,int _r=0):l(_l),p(_p),r(_r){}
}stk[N];
inline void getback(){
	while(top){
		int l=stk[top].l,p=stk[top].p,r=stk[top].r;
		ans-=S(r-l-1),ans+=S(p-l-1),ans+=S(r-p-1);
		pre[r]=p,suf[l]=p,top--;
	}
}
inline void delet(int p){
	int l=pre[p],r=suf[p];
	ans-=S(p-l-1)+S(r-p-1),ans+=S(r-l-1);
	stk[++top]=node(l,p,r);
	pre[r]=l,suf[l]=r;
}
inline void build(){
	for(int i=1,last=0;i<=L+1;i++){
		if(vis[i])suf[last]=i,pre[i]=last,ans+=S(i-last-1),last=i;
	}
}
int dep[N],val[N],nxt[N][26],tot;
inline void insert(int id,int coef){
	int p=0;coef*=v[id];
	for(int i=bg[id];i<=ed[id];i++){
		int c=a[i];
		if(!nxt[p][c])nxt[p][c]=++tot,dep[tot]=dep[p]+1;
		p=nxt[p][c];
		int pre=(val[p]>0&&1ll*B*val[p]+1ll*A*dep[p]>=C);
		val[p]+=coef;
		int nxt=(val[p]>0&&1ll*B*val[p]+1ll*A*dep[p]>=C);
		if(!pre&&nxt){vis[g[dep[p]]]++;}
		if(pre&&!nxt){vis[g[dep[p]]]--;if(!vis[g[dep[p]]])delet(g[dep[p]]);}
	}
}
signed main(){
	n=read(),A=read(),B=read(),C=read();
	for(int i=1;i<=n;i++)v[i]=read();
	for(int i=1;i<=n;i++){
		scanf("%s",s+1);
		len[i]=strlen(s+1);
		bg[i]=num+1;
		for(int j=1;j<=len[i];j++)
		a[++num]=s[j]-'a';
		ed[i]=num;
		chemx(L,len[i]);
	}
	for(int i=1;i<=L;i++)g[i]=read();
	m=read();vis[L+1]=1;
	for(int i=1;i<=m;i++)p[i].l=read(),p[i].r=read(),p[i].id=i;
	blo=sqrt(num);
	bel[1]=1;V=S(L);
	for(int i=1,j=2;j<=n;j++){
		bel[j]=bel[j-1];
		if(ed[j]-bg[i]+1>blo)
		bel[j]++,i=j;
	}
	cnt=bel[n];
	sort(p+1,p+m+1);
	for(int i=1;i<=n;i++)
		if(!fir[bel[i]])fir[bel[i]]=i;
	for(int i=1,j,pos=1;i<=cnt;i++){
		ans=0;
		for(j=fir[i];j<=n;j++)insert(j,1);
		build();
		for(j=n;pos<=m&&bel[p[pos].l]==i;pos++){
			while(j>p[pos].r)insert(j,-1),j--;
			top=0;
			for(int k=fir[i];k<p[pos].l;k++)insert(k,-1);
			anc[p[pos].id]=V-ans;
			for(int k=fir[i];k<p[pos].l;k++)insert(k,1);
			getback();
		}
		while(j>=fir[i])insert(j,-1),j--;
	}
	for(int i=1;i<=m;i++){
		ll gg=gcd(V,anc[i]);
		cout<<anc[i]/gg<<"/"<<V/gg<<'\n';
	}
}

T3：

据说出题人用屁眼造数据
造一个半径为小数点后不知道多少位的圆，平方一下就没了的那种

于是咕咕咕