【NOIP2016】补题

今天突然想到自己居然还没把NOIP2016补完

简直是傻逼、、、

所以开始写

D1T1:模拟

D1T2：NOIP最难的一道题，首先求LCA

离线下，把观察员单独提出来

然后可以维护一个类似桶排序的东西

#include<bits/stdc++.h>
#define N 400010
using namespace std;
int n,m,top[N],wson[N],size[N],tot=,head[N],w[N],fa[N],d[N];
struct Edge{int u,v,next;}G[N<<],E[N<<];
int f[N<<],g[N<<],ans[N],ind[N][],out[N][];
inline void addedge(int u,int v){
    G[++tot].u=u;G[tot].v=v;G[tot].next=head[u];head[u]=tot;
    G[++tot].u=v;G[tot].v=u;G[tot].next=head[v];head[v]=tot;
}
int head1[N],tot2=;
inline void addedge2(int &u,int v){
    E[++tot2].u=u;E[tot2].v=v;E[tot2].next=u;u=tot2;
}
inline void addquery(int t,int x,int y,int lca){
    addedge2(ind[x][t],lca);
    if(y>)addedge2(out[fa[y]][t],lca);
}
inline int read(){
    int f=,x=;char ch;
    do{ch=getchar();if(ch=='-')f=-;}while(ch<''||ch>'');
    do{x=x*+ch-'';ch=getchar();}while(ch>=''&&ch<='');
    return f*x;
}
void dfs1(int u,int f){
    size[u]=;
    for(int i=head[u];i;i=G[i].next){
        int v=G[i].v;if(v==f)continue;
        fa[v]=u;d[v]=d[u]+;
        dfs1(v,u);
        size[u]+=size[v];
        if(size[v]>size[wson[u]])wson[u]=v;
    }
}
void dfs2(int u,int tp){
    top[u]=tp;if(wson[u])dfs2(wson[u],tp);
    for(int i=head[u];i;i=G[i].next){
        int v=G[i].v;if(v==fa[u]||v==wson[u])continue;
        dfs2(v,v);
    }
}
inline int qlca(int x,int y){
    while (top[x]!=top[y]){
        if(d[top[x]]<d[top[y]])swap(x,y);
        x=fa[top[x]];
    }
    if(d[x]>d[y])swap(x,y);return x;
}
void dfs3(int u){
    ans[u]-=f[w[u]+d[u]]+g[w[u]-d[u]+N];
    for(int i=ind[u][];i;i=E[i].next)f[E[i].v]++;
    for(int i=out[u][];i;i=E[i].next)f[E[i].v]--;
    for(int i=ind[u][];i;i=E[i].next)g[E[i].v+N]++;
    for(int i=out[u][];i;i=E[i].next)g[E[i].v+N]--;
    for(int i=head[u];i;i=G[i].next)if(G[i].v!=fa[u])dfs3(G[i].v);
    ans[u]+=f[w[u]+d[u]]+g[w[u]-d[u]+N];
}
int main(){
    n=read();m=read();
    for(int i=;i<n;i++){addedge(read(),read());}
    for(int i=;i<=n;i++)w[i]=read();
    dfs1(,);dfs2(,);
    for(int i=;i<=m;i++){
        int x=read(),y=read(),lca=qlca(x,y);
        if(d[x]-d[lca]==w[lca])ans[lca]--;
        addquery(,x,lca,d[x]);addquery(,y,lca,d[x]-*d[lca]);
    }
    dfs3();
    for(int i=;i<=n;i++)printf("%d ",ans[i]);puts("");
}

D1T3：其实是很简单的概率dp，但是转移的时候讨论比较恶心。

