noip2016代码

----------------------------------------------------------------------------------

以下均为AC代码

-----------------------------------------------------------------------------------

Day1 T1

#include<cstdio>

#include<cstring>

#include<iostream>

using namespace std;

inline int read(){

    register int x=;bool f=;

    register char ch=getchar();

    while(ch<''||ch>''){if(ch=='-')f=;ch=getchar();}

    while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}

    return f?x:-x;

}

const int N=1e5+;

struct node{

    int f;

    char s[];

    node(){

        f=;

        memset(s,,sizeof s);

    }

}e[N<<];

int n,m;

int main(){

    n=read();m=read();

    for(int i=;i<=n;i++){

        e[i].f=read();

        scanf("%s",&e[i].s);

    }

    int zj=;

    for(int i=,opt,p;i<=m;i++){

        opt=read();p=read();p%=n;

        if(!e[zj].f){

            if(opt){

                if(zj+p<=n) zj+=p;

                else zj=zj+p-n;

            }

            else{

                if(zj<=p) zj=n-p+zj;

                else zj-=p;

            }

        }

        else{

            if(!opt){

                if(zj+p<=n) zj+=p;

                else zj=zj+p-n;

            }

            else{

                if(zj<=p) zj=n-p+zj;

                else zj-=p;

            }

        }

    }

    printf("%s",e[zj].s);

    return ;

}

Day1 T2

/*

题解：by shenben

如果观察点x可以观察到u->v路径上的点：

1、u->fa：t[x]+de[x]=de[u]

2、fa->v: t[x]-de[x]=dis[u,v]-de[v]

对于每条路径：用差分去维护每个点的信息

对于答案统计：可以做dfs序，然后用主席树解决“区间有多少个某数”的问题

*/

#include<cstdio>

#include<cstring>

#include<iostream>

#define me(a) memset(a,0,sizeof a)

using namespace std;

inline int read(){

    register int x=;bool f=;

    register char ch=getchar();

    while(ch<''||ch>''){if(ch=='-')f=;ch=getchar();}

    while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}

    return f?x:-x;

}

const int N=3e5+;

const int M=N*;

int le[N],ri[N],dep[N],f[N][],ans[N],t[N],tt[N],st[N],en[N],lca[N];

int sum[M],ls[M],rs[M],root[N*];

struct node{

    int v,next;

}e[N<<];

int n,m,now,num,tot,cnt,head[N];

void add(int x,int y){

    e[++tot].v=y;e[tot].next=head[x];head[x]=tot;

}

void dfs(int x,int fa){

    le[x]=++cnt;

    for(int i=head[x],v;i;i=e[i].next){

        if((v=e[i].v)!=fa){

            dep[v]=dep[x]+;

            f[v][]=x;

            dfs(v,x);

        }

    }

    ri[x]=cnt;

}

int calc_lca(int a,int b){

    if(dep[a]<dep[b]) swap(a,b);

    int t=dep[a]-dep[b];

    for(int i=;i<=;i++){

        if(<<i&t){

            a=f[a][i];

        }

    }

    if(a==b) return a;

    for(int i=;i>=;i--){

        if(f[a][i]!=f[b][i]){

            a=f[a][i];

            b=f[b][i];

        }

    }

    return f[a][];

}

void change(int &rt,int l,int r,int pos,int val){

    if(!pos) return ;

    if(!rt) rt=++num;

    sum[rt]+=val;

    if(l==r) return ;

    int mid=l+r>>;

    if(pos<=mid) change(ls[rt],l,mid,pos,val);

    else change(rs[rt],mid+,r,pos,val);

}

int query(int rt,int l,int r,int x,int y){

    if(!rt) return ;

    if(l==x&&r==y) return sum[rt];

    int mid=l+r>>;

    if(y<=mid) return query(ls[rt],l,mid,x,y);

    else if(x>mid) return query(rs[rt],mid+,r,x,y);

    else return query(ls[rt],l,mid,x,mid)+query(rs[rt],mid+,r,mid+,y);

}

void Cl(){

    num=;me(ls);me(rs);me(sum);me(root);

}

