Noip2016滚粗记QAQ

Posted zhouyis

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了Noip2016滚粗记QAQ相关的知识，希望对你有一定的参考价值。

day1

t1

XBG

#include<map>
#include<cstdio>
#include<string>
#include<string.h>
#include<iostream>
#include<algorithm>
using namespace std;
int d[123456];
char ty[123456][20];
int n,m;
//... 
int main(){
    freopen("toy.in","r",stdin);
    freopen("toy.out","w",stdout);
    scanf("%d %d",&n,&m);
    for(int i=0;i<n;i++)scanf("%d%s",d+i,ty[i]);
    //for(int i=1;i<=n;i++)d[i]=d[i]?1:-1;
    int now=0;
    while(m--){
        int op,s;
        scanf("%d %d",&op,&s);
        if(!d[now]){
            if(op)now=(now+s)%n;
            else now=((now-s)%n+n)%n;    
        }else{
            if(op)now=((now-s)%n+n)%n;
            else now=(now+s)%n;
        }
    }
    cout<<ty[now];
    return 0;
}

t2

暴力30

子任务Si=1

以1为根统计每个子树里面1的个数

子任务链的写错了。

45分

#include<map>
#include<cstdio>
#include<string>
#include<vector>
#include<string.h>
#include<iostream>
#include<algorithm>
using namespace std;
struct edge{
    int to,next;
}e[601234];
int cnt,last[301234],dep[301234],w[301234],ans[301234];
int n,m;
int s[301234],t[301234];
void insert(int a,int b){
    e[++cnt]=(edge){b,last[a]};last[a]=cnt;
    e[++cnt]=(edge){a,last[b]};last[b]=cnt;
}
int fa[301234][20];
void dfs(int v,int f){
    //fa[v][0]=f;
    //for(int i=1;i<=20;i++)fa[v][i]=fa[fa[v][i-1]][i-1];
    for(int i=last[v];i;i=e[i].next){
        int b=e[i].to;if(b!=f){
            dep[b]=dep[v]+1;
            dfs(b,v);
        }
    }
}
int an[301234];
void dfs2(int v,int t,int f){
    fa[v][0]=f;
    //for(int i=1;i<=20;i++)fa[v][i]=fa[fa[v][i-1]][i-1];
    for(int i=last[v];i;i=e[i].next){
        int b=e[i].to;if(b!=f){
            dfs2(b,t,v);
            if(b==t){
                an[v]=i;
            }
            if(an[b]){
                an[v]=i;
            }
        }
    }
}
void dfs3(int s,int t){
    int v=s,ti=0;
    ans[s]+=(ti==w[s]);
    while(s!=t){
        s=e[an[s]].to;
        ++ti;
        ans[s]+=(ti==w[s]);
    }
}
int siz[301234];
void dfs4(int v,int f){
    //fa[v][0]=f;
    //for(int i=1;i<=20;i++)fa[v][i]=fa[fa[v][i-1]][i-1];
    for(int i=last[v];i;i=e[i].next){
        int b=e[i].to;if(b!=f){
            dep[b]=dep[v]+1;
            dfs4(b,v);
            siz[v]+=siz[b];
            //fprintf(stderr,"Debug (%d)=%d\n",v,siz[v]);
        }
    }
}
vector<int>vec[301234];
int main(){
    freopen("running.in","r",stdin);
    freopen("running.out","w",stdout);
    scanf("%d %d",&n,&m);
    for(int i=1,u,v;i<n;i++){
        scanf("%d %d",&u,&v);
        insert(u,v);
    }
    //dfs(1,0);
    for(int i=1;i<=n;i++)scanf("%d",w+i);
    for(int i=1;i<=m;i++)scanf("%d %d",s+i,t+i);
    if(n<=1000){
        for(int i=1;i<=m;i++){
            memset(an,0,sizeof(an));
            dfs2(s[i],t[i],0);
            dfs3(s[i],t[i]);
        }
        for(int i=1;i<=n;i++)printf("%d ",ans[i]);
    }else
    if(n%10==4){
        for(int i=1;i<=n;i++)if(t[i]>=s[i])vec[i].push_back(s[i]);
        for(int i=1;i<=n;i++){
            if(i-w[i]>=0){
                int Ans=vec[i-w[i]].size(),cnt=0;
                for(int x=0;x<Ans;x++){
                    if(t[vec[i-w[i]][x]]<i)cnt++;
                }
                ans[i]+=(Ans-cnt);
            }
        }
        
        //for(int i=1;i<=n;i++)printf("%d ",ans[i]);
        for(int i=1;i<=n;i++)vec[i].clear();
        for(int i=1;i<=n;i++)if(t[i]<s[i])vec[i].push_back(s[i]);
        for(int i=1;i<=n;i++){
            if(i+w[i]<=n){
                int Ans=vec[i+w[i]].size(),cnt=0;
                //fprintf(stderr,"at %d\n",Ans);
                for(int x=0;x<Ans;x++){
                    if(t[vec[i+w[i]][x]]>i){
                        cnt++;
                    }
                }
                ans[i]+=(Ans-cnt);
            }
        }
        for(int i=1;i<=n;i++)printf("%d ",ans[i]);
    }else {
        for(int i=1;i<=n;i++)siz[t[i]]++;
        dep[1]=0; 
        dfs4(1,0);
        for(int i=1;i<=n;i++)printf("%d ",(w[i]==dep[i])*siz[i]);
    }
    return 0;//... //... //... 
}

