2017北京国庆刷题Day6 morning

期望得分：100+100+20=220

实际得分：100+100+20=220

模拟栈

#include<cstdio>
#include<cstring>
using namespace std;
#define N 10002
char s[N],st[N];
int top;
int main()
{
    freopen("kakutani.in","r",stdin);
    freopen("kakutani.out","w",stdout);
    int n,len,lt;
    scanf("%d",&n);
    while(n--)
    {
        scanf("%s",s);
        len=strlen(s);
        top=;
        for(int i=;i<len;i++)
            if(s[i]=='' || s[i]=='') continue;
            else if(s[i]!='') st[++top]=s[i];
            else
            {
                if(st[top]=='') top--;
                else st[++top]='';
            }
        if(top) for(int i=;i<=top;i++) putchar(st[i]);
        else putchar('');
        printf("\n");
    }
}

枚举i为间隔K个录像的左端点，那么间隔录像为[i,i+k-1]

设第一段为[Sa,Ta]，第二段为[Sb,Tb]，Ma为min[Sa,Ta],Mb为min[Sb,Tb]

随着Sa的左移，Ma单调不增

随着Tb的右移，Mb单调不增

如果枚举Sa，则有以下式子： Ta-Sa+1+Tb-Sb+1>=B

即Tb>=B+Sa+Sb-Ta-2

因为Tb右移，Mb单调不增，所以Tb取等号最优

所以二分Sa的位置，用st表查询Ma，Mb

#include<cmath>
#include<cstdio>
#include<iostream>
#include<algorithm>
#define N 1000001
using namespace std;
int st[N][];
int p,n;
int logg2[N];
void read(int &x)
{
    x=; char c=getchar();
    while(!isdigit(c)) c=getchar();
    while(isdigit(c)) { x=x*+c-''; c=getchar(); }
}
void stpre()
{
    p=log2(n);
    for(int j=,k=;j<=p;j++,k<<=)
        for(int i=;i+k-<=n;i++)
            st[i][j]=min(st[i][j-],st[i+(k>>)][j-]);
    for(int i=;i<=n;i++) logg2[i]=log2(i);
}
int getmin(int s,int t)
{
    p=logg2[t-s+];
    int len=<<p;
    return min(st[s][p],st[t-len+][p]);
}
int main()
{
    freopen("dota.in","r",stdin);
    freopen("dota.out","w",stdout);
    int A,B,k;
    read(n); read(A); read(B); read(k);
    for(int i=;i<=n;i++) read(st[i][]);
    stpre();
    int Sa,Sb,Ta,Tb,Ma,Mb,Mi;
    int ans=;
    int l,r,mid,tmp;
    for(int i=A+;i<=n-A-k+;i++)
    {
        Ta=i-; Sb=i+k;
        Tb=max(Sb+A-,B+i-A+Sb-Ta-);
        Ma=getmin(i-A,Ta); Mb=getmin(Sb,Tb);
        l=; r=i-A; tmp=min(Ma,Mb);
        while(l<=r)
        {
            mid=l+r>>;
            Tb=max(Sb+A-,B+mid+Sb-Ta-);
            Ma=getmin(mid,Ta); Mb=getmin(Sb,Tb);
            tmp=max(tmp,min(Ma,Mb));
            if(Ma<Mb) l=mid+;
            else if(Ma==Mb) break;
            else r=mid-;
        }
        ans=max(ans,tmp);
    }
    printf("%d",ans);
}

数位DP

dp[i][j][0/1] 前i位匹配到X的第j位，是否已经包含1个X的数的个数

二分，计算<=mid的数里的答案

其中的匹配用kmp

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long LL;
LL L,R,K;
LL dp[][][];
int len,f[];
char X[];
int a[],cnt;
void kmp()
{
    len=strlen(X);
    int j;
    for(int i=;i<len;i++)
    {
        j=f[i];
        while(j && X[i]!=X[j]) j=f[j];
        f[i+]=X[i]==X[j] ? j+ : ;
    }
}
LL dfs(int dep,int w,bool lim,bool get)
{
    if(!dep) return get;
    if(!lim && ~dp[dep][w][get]) return dp[dep][w][get];
    LL res=;
    int k=lim ? a[dep] : ;
    int j;
    for(int i=;i<=k;i++)
    {
        j=w;
        while(j && X[j]-''!=i) j=f[j];
        if(X[j]-''==i) j++;
        res+=dfs(dep-,j,lim&&(i==k),get||(j==len));
    }
    return lim ? res : dp[dep][w][get]=res;
}
LL query(LL lim)
{
    int cnt=;
    while(lim)
    {
        a[++cnt]=lim%;
        lim/=;
    }
    memset(dp,-,sizeof(dp));
    return dfs(cnt,,true,false);
}
int main()
{
    freopen("spenum.in","r",stdin);
    freopen("spenum.out","w",stdout);
    scanf("%I64d%I64d%s%I64d",&L,&R,X,&K);
    kmp();
    LL tmp=query(L-);
    if(query(R)<K) { puts("Hey,wake up!"); return ; }
    LL l=L,r=R,ans,mid;
    while(l<=r)
    {
        mid=l+r>>;
        if(query(mid)-tmp<K) l=mid+;
        else ans=mid,r=mid-;
    }
    printf("%I64d",ans);
}