int main(){

    n=read();m=read();

    for(int i=,x,y;i<n;i++){

        x=read();y=read();

        add(x,y);add(y,x);

    }

    dep[]=;

    dfs(,);

    for(int j=;j<=;j++){

        for(int i=;i<=n;i++){

            f[i][j]=f[f[i][j-]][j-];

        }

    }

    for(int i=;i<=n;i++) t[i]=read();

    for(int i=;i<=m;i++) st[i]=read(),en[i]=read(),lca[i]=calc_lca(st[i],en[i]);

    for(int i=;i<=n;i++) tt[i]=t[i]+dep[i];

    for(int i=;i<=m;i++){

        now=dep[st[i]];

        change(root[now],,n,le[st[i]],);

        change(root[now],,n,le[f[lca[i]][]],-);

    }

    for(int i=;i<=n;i++) ans[i]+=query(root[tt[i]],,n,le[i],ri[i]);

    Cl();

    for(int i=;i<=n;i++) tt[i]=t[i]-dep[i]+n+;

    for(int i=;i<=m;i++){

        now=dep[st[i]]-(dep[lca[i]]<<)+n+;

        change(root[now],,n,le[en[i]],);

        change(root[now],,n,le[lca[i]],-);

    }

    for(int i=;i<=n;i++) ans[i]+=query(root[tt[i]],,n,le[i],ri[i]);

    for(int i=;i<=n;i++) printf("%d ",ans[i]);

    return ;

}

Day1 T3

#include<cstdio>

#include<cstring>

#include<iostream>

#define DB double

#define inf 0x3f3f3f3f

using namespace std;

inline int read(){

    register int x=;bool f=;

    register char ch=getchar();

    while(ch<''||ch>''){if(ch=='-')f=;ch=getchar();}

    while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}

    return f?x:-x;

}

const int N=;

const int M=;

int n,m,V,E,f[N][N];

int c[M],d[M];

DB K[M],dp[][M][];

void init(){

    memset(f,inf,sizeof f);

    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",&K[i]);

    for(int i=,x,y,z;i<=E;i++){

        x=read();y=read();z=read();

        f[x][y]=min(f[x][y],z);

        f[y][x]=min(f[y][x],z);

    }

}

void floyed(){

    for(int i=;i<=V;i++) f[i][i]=;

    for(int k=;k<=V;k++){

        for(int i=;i<=V;i++){

            for(int j=;j<=V;j++){

                if(i!=j&&i!=k&&j!=k){

                    f[i][j]=min(f[i][j],f[i][k]+f[k][j]);

                }

            }

        }

    }

}

void work(){

    //memset(dp,0x7fffffff,sizeof dp);对double类型慎用

    for(int i=;i<=;i++){

        for(int j=;j<=m;j++){

            dp[i][j][]=dp[i][j][]=inf;

        }

    }

    int now=;

    dp[][][]=;dp[][][]=;

    for(int i=;i<=n;i++){

        now^=;

        dp[now][][]=dp[now^][][]+f[c[i-]][c[i]];

        for(int j=;j<=min(i,m);j++){

            dp[now][j][]=min(dp[now^][j][]+f[c[i-]][c[i]],dp[now^][j][]+K[i-]*f[d[i-]][c[i]]+(-K[i-])*f[c[i-]][c[i]]);

            dp[now][j][]=min(dp[now^][j-][]+K[i]*f[c[i-]][d[i]]+(-K[i])*f[c[i-]][c[i]],dp[now^][j-][]+K[i-]*K[i]*f[d[i-]][d[i]]+(-K[i-])*K[i]*f[c[i-]][d[i]]+K[i-]*(-K[i])*f[d[i-]][c[i]]+(-K[i-])*(-K[i])*f[c[i-]][c[i]]);

        }

    }

    DB ans=inf;

    for(int i=;i<=m;i++) ans=min(ans,min(dp[now][i][],dp[now][i][]));

    printf("%.2lf",ans);

}

int main(){

    init();

    floyed();

    work();

    return ;

}

Day2 T1

#include<cstdio>

#include<ios>

using namespace std;

inline int read(){

    register int x=;bool f=;

    register char ch=getchar();

    while(ch<''||ch>''){if(ch=='-')f=;ch=getchar();}

    while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}

    return f?x:-x;

}

const int N=;

int n,m,k,T;

int C[N][N],S[N][N];

void prepare(){

    for(int i=;i<=;i++){

        C[i][]=%k;

        for(int j=;j<=i;j++){

            C[i][j]=(C[i-][j]+C[i-][j-])%k;

        }

    }

    S[][]=(C[][]==);

    for(int i=;i<=;i++){

        for(int j=;j<=i;j++){

            if(!j){

                S[i][j]=S[i-][j]+(C[i][j]==);

            }

            else if(i==j){

                S[i][j]=S[i][j-]+(C[i][j]==);

            }

            else{

                S[i][j]=S[i-][j]+S[i][j-]-S[i-][j-]+(C[i][j]==);

            }

        }

    }

}

int main(){

    T=read(),k=read();

    prepare();

    for(;T--;){

        n=read();m=read();

        printf("%d\n",S[n][min(n,m)]);

    }

    return ;

}

Day2 T2

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#define DB double

using namespace std;

inline int read(){

    register int x=;bool f=;

    register char ch=getchar();

    while(ch<''||ch>''){if(ch=='-')f=;ch=getchar();}

    while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}

    return f?x:-x;

}

const int N=1e5+,M=7e6+;

int n,m,q,u,v,t;

int mark,q0[N],q1[M],q2[M];

int main(){

    n=read();m=read();q=read();u=read();v=read();t=read();

    memset(q0,-0x3f,sizeof q0);

    memset(q1,-0x3f,sizeof q1);

    memset(q2,-0x3f,sizeof q2);

    int h0=,t0=n,h1=,t1=,h2=,t2=;

    for(int i=;i<n;i++) q0[i]=read();

    sort(q0,q0+n,greater<int>());

    DB p=(DB)u/(DB)v;

    for(int i=,x,x1;i<=m;i++){

        x=max(q0[h0],max(q1[h1],q2[h2]));

        if(x==q0[h0]) h0++;

        else if(x==q1[h1]) h1++;

        else h2++;

        x1=(x+=mark)*p;mark+=q;

        q1[t1++]=x1-mark;

        q2[t2++]=x-x1-mark;

        if(i%t) continue;

        printf("%d ",x);

    }

    putchar('\n');

    for(int i=,x;i<=n+m;i++){

        x=max(q0[h0],max(q1[h1],q2[h2]));

        if(x==q0[h0]) h0++;

        else if(x==q1[h1]) h1++;

        else h2++;

        if(i%t) continue;

        printf("%d ",x+mark);

    }

    return ;

}

Day2 T3

#include<cstdio>

#include<cstring>

#include<cmath>

#include<iostream>

#define DB double

using namespace std;

const int N=;

const DB eps=1e-;

int T,n,s[N][N],dp[<<N-];

DB fm,a,b,x[N],y[N];

int main(){

    ios::sync_with_stdio(false);

    for(cin>>T;T--;){

        memset(dp,0x3f,sizeof dp);

        int i,j,k,sta;

        cin>>n>>i;

        for(i=;i<=n;i++) cin>>x[i]>>y[i];

        for(i=;i<=n;i++){

            for(j=i+;j<=n;j++){

                fm=x[i]*x[j]*(x[i]-x[j]);

                a=x[j]*y[i]-x[i]*y[j];

                b=x[i]*x[i]*y[j]-x[j]*x[j]*y[i];

                s[i][j]=;

                if(a*fm<){

                    for(k=;k<=n;k++){

                        if(fabs(a*x[k]*x[k]+b*x[k]-y[k]*fm)<eps){

                            s[i][j]|=(<<k-);

                        }

                    }

                }

            }

        }

        dp[]=;

        for(sta=;sta<(<<n);sta++){

            if(sta<0x3f3f3f3f){

                for(i=;<<i&sta;i++);i++;

                dp[<<i-|sta]=min(dp[<<i-|sta],dp[sta]+);

                for(j=i+;j<=n;j++){

                    dp[sta|s[i][j]]=min(dp[sta|s[i][j]],dp[sta]+);

                }

            }

        }

        printf("%d\n",dp[sta-]);

    }

    return ;

}