2018.7.31 Noip2018模拟测试赛（十六）

日期：	七月最后一天	总分：	300分
难度：	提高 ~ 省选	得分：	30分（少的可怜）

我太弱了：（题目目录）

T1：Mushroom追妹纸

T2：抵制克苏恩

T3：美味

失分分析：（QAQ）

开始全部题目看了一遍，第二题期望dp，果断放弃……

看到T3，感觉像是线性基，但是要修改，似乎不可做……（QAQ）

只剩下T1，想到正解——后缀数组+KMP，结果忘记怎么打，耗了整个比赛……

最后T1打炸了，成功炸成屎……

剩下的二十分钟，赶紧打一个T3线性基暴力，但样例竟然过不了……

定睛一看，[○･｀Д´･ ○]，我TM竟然理解错题意了！！MMPPP！！

最后T3暴力30分………………

题解：

T1：Mushroom追妹纸

字符串题，前两个条件后缀数组求最长公共字串。

第三个条件，用 KMP 求出$s_3$在前两个串出现的位置，在统计答案的时候，不要选到$s_3$就好了

CODE:

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 using namespace std;

 int n,n1,n2,n3,m,ans=;
 int nxt[],pos[],f[];
 int tp[],rak[],tax[];
 char s[],s1[],s2[],s3[];
 int height[],sa[];

 void get_next(){
     int i=,j=;
     for(int i=;i<=n3;i++)nxt[i]=;
     nxt[]=;
     while(i<=n3){
         if(j==||s3[i]==s3[j])
             if(s3[++i]!=s3[++j])nxt[i]=j;
             else nxt[i]=nxt[j];
         else j=nxt[j];
     }
 }

 void solve(){
     int i=,j=;
     while(i<=n){
         if(j==||s[i]==s3[j])i++,j++;
         else j=nxt[j];
         if(j>n3){
             pos[++pos[]]=i-;
             j=nxt[j];
         }
     }
 }

 void sort(int a[],int b[]){
     for(int i=;i<=m;i++)tax[i]=;
     for(int i=;i<=n;i++)tax[a[i]]++;
     for(int i=;i<=m;i++)tax[i]+=tax[i-];
     for(int i=n;i>=;i--)sa[tax[a[b[i]]]--]=b[i];
 }

 bool comp(int r[],int a,int b,int k){
     return r[a]==r[b]&&r[a+k]==r[b+k];
 }

 void get_sa(int a[],int b[]){
     for(int i=;i<=n1;i++)s[i]=s1[i];
     s[n1+]='$';
     for(int i=;i<=n2;i++)s[n1+i+]=s2[i];
     n=n1+n2+;
     for(int i=;i<=n;i++)a[i]=s[i],b[i]=i;
     m=,sort(a,b);
     for(int j=,p=;p<n;j<<=,m=p){
         p=;
         for(int i=;i<=j;i++)b[++p]=n-j+i;
         for(int i=;i<=n;i++)if(sa[i]>j)b[++p]=sa[i]-j;
         sort(a,b);
         int *t=a;a=b;b=t;
         a[sa[]]=p=;
         for(int i=;i<=n;i++)
         a[sa[i]]=comp(b,sa[i],sa[i-],j)?p:++p;
     }
     for(int i=;i<=n;i++)rak[sa[i]]=i;
 }

 void get_height(){
     int i=,j=,k=;
     for(i=;i<=n;height[rak[i++]]=k){
         j=sa[rak[i]-];
         if(k)k--;
         while(s[i+k]==s[j+k])k++;
     }
 }

 int main(){
     scanf("%s%s%s",s1+,s2+,s3+);
     n1=strlen(s1+),n2=strlen(s2+),n3=strlen(s3+);
     get_sa(rak,tp),get_height();
     get_next(),solve();
     for(int i=;i<=n;i++){
         int *k=pos,s;
         do{
             k=lower_bound(k+,pos+pos[]+,i);
             s=*k;
             if(k-pos>pos[]){s=2e9;break;}
         }while(*k-n3+<i);
         f[i]=s-i;
     }
     for(int i=;i<=n;i++)
         if((sa[i]<=n1&&sa[i-]>n1)||(sa[i]>n1&&sa[i-]<=n1))
             ans=max(ans,min(height[i],min(f[sa[i]],f[sa[i-]])));
     printf("%d",ans);
 }

