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=0;
        for(int i=0;i<len;i++)
            if(s[i]=='4' || s[i]=='7') continue;
            else if(s[i]!='3') st[++top]=s[i];
            else 
            {
                if(st[top]=='1') top--; 
                else st[++top]='3';
            }
        if(top) for(int i=1;i<=top;i++) putchar(st[i]);
        else putchar('0');
        printf("
");
    }
}

枚举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][21];
int p,n;
int logg2[N];
void read(int &x)
{
    x=0; char c=getchar();
    while(!isdigit(c)) c=getchar();
    while(isdigit(c)) { x=x*10+c-'0'; c=getchar(); }
}
void stpre()
{
    p=log2(n);
    for(int j=1,k=2;j<=p;j++,k<<=1)
        for(int i=1;i+k-1<=n;i++)
            st[i][j]=min(st[i][j-1],st[i+(k>>1)][j-1]);
    for(int i=1;i<=n;i++) logg2[i]=log2(i);
}
int getmin(int s,int t)
{
    p=logg2[t-s+1];
    int len=1<<p;
    return min(st[s][p],st[t-len+1][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=1;i<=n;i++) read(st[i][0]);
    stpre();
    int Sa,Sb,Ta,Tb,Ma,Mb,Mi;
    int ans=0;
    int l,r,mid,tmp;
    for(int i=A+1;i<=n-A-k+1;i++)
    {
        Ta=i-1; Sb=i+k; 
        Tb=max(Sb+A-1,B+i-A+Sb-Ta-2);
        Ma=getmin(i-A,Ta); Mb=getmin(Sb,Tb); 
        l=1; r=i-A; tmp=min(Ma,Mb);
        while(l<=r)
        {
            mid=l+r>>1;
            Tb=max(Sb+A-1,B+mid+Sb-Ta-2);
            Ma=getmin(mid,Ta); Mb=getmin(Sb,Tb);
            tmp=max(tmp,min(Ma,Mb));
            if(Ma<Mb) l=mid+1;                
            else if(Ma==Mb) break;
            else r=mid-1;
        }
        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[19][19][2];
int len,f[19];
char X[19];
int a[19],cnt;
void kmp()
{
    len=strlen(X);
    int j;
    for(int i=1;i<len;i++)
    {
        j=f[i];
        while(j && X[i]!=X[j]) j=f[j];
        f[i+1]=X[i]==X[j] ? j+1 : 0;
    }
}
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=0;
    int k=lim ? a[dep] : 9;
    int j;
    for(int i=0;i<=k;i++)
    {
        j=w;
        while(j && X[j]-'0'!=i) j=f[j];
        if(X[j]-'0'==i) j++;
        res+=dfs(dep-1,j,lim&&(i==k),get||(j==len));
    }
    return lim ? res : dp[dep][w][get]=res;
}
LL query(LL lim)
{
    int cnt=0;
    while(lim)
    {
        a[++cnt]=lim%10;
        lim/=10;
    }
    memset(dp,-1,sizeof(dp));
    return dfs(cnt,0,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-1);
    if(query(R)<K) { puts("Hey,wake up!"); return 0; }
    LL l=L,r=R,ans,mid;
    while(l<=r)
    {
        mid=l+r>>1;
        if(query(mid)-tmp<K) l=mid+1;
        else ans=mid,r=mid-1;
    }
    printf("%I64d",ans);
}