NOIP2016DAY1题解

https://www.luogu.org/problem/lists?name=&orderitem=pid&tag=33%2C83


 

T1:玩具谜题

题解:

沙茶模拟

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int N=1e5+5,INF=1e9+5;
typedef long long ll;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
int n,m,a,s;
struct toy{
    int d;
    char s[15];
}t[N];
int main(){
    n=read();m=read();    
    for(int i=0;i<n;i++){t[i].d=read();scanf("%s",t[i].s);}
    int now=0;
    while(m--){ //printf("now %d\n",now);
        a=read();s=read();
        if(t[now].d==0){
            if(a==0) now=(now-s+n)%n;
            else now=(now+s)%n;
        }else{
            if(a==0) now=(now+s)%n;
            else now=(now-s+n)%n;
        }
    }
    printf("%s",t[now].s);
    return 0;
}
mengbier

T2:天天爱跑步

题解:

现在基本上能想出思路,但是写了三节课...............

对于一个点j,能被他观察到的也就是他的子树中 ds-dj=wj或者dt-dj=len(s,t)-wj

就是 ds=dj+wj或dt-len(s,t)=dj-wj

只有j是改变的,然后只要找每个点的子树中满足那两个条件的端点个数就可以了

统计个数直接开个权值数组,遇到s或t就加,遇到lca就减(之外的节点不可能被这个lca的s和t贡献了)

对于一条路线,在S和T分别打两个tag,并在lca处回收tag

因为一个节点会被多次,所以每个节点u开一个链表,保存有那些s是u,那些t是u,那些lca是u

但是这样lca(s,t)这个点就不会受到s或t的贡献了(加了后又回收了),所以把s或t的回收放在fa[lca]就可以了

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int N=3e5+5,M=3e5+5,INF=1e9+5,B=3e5+2;
typedef long long ll;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

int n,m,u,v,w[N],s,t;
struct edge{
    int v,ne;
}e[M<<1];
int cnt,h[N];
inline void ins(int u,int v){
    cnt++;
    e[cnt].v=v;e[cnt].ne=h[u];h[u]=cnt;
    cnt++;
    e[cnt].v=u;e[cnt].ne=h[v];h[v]=cnt;
}
struct thing{
    int w,ne;
}e1[N],e2[N],e3[N],e4[N];
int cnt1,cnt2,cnt3,cnt4,h1[N],h2[N],h3[N],h4[N];
inline void ins1(int u,int v,int w){
    cnt1++;
    e1[cnt1].w=w;e1[cnt1].ne=h1[u];h1[u]=cnt1;
    cnt3++;
    e3[cnt3].w=w;e3[cnt3].ne=h3[v];h3[v]=cnt3;
}
inline void ins2(int u,int v,int w){
    cnt2++;
    e2[cnt2].w=w;e2[cnt2].ne=h2[u];h2[u]=cnt2;
    cnt4++;
    e4[cnt4].w=w;e4[cnt4].ne=h4[v];h4[v]=cnt4;
}
int vis[N],deep[N],fa[N][21];
void dfs(int u){
    vis[u]=1;
    for(int j=1;j<=20;j++)
        fa[u][j]=fa[fa[u][j-1]][j-1];
    
    for(int i=h[u];i;i=e[i].ne){
        int v=e[i].v;
        if(vis[v]) continue;
        fa[v][0]=u;
        deep[v]=deep[u]+1;
        dfs(v);    
    }
}
int lca(int x,int y){
    if(deep[x]<deep[y]) swap(x,y);
    int bin=deep[x]-deep[y];
    for(int i=0;i<=20;i++) 
        if((1<<i)&bin) x=fa[x][i];
    
    for(int i=20;i>=0;i--)
        if(fa[x][i]!=fa[y][i]){
            x=fa[x][i];
            y=fa[y][i];
        }
    if(x==y) return x;
    else return fa[x][0];
}
int c1[N<<1],c2[N<<1];
int ans[N];
void dfsSol(int u,int fa){
    int p1=w[u]+deep[u],p2=deep[u]-w[u]+B;
    ans[u]-=c1[p1]+c2[p2];
    for(int i=h[u];i;i=e[i].ne) 
        if(e[i].v!=fa) dfsSol(e[i].v,u);
    for(int i=h1[u];i;i=e1[i].ne) c1[e1[i].w]++;
    for(int i=h3[u];i;i=e3[i].ne) c1[e3[i].w]--;
    for(int i=h2[u];i;i=e2[i].ne) c2[e2[i].w]++;
    for(int i=h4[u];i;i=e4[i].ne) c2[e4[i].w]--;
    ans[u]+=c1[p1]+c2[p2];
}
int main(){
    freopen("in.txt","r",stdin);
    n=read();m=read();
    for(int i=1;i<=n-1;i++) u=read(),v=read(),ins(u,v);
    for(int i=1;i<=n;i++) w[i]=read();
    dfs(1);
    for(int i=1;i<=m;i++){
        s=read();t=read();
        int a=lca(s,t),len=deep[s]+deep[t]-deep[a]-deep[a];
        ins1(s,fa[a][0],deep[s]);
        ins2(t,a,deep[t]-len+B);
        //printf("lca %d %d %d %d %d\n",i,s,t,a,len);
    }
    dfsSol(1,0);
    for(int i=1;i<n;i++) printf("%d ",ans[i]);printf("%d",ans[n]);
}
天天是谁