长啊！！（别人都是10几20毫秒，我竟然500ms+ MMP (▼ヘ▼#)）

T2：抵制克苏恩

期望dp，$f[k][a][b][c]$ 表示攻击$k$次后，还剩$a$个 1 血，$b$个2血，$c$个3血，英雄受伤害的期望，直接递推！

CODE：

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 using namespace std;

 int T,n,m,K,A,B,C;
 double f[][][][];

 double dp(int k,int a,int b,int c){
     if(k==)return ;
     c=min(c,-a-b);
     if(f[k][a][b][c])return f[k][a][b][c];
     f[k][a][b][c]+=1.0/(a+b+c+)*(dp(k-,a,b,c)+);
     if(a)f[k][a][b][c]+=1.0*a/(a+b+c+)*dp(k-,a-,b,c);
     if(b)f[k][a][b][c]+=1.0*b/(a+b+c+)*dp(k-,a+,b-,c+);
     if(c)f[k][a][b][c]+=1.0*c/(a+b+c+)*dp(k-,a,b+,c);
     return f[k][a][b][c];
 }

 int main(){
     scanf("%d",&T);
     while(T--){
         scanf("%d%d%d%d",&K,&A,&B,&C);
         printf("%.2f\n",dp(K,A,B,C));
     }
 }

T3：美味

主席树，类似 Trie 树的思想，按位贪心；

用一个变量$a$从高位走，0就走1,1就走0；

主席树维护值域和区间，判断，具体看代码：

 #include<iostream>
 #include<cstdio>
 using namespace std;

 int n,m,b,x,l,r,y,tot=,root[];
 int sum[],lch[],rch[];

 int read(){
     int x=;char ch=;
     while(ch<''||ch>''){ch=getchar();}
     while(ch>=''&&ch<=''){x=(x<<)+(x<<)+(ch^);ch=getchar();}
     return x;
 }

 void update(int &o,int last,int l,int r,int x){
     o=++tot;
     lch[o]=lch[last],rch[o]=rch[last],sum[o]=sum[last]+;
     if(r-l==)return;
     else{
         int mid=l+r>>;
         if(x<mid)update(lch[o],lch[last],l,mid,x);
         else update(rch[o],rch[last],mid,r,x);
     }
 }

 int query(int o1,int o2,int l,int r,int x,int y){
     if(l>=y||r<=x)return ;
     if(l>=x&&r<=y)return sum[o2]>sum[o1];
     int mid=l+r>>;
     int ans=;
     if(x<mid)ans|=query(lch[o1],lch[o2],l,mid,x,y);
     if(y>mid)ans|=query(rch[o1],rch[o2],mid,r,x,y);
     return ans;
 }

 int main(){
     n=read(),m=read();
     for(int i=;i<=n;i+=){
         x=read();
         update(root[i],root[i-],,,x);
         x=read();
         update(root[i+],root[i],,,x);
     }
     register int i,j,a;
     for(i=,a=;i<=m;++i,a=){
         b=read(),x=read(),l=read(),r=read();
         for(register int j=;~j;--j){
             if(b&(<<j)){
                 int L=max(a-x,),R=a+(<<j)--x;
                 if(!query(root[l-],root[r],,,L,R+))a^=(<<j);//没法走0就走1
             }else{
                 int L=max(a+(<<j)-x,),R=a+(<<j+)--x;
                 if(query(root[l-],root[r],,,L,R+))a^=(<<j);//能走1就走1
             }
         }
         printf("%d\n",a^b);
     }
 }

成功 1960ms 卡过！！呵呵……  ┐(´∇｀)┌