2017/7/24 离线赛

T1 插队

    很显然，题目的要求是对于第$i$只猴子，从他开始从后好往前数找到第二个比他高的猴子，然后再+1。比赛时，我想，既然要求第二个，那我们就要先求第一个，于是我敲了一个二分+线段树的$O(nlog^2n)$，求一次还行，但是求两次就爆了。比赛完，我看了别人的代码，发现可以用单调栈求比自己大的第一个数，只要$O(n)$就行了。

vector<int>e[M];
int n,top;
int h[M],stk[M],ans[M];
int main(){
    Rd(n);
    for(int i=1;i<=n;i++){
        Rd(h[i]);
        while(top&&h[i]>h[stk[top-1]])top--;
        e[stk[top-1]].push_back(i);
        stk[top++]=i;
    }
    top=0;
    for(int i=1;i<=n;i++){
        for(int j=0;j<e[i].size();j++){
            int t=e[i][j];
            while(top&&h[t]>h[stk[top-1]])top--;
            ans[t]=stk[top-1];
        }
        while(top&&h[i]>h[stk[top-1]])top--;
        stk[top++]=i;
    }
    for(int i=1;i<=n;i++){
        Pt(ans[i]+1);
        putchar('\n');
    }
    return 0;
}

T2 比赛

    比赛的时候想到要用$dp[i]$表示点$i$到达根节点的最小花费，但一直被一个点卡住，所以就不往下想了，后来发现自己题目看错了QAQ。我以为一张地图可以指明$i$到根的任意$k$条路径，然后就炸了，而实际是指明从$i$向上$k$条路径，然后用倍增维护一下dp值即可。

typedef pair<int,int>P;
vector<P>mp[N];
struct edge{
    int v,nxt;
}e[N<<1];
int n,m,q;
int stk[N];
int head[N],cnt;
ll ans[N],BZ[N][Max];
void add_edge(int u,int v){
    e[++cnt].v=v;e[cnt].nxt=head[u];head[u]=cnt;
}
void dfs(int x,int t,int d){
    stk[d]=x;
    if(x>1){
        for(int i=0;i<mp[x].size();i++){
            int k=mp[x][i].first,v=mp[x][i].second,t=0;
            k=min(k,d-1);
            while((1<<(t+1))<=k)t++;
            ans[x]=min(ans[x],min(BZ[stk[d-1]][t],BZ[stk[d-k-1+(1<<t)]][t])+v);
        }
    }else ans[x]=0;
    BZ[x][0]=ans[x];
    for(int i=1;(1<<i)<=d;i++)
    BZ[x][i]=min(BZ[x][i-1],BZ[stk[d-(1<<(i-1))]][i-1]);
    for(int i=head[x];~i;i=e[i].nxt)
    if(e[i].v!=t)dfs(e[i].v,x,d+1);
}
int main(){
    int a,b,x,k,v;
    memset(head,-1,sizeof(head));
    Rd(n);Rd(m);Rd(q);
    for(int i=1;i<n;i++){
        Rd(a);Rd(b);
        add_edge(a,b);
        add_edge(b,a);
    }
    for(int i=1;i<=m;i++){
        Rd(x);Rd(k);Rd(v);
        mp[x].push_back(P(k,v));
    }
    memset(ans,63,sizeof(ans));
    dfs(1,-1,1);
    for(int i=1;i<=q;i++){
        Rd(x);
        Pt(ans[x]);
        putchar('\n');
    }
    return 0;
}

T3 采油区域

    这道题刚看到时没有啥思路，于是维护了个二维前缀和，求出所有正方形的值，不知道该怎么版，于是将值从大到小sort了一波，水到了75分_(:зゝ∠)_，我也是很服自己。看了题解才知道这要用$dp$，还要分类讨论，于是就有了以下代码

int n,m,k,ans;
int sum[M][M];
int res[M][M];
int a[4][M][M];//0为左上,1为右上,2为左下,3为右下 
int rx[]={-1,-1,1,1},ry[]={-1,1,-1,1};
void calc(int x){
    if(x==0){
        for(int i=1;i+k-1<=n;i++)
        for(int j=1;j+k-1<=m;j++)
        a[x][i][j]=max(res[i][j],max(a[x][i+rx[x]][j],a[x][i][j+ry[x]]));
    }else if(x==1){
        for(int i=1;i+k-1<=n;i++)
        for(int j=m-k+1;j>=1;j--)
        a[x][i][j]=max(res[i][j],max(a[x][i+rx[x]][j],a[x][i][j+ry[x]]));
    }else if(x==2){
        for(int i=n-k+1;i>=1;i--)
        for(int j=1;j+k-1<=m;j++)
        a[x][i][j]=max(res[i][j],max(a[x][i+rx[x]][j],a[x][i][j+ry[x]]));
    }else{
        for(int i=n-k+1;i>=1;i--)
        for(int j=m-k+1;j>=1;j--)
        a[x][i][j]=max(res[i][j],max(a[x][i+rx[x]][j],a[x][i][j+ry[x]]));
    }
}
int main(){
    int x;
    Rd(n);Rd(m);Rd(k);
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++){
        Rd(x);
        sum[i][j]=sum[i-1][j-1]+(sum[i-1][j]-sum[i-1][j-1])+(sum[i][j-1]-sum[i-1][j-1])+x;
    }
    for(int i=1;i+k-1<=n;i++)
    for(int j=1;j+k-1<=m;j++){
        int x=i+k-1,y=j+k-1;
        res[i][j]=sum[x][y]-sum[i-1][j-1]-(sum[x][j-1]-sum[i-1][j-1])-(sum[i-1][y]-sum[i-1][j-1]);
    }
    for(int i=0;i<4;i++)calc(i);
    for(int i=1;i+k-1<=n;i++)
    for(int j=1;j+k+k-1<=m;j++)
    ans=max(ans,a[0][i][j]+a[1][i][j+k]+a[3][i+k][1]);
    for(int i=k+1;i+k-1<=n;i++)
    for(int j=1;j+k+k-1<=m;j++)
    ans=max(ans,a[2][i][j]+a[3][i][j+k]+a[1][i-k][1]);
    for(int i=1;i+k+k-1<=n;i++)
    for(int j=1;j+k-1<=m;j++)
    ans=max(ans,a[0][i][j]+a[2][i+k][j]+a[3][1][j+k]);
    for(int i=1;i+k+k-1<=n;i++)
    for(int j=k+1;j+k-1<=m;j++)
    ans=max(ans,a[1][i][j]+a[3][i+k][j]+a[2][1][j-k]);
    for(int i=k+1;i+k+k-1<=n;i++)
    for(int j=1;j+k-1<=m;j++)
    ans=max(ans,res[i][j]+a[1][i-k][1]+a[3][i+k][1]);
    for(int i=1;i+k-1<=n;i++)
    for(int j=k+1;j+k+k-1<=m;j++)
    ans=max(ans,res[i][j]+a[2][1][j-k]+a[3][1][j+k]);
    Pt(ans);
    putchar('\n');
    return 0;
}

总之呢，这次比赛学到了很多，一个是单调栈，还有一个是倍增转移dp。