T3:换教室

题解:

真的是道好水的DP:f[i][j][0/1]表示到了时间段i申请了j次,转移看这次和上次申请成不成功行了

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int NN=2005,N=305;
const double INF=1e9+5;
typedef long long ll;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
int tn,tm,n,m,c[NN],d[NN],g[N][N],u,v,w;
double k[NN];
void floyd(){
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                g[i][j]=min(g[i][j],g[i][k]+g[k][j]);
}

double f[NN][NN][2];
void dp(){
    for(int i=1;i<=tn;i++)
        for(int j=0;j<=tm;j++) f[i][j][0]=f[i][j][1]=INF;
    f[0][0][0]=0;
    for(int i=1;i<=tn;i++)
        for(int j=0;j<=tm;j++){//printf("hi %d %d %.2f %.2f\n",i,j,f[i][j][0],f[i][j][1]);
            f[i][j][0]=min(f[i-1][j][0]+g[c[i-1]][c[i]],
                            f[i-1][j][1]+g[d[i-1]][c[i]]*k[i-1]+g[c[i-1]][c[i]]*(1-k[i-1]));
            if(j>0) f[i][j][1]=f[i-1][j-1][0]+g[c[i-1]][d[i]]*k[i]+g[c[i-1]][c[i]]*(1-k[i]);
            if(j>1) f[i][j][1]=min(f[i][j][1],f[i-1][j-1][1]+k[i-1]*(g[d[i-1]][d[i]]*k[i]+g[d[i-1]][c[i]]*(1-k[i]))+
                                (1-k[i-1])*(g[c[i-1]][d[i]]*k[i]+g[c[i-1]][c[i]]*(1-k[i])));
        }
    double ans=INF;
    for(int j=0;j<=tm;j++) ans=min(ans,min(f[tn][j][0],f[tn][j][1])); 
    printf("%.2f",ans);
}
int main(){
    //freopen("in.txt","r",stdin);
    tn=read();tm=read();n=read();m=read();
    for(int i=1;i<=tn;i++) c[i]=read();
    for(int i=1;i<=tn;i++) d[i]=read();
    for(int i=1;i<=tn;i++) scanf("%lf",&k[i]);
    for(int i=1;i<=n;i++) for(int j=i+1;j<=n;j++) g[i][j]=g[j][i]=INF;
    for(int i=1;i<=m;i++){
        u=read();v=read();w=read();
        if(u!=v) g[u][v]=g[v][u]=min(g[u][v],w);
    }
    floyd();
    dp();
}
啦啦啦期望

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值