#include<bits/stdc++.h>
#define inf 1000000007
using namespace std;
int n,m,v,e,c[],d[];double dis[][];
double p[],dp[][][];
inline void floyed(){
    for(int k=;k<=v;k++)
    for(int i=;i<=v;i++)
    for(int j=;j<=v;j++)dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
}
inline int read(){
    int f=,x=;char ch;
    do{ch=getchar();if(ch=='-')f=-;}while(ch<''||ch>'');
    do{x=x*+ch-'';ch=getchar();}while(ch>=''&&ch<='');
    return f*x;
}
int main(){
    n=read();m=read();v=read();e=read();
    for(int i=;i<=n;i++)c[i]=read();for(int i=;i<=n;i++)d[i]=read();
    for(int i=;i<=n;i++)scanf("%lf",&p[i]);
    for(int i=;i<=v;i++)for(int j=;j<=v;j++)if(i==j)dis[i][j]=;else dis[i][j]=inf;
    for(int i=;i<=e;i++){
        int u,v;double w;scanf("%d%d%lf",&u,&v,&w);
        dis[u][v]=dis[v][u]=min(dis[u][v],w);
    }
    floyed();
for(int i=;i<=n;i++)for(int j=;j<=m;j++)dp[i][j][]=dp[i][j][]=inf;
    dp[][][]=;dp[][][]=;
    for(int i=;i<=n;i++)
        for(int j=;j<=m;j++){
            dp[i][j][]=min(dp[i-][j][]+dis[c[i]][c[i-]],dp[i-][j][]+dis[c[i]][c[i-]]*(-p[i-])+dis[c[i]][d[i-]]*p[i-]);
            if(j>){
                double tmp1=,tmp2=;
                tmp1=(dp[i-][j-][]+dis[d[i]][c[i-]])*p[i]+(dp[i-][j-][]+dis[c[i]][c[i-]])*(-p[i]);
                tmp2=(dp[i-][j-][]+dis[d[i]][d[i-]]*p[i-]+dis[d[i]][c[i-]]*(-p[i-]))*p[i]+(dp[i-][j-][]+dis[c[i]][d[i-]]*p[i-]+dis[c[i]][c[i-]]*(-p[i-]))*(-p[i]);
                dp[i][j][]=min(tmp1,tmp2);
            }
    }
    double ans=inf;for(int i=;i<=m;i++)ans=min(ans,min(dp[n][i][],dp[n][i][]));
    printf("%.2lf\n",ans);
}

D2T1：参见之前blog

D2T2：比较妙的一道题，证明完之后三个单调队列即可。

#include<bits/stdc++.h>
#define N 7100007
#define inf 0x7fffffff
using namespace std;
typedef long long ll;
int n,m,q,u,v,tot;
ll k,t,ans,x,y,delt;
ll Q[][N],l[],r[];
bool cmp(int x,int y){return x>y;}
inline int read(){
    int f=,x=;char ch;
    do{ch=getchar();if(ch=='-')f=-;}while(ch<''||ch>'');
    do{x=x*+ch-'';ch=getchar();}while(ch>=''&&ch<='');
    return f*x;
}
int main(){
    n=read();m=read();q=read();u=read();v=read();tot=read();
    for(int i=;i<=m+n;i++)Q[][i]=Q[][i]=Q[][i]=-inf;
    for(int i=;i<=n;i++)scanf("%lld",&Q[][i]);
    sort(Q[]+,Q[]+n+,cmp);
    l[]=,r[]=n;l[]=l[]=;r[]=r[]=;
    for(int i=;i<=m;i++){
        delt=-inf;for(int j=;j<=;j++)if(Q[j][l[j]]>delt)delt=Q[j][l[j]],k=j;
        ll t=delt+(i-)*q;l[k]++;
        if(i%tot==)printf("%lld ",t);
        x=t*u/v;y=t-x;
        Q[][++r[]]=x-i*q;
        Q[][++r[]]=y-i*q;
    }
    puts("");
    for(int i=;i<=m+n;i++){
        delt=-inf;for(int j=;j<=;j++)if(Q[j][l[j]]>delt)delt=Q[j][l[j]],k=j;
        t=delt+m*q;l[k]++;
        if(i%tot==)printf("%lld ",t);
    }
}

D2T3：状压dp，参考之前blog。