t3

2^n*2^m暴力枚举所有申请和批准情况

 

 

暴力64，数组小了不然是72

要特判m<2

#include<map>
#include<cstdio>
#include<string>
#include<vector>
#include<string.h>
#include<iostream>
#include<algorithm>
using namespace std;
int n,m,v,e;
int c[333],d[333];
double k[333];/*
int min(int a,int b){
    return a<b?a:b;
}*/
struct zy{
    int i,cost;
}F[1234];
bool cmp(zy a,zy b){
    return a.cost>b.cost;
}
double brutef(){
    double ans=1e23;
    int _F=0;
    for(int i=1;i<=n;i++){
        int tmp=0;
        if(tmp>1)tmp=tmp+g[c[i-1]][c[i]]-g[c[i-1]][d[i]];
        if(tmp<n)tmp=tmp+g[c[i+1]][c[i]]-g[c[i+1]][d[i]];
        F[++_F]=(zy){i,-tmp};
    }
    sort(F+1,F+_F+1,cmp);
    int q[1234];
    memset(q,0,sizeof(q));
    for(int i=0;i<=m;i++){
        q[F[i].i]=1;double tmp2=0;
        for(int t=1;t<n;t++){
            int a=c[t],b=c[t+1];
            if(q[t])a=d[t];
            if(q[t+1])b=d[t+1];
            tmp2+=g[a][b];
        }
        //cout<<F[i].i<<","<<tmp2<<endl;
        ans=min(ans,tmp2);
    }
    return ans;
}
int main() {
    freopen("classroom.in","r",stdin);
    freopen("classroom.out","w",stdout);
    scanf("%d%d%d%d",&n,&m,&v,&e);
    for(int i=1;i<=n;i++)scanf("%d",c+i);
    for(int i=1;i<=n;i++)scanf("%d",d+i);
    bool xz2=1;
    for(int i=1;i<=n;i++)scanf("%lf",k+i),xz2=xz2&&(k[i]==1.0);
    
    for(int i=1;i<=v;i++)for(int j=1;j<=v;j++)g[i][j]=(i!=j)*(1<<20);
    for(int i=1;i<=e;i++){
        int a,b,w;
        scanf("%d %d %d",&a,&b,&w);
        g[a][b]=min(g[a][b],w);
        g[b][a]=min(g[b][a],w);
    }
    
    for(int p=1;p<=v;p++)for(int i=1;i<=v;i++)for(int j=1;j<=v;j++)g[i][j]=min(g[i][j],g[i][p]+g[p][j]);
    if(m<2){
        double ans=0,tmp=0;
        for(int i=1;i<n;i++){
            ans+=(g[c[i]][c[i+1]]);
        }
        tmp=ans;
        if(m)
        for(int x=1;x<=n;x++){
            //cout<<x<<" :\n";
            int res=0;
            for(int i=1;i<n;i++){
                int a=c[i],b=c[i+1];
                if(i==x)a=d[x];
                if(i+1==x)b=d[x];
                res+=g[a][b];
                
            }
            //cout<<res*k[x]+tmp*(1-k[x])<<endl;
            ans=min(ans,res*k[x]+tmp*(1-k[x]));
        }
        printf("%.2lf",(double)ans);
    }else if(n<=20&&m<=10){
        double ans=1e23;
        for(int i=0;i<(1<<n);i++){
            int p[30],q[33],cnt=0;
            for(int j=1;j<=n;j++)if(i&(1<<(j-1)))q[++cnt]=j;
            if(cnt<=m){
            //cerr<<"safe";
            memset(p,0,sizeof(p));
                double tmp=0;
                for(int j=0;j<(1<<cnt);j++){
                    double prob=1,tmp2=0;
                    for(int t=1;t<=cnt;t++){
                        if(j&(1<<(t-1))){
                            prob=prob*k[q[t]];
                            p[q[t]]=1;
                        }else {
                            prob=prob*(1-k[q[t]]);
                            p[q[t]]=0;
                        }
                    }
                    /*printf("(%d,%d)\n",i,j);
                    for(int t=1;t<=cnt;t++)cout<<q[t]<<"!!";cout<<endl;*/
                    for(int t=1;t<n;t++){
                        int a=c[t],b=c[t+1];
                        if(p[t])a=d[t];
                        if(p[t+1])b=d[t+1];
                        tmp2+=g[a][b];
                    }
                    //cout<<tmp2<<endl;;
                    tmp+=tmp2*prob;
                }
                //cout<<"---"<<tmp<<endl;
                ans=min(ans,tmp);
                }
            }//... 
        
        printf("%.2lf",(double)ans);
    }else{
        printf("%.2lf",(double)brutef());
    }
    
}

day2

t1

XBG

#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
int t,k,n,m,c[2333][2333],s[2333][2333];
int main(){
    freopen("problem.in","r",stdin);
    freopen("problem.out","w",stdout);
    scanf("%d%d",&t,&k);
    for(int i=0;i<=2000;i++)c[i][1]=i%k;
    for(int i=1;i<=2000;i++)
        for(int j=2;j<=i;j++)
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%k;
    for(int i=1;i<=2000;i++)
        for(int j=1;j<=2000;j++){
            s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+(j<=i)*(!c[i][j]);
        }
    while(t--){
        scanf("%d %d",&n,&m);
        /*for(int i=1;i<=n;i++,puts(""))
            for(int j=1;j<=i;j++)printf("%d ",s[i][j]);
        */printf("%d\n",s[n][m]);
    }
    return 0;
}

t2

拿三个队列，一个存x，一个存px，一个存x-px。

然后这三个队列是单调递减的，具体证明较简单。

#include<queue>
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
int n,m,q,u,v,t;double p;
int a[123456];
int PQ[123456],_PQ,od[7123456];
namespace solution{
    int x[7123456][2],y[7123456][2],z[7123456][2];
    int hx,hy,hz,tx,ty,tz;
    
    void init(){hx=hy=hz=1;tx=n;ty=tz=0;}
    void debug(int tim){
        for(int i=hx;i<=n;i++)printf("%d ",x[i][0]+q*(tim-x[i][1]));
        printf("!");for(int i=hy;i<=ty;i++)printf("%d ",y[i][0]+q*(tim-y[i][1]));
        printf("!");for(int i=hz;i<=tz;i++)printf("%d ",z[i][0]+q*(tim-z[i][1]));
        printf("!");
        puts("");
    }
    //wo hui zuo d2t2 ei
    void getres(){//debug(m);
        for(int i=1;i<=n+m;i++){
            int c=-12345;int *I;    
            if(hx<=n)if(x[hx][0]+q*(m-x[hx][1])>c)c=x[hx][0]+q*(m-x[hx][1]),I=&hx;
            if(hy<=ty)if(y[hy][0]+q*(m-y[hy][1])>c)c=y[hy][0]+q*(m-y[hy][1]),I=&hy;
            if(hz<=tz)if(z[hz][0]+q*(m-z[hz][1])>c)c=z[hz][0]+q*(m-z[hz][1]),I=&hz;
            ++*I;
            if(!(i%t))printf("%d ",c);
        }
    }
    int main(){
        sort(a+1,a+n+1,greater<int>());//nlogn
        for(int i=1;i<=n;i++)x[i][0]=a[i];
        p=u;p/=v;
        init();
        for(int i=1;i<=m;i++){
            int c=-12345;
            int *I;
            if(hx<=n)if(x[hx][0]+q*(i-1-x[hx][1])>c)c=x[hx][0]+q*(i-1-x[hx][1]),I=&hx;
            if(hy<=ty)if(y[hy][0]+q*(i-1-y[hy][1])>c)c=y[hy][0]+q*(i-1-y[hy][1]),I=&hy;
            if(hz<=tz)if(z[hz][0]+q*(i-1-z[hz][1])>c)c=z[hz][0]+q*(i-1-z[hz][1]),I=&hz;
            int A,B;
            A=p*c;
            B=c-A;
            //printf("cut %d into(%d,%d)\n",c,A,B);
            y[++ty][0]=A;
            z[++tz][0]=B;
            z[tz][1]=y[ty][1]=i;
            ++*I;
            od[i]=c;
            //cerr<<"h..="<<hx<<endl; 
            //debug(i);
        }
        for(int i=t;i<=m;i+=t)printf("%d ",od[i]);puts("");
        
        getres();
        return 0;
    }
}
int main(){
    freopen("earthworm.in","r",stdin);
    freopen("earthworm.out","w",stdout);
    //cerr<<sizeof(solution::x)*3;
    scanf("%d%d%d%d%d%d",&n,&m,&q,&u,&v,&t);
    for(int i=1;i<=n;i++)scanf("%d",a+i);
    if(n<=1000&&m<=1000){
        for(int i=1;i<=n;i++)PQ[i]=a[i];_PQ=n;
        for(int _t=1;_t<=m;_t++){
            int mx=1;
            for(int i=2;i<=_PQ;i++)if(PQ[i]>PQ[mx])mx=i;
            int x=PQ[mx];
            PQ[mx]=(int)u*1.0/v*x;
            PQ[++_PQ]=x-PQ[mx];
            //printf("cut %d into(%d,%d)\n",x,PQ[mx],PQ[_PQ]);
            for(int i=1;i<_PQ;i++)if(i!=mx)PQ[i]+=q;
            od[_t]=x;
        }
        for(int i=t;i<=m;i+=t)printf("%d ",od[i]);
        sort(PQ+1,PQ+n+m+1,greater<int>());
        puts("");
        for(int i=t;i<=n+m;i+=t)printf("%d ",PQ[i]);
    }else 
        solution::main();
    return 0;
}

t3

枚举两个点作抛物线，预处理能过什么点

状压（即或的值为2^n-1)

注意处理单点，分数未知

riio写错了

 

总结：d1t2 d2t3暴力写错了很伤

d2t2很像我之前写的QingdaoRegional的G

最后100+45+64+100+100+?未完待续

以上是关于Noip2016滚粗记QAQ的主要内容，如果未能解决你的问题，请参考以下文章

[日常] NOIP 2017滚粗记

NOIP2016滚粗记

noip2016滚粗记

NOIP2017 滚粗记

NOI WC2016滚粗记

noip2017滚粗记