[NOIP] 记NOIP2017

最新推荐文章于 2021-01-03 02:25:05 发布
最新推荐文章于 2021-01-03 02:25:05 发布
阅读量538 收藏
点赞数 1
分类专栏： Presentation
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Coldef/article/details/78607956
版权

我好像没什么长进啊

D0

和去年一样在浪

D1

T1一眼不可做，两眼傻逼题？裴蜀定理推一推就是a(b-1)-b
T2大模拟
T3先判零环的情况，然后令 fi,j 表示到第 i 个点，当前距离为 disi+j 的方案数， disi 表示起点到 i 的最短路
重构一张图，如果 disi+wi,j=disj 那么连一条 i→j 的边，按照拓扑序转移就行了

表示度数数组没有清空，被卡了30分啊啊啊啊啊啊啊啊！！！

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <queue>

using namespace std;

typedef pair<int,int> par;

const int N=100010,inf=1<<30;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}

inline void rea(int &x){
  char c=nc(); x=0;
  for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}

int t,n,m,k,p,cnt,G[N];
struct iedge{
  int t,nx,w;
}E[N<<2],iE[N<<2];

int iG[N],icnt,rG[N];

void iclear(){
  memset(iG,0,sizeof(iG)); icnt=0;
  memset(G,0,sizeof(G)); cnt=0;
  memset(rG,0,sizeof(rG));
}

inline void addedge(int x,int y,int z){
  E[++cnt].t=y; E[cnt].nx=G[x]; E[cnt].w=z; G[x]=cnt;
  E[++cnt].t=x; E[cnt].nx=rG[y]; rG[y]=cnt; E[cnt].w=z;
}

inline void iaddedge(int x,int y){
  iE[++icnt].t=y; iE[icnt].nx=iG[x]; iG[x]=icnt;
}

int dis[N],rdis[N],du[N];
priority_queue<par> Q;

inline void dij(int s,int *cG,int *cdis){
  for(int i=1;i<=n;i++) cdis[i]=inf;
  cdis[s]=0; Q.push(par(0,s));
  while(!Q.empty()){
    int x=Q.top().second,y=-Q.top().first; Q.pop();
    if(y!=cdis[x]) continue;
    for(int i=cG[x];i;i=E[i].nx)
      if(cdis[E[i].t]>cdis[x]+E[i].w){
    cdis[E[i].t]=cdis[x]+E[i].w;
    Q.push(par(-cdis[E[i].t],E[i].t));
      }
  }
}

int vis[N],low[N],dfn[N],sta[N],tp,tms,imx;

void tarjan(int x){
  low[x]=dfn[x]=++tms;
  vis[x]=1; sta[++tp]=x;
  for(int i=iG[x];i;i=iE[i].nx){
    if(!vis[iE[i].t]) tarjan(iE[i].t);
    if(vis[iE[i].t]==1) low[x]=min(low[x],low[iE[i].t]);
  }
  if(low[x]==dfn[x]){
    int cura=inf,curb=inf,cur,ccnt=0;
    while(1){
      cur=sta[tp--]; vis[cur]=2; 
      ccnt++;
      cura=min(cura,dis[cur]); curb=min(curb,rdis[cur]);
      if(cur==x || !tp) break;
    }
    if(ccnt>1) 
      imx=min(imx,cura+curb);
  }
}

inline bool nooo(){
  dij(n,rG,rdis);
  for(int i=1;i<=n;i++) vis[i]=low[i]=dfn[i]=0;
  tms=tp=0; imx=inf;
  for(int i=1;i<=n;i++)
    if(!vis[i]) tarjan(i);
  return imx<=dis[n]+k;
}

int f[N][55];
int q[N],l,r;

int main(){
  rea(t);
  while(t--){
    iclear();
    rea(n); rea(m); rea(k); rea(p);
    for(int i=1,x,y,z;i<=m;i++){
      rea(x); rea(y); rea(z); addedge(x,y,z);
      if(z==0) iaddedge(x,y);
    }
    dij(1,G,dis); 
    if(nooo()){
      puts("-1"); continue;
    }
    for(int i=1;i<=n;i++)
      for(int j=0;j<=k;j++){
    if(dis[i]+j+rdis[i]>dis[n]+k) break;
    f[i][j]=0;
      }
    for(int i=1;i<=n;i++) du[i]=0;
    for(int x=1;x<=n;x++)
      for(int i=G[x];i;i=E[i].nx)
    if(dis[E[i].t]==dis[x]+E[i].w) du[E[i].t]++;
    l=1; r=0;
    q[r=1]=1;
    while(l<=r){
      int x=q[l++];
      for(int i=G[x];i;i=E[i].nx)
    if(dis[E[i].t]==dis[x]+E[i].w){
      du[E[i].t]--;
      if(!du[E[i].t]) q[++r]=E[i].t;
    }
    }
    f[1][0]=1;
    for(int j=0;j<=k;j++){
      for(int ii=1;ii<=r;ii++){
    int x=q[ii];
    int cur=f[x][j];
    if(!cur || dis[x]+j+rdis[x]>dis[n]+k) continue;
    for(int i=G[x];i;i=E[i].nx){
      int dd=dis[x]+j+E[i].w-dis[E[i].t];
      if(dd<=k) f[E[i].t][dd]=(f[E[i].t][dd]+f[x][j])%p;
    }
      }
    }
    int ans=0;
    for(int i=0;i<=k;i++) ans=(ans+f[n][i])%p;
    printf("%d\n",ans);
  }
  return 0;
}

D1就这样吧…

D2

T1并查集送分
T2好像写了个三进制的状压DP，复杂度是 O(4n)
按层DP， fi,S
令三进制状态 S 中，为0的位置表示这个点不在前 i 层中，为 1 表示在 1∼i−1 层中， 2 表示在第 i 层，这样转移就好了……

数组开小-10

#include <cstdio>
#include <iostream>
#include <algorithm>

using namespace std;

int n,m;
int lk[15][15],pw[15],val[15];
int f[15][531445];
int ttt[531445],lst[4150],ttt1[4150],sz[531445];

inline int iinn(int x,int y){
    return (x/pw[y-1])%3;
}

inline int cst(int S){
    int ret=0;
    for(;S;S-=(1<<(lst[S]-1))){
        if(val[lst[S]]==(1<<30)) return 1<<30;
        ret+=val[lst[S]];
    }
    return ret;
}

int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++) lk[i][j]=1<<30;
    for(int i=1,x,y,z;i<=m;i++){
        scanf("%d%d%d",&x,&y,&z);
        lk[x][y]=lk[y][x]=min(lk[x][y],z);
    }
    pw[0]=1; for(int i=1;i<=n;i++) pw[i]=3*pw[i-1];
    for(int i=1;i<=n;i++)
        for(int j=0;j<pw[n];j++) f[i][j]=1<<30;

    for(int i=0;i<pw[n];i++)
        for(int j=1;j<=n;j++)
            if(!iinn(i,j)) ttt[i]|=1<<(j-1);
            else sz[i]++;

    for(int i=1;i<(1<<n);i++)
        for(int j=1;j<=n;j++)
            if((i>>(j-1))&1) ttt1[i]+=pw[j-1]*2;

    lst[1]=1;
    for(int i=2;i<(1<<n);i++)
        lst[i]=lst[i>>1]+1;

    for(int i=1;i<=n;i++) f[1][pw[i-1]*2]=0;
    int ans=1<<30;
    for(int i=1;i<=n;i++)
        for(int S=0;S<pw[n];S++){
            if(f[i][S]>=ans) continue;
            if(sz[S]==n) ans=min(ans,f[i][S]);
            for(int j=1;j<=n;j++)
                if(!iinn(S,j)){
                    val[j]=1<<30;
                    for(int k=1;k<=n;k++)
                        if(iinn(S,k)==2) val[j]=min(val[j],lk[k][j]);
                }
            int nxt=0,SS=ttt[S];
            for(int j=1;j<=n;j++)
                if(iinn(S,j)) nxt+=pw[j-1];
            for(int ss=SS;ss;ss=(ss-1)&SS){
                int cur=cst(ss);
                if(cur==(1<<30)) continue;
                f[i+1][nxt+ttt1[ss]]=min(f[i+1][nxt+ttt1[ss]],f[i][S]+cur*i);
            }
        }
    printf("%d\n",ans);
    return 0;
}

T3我是个傻逼啊
暴力打挂，没开longlong，-30

大力平衡树模拟就好了…

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstdlib>
#define fi first
#define se second

using namespace std;

typedef long long ll;

const int N=2000010;

struct node{
  node *ls,*rs;
  int x,y,val,size,key,tag;
}pool[N],*t;

int n,m,q;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}

inline void read(int &x){
  char c=nc(); x=0;
  for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}

node *ro[N],*L;

inline node *crt(int x,int y,int k,int v){
  node *ret=t++;
  ret->x=x; ret->y=y; ret->val=v; ret->key=k;
  return ret;
}

inline void Up(node *x){
  x->size=(x->ls?x->ls->size:0)+(x->rs?x->rs->size:0)+1;
}

inline void Push(node *x){
  if(!x->tag) return ;
  if(x->ls)
    x->ls->tag+=x->tag,x->ls->key+=x->tag;
  if(x->rs)
    x->rs->tag+=x->tag,x->rs->key+=x->tag;
  x->tag=0;
}

inline node *Merge(node *x,node *y){
  if(!x || !y) return x?x:y;
  Push(x); Push(y);
  if(x->val>y->val){
    x->rs=Merge(x->rs,y); Up(x); return x;
  }
  else{
    y->ls=Merge(x,y->ls); Up(y); return y;
  }
}

void PutAns(ll x){
  if(x>=10) PutAns(x/10);
  putchar(x%10+'0');
}

typedef pair<node*,node*> Dnode;

Dnode Split(node *cur,int x){
  if(!cur) return Dnode(NULL,NULL);
  Push(cur); Dnode ret;
  if(cur->key<x){
    ret=Split(cur->rs,x);
    cur->rs=ret.fi; ret.fi=cur;
  }
  else{
    ret=Split(cur->ls,x);
    cur->ls=ret.se; ret.se=cur;
  }
  Up(cur); return ret;
}

inline node *lst(node *cur){
  while(cur->rs) cur=cur->rs;
  return cur;
}

inline node *fst(node *cur){
  while(cur->ls) cur=cur->ls;
  return cur;
}

int main(){
  t=pool;
  read(n); read(m); read(q);
  for(int i=1;i<=n;i++)
    L=Merge(L,crt(i,m,i,rand()));
  for(int i=1;i<=n;i++)
    ro[i]=Merge(ro[i],crt(i,1,1,rand()));
  while(q--){
    int x,y;
    read(x); read(y);
    if(y==m){
      Dnode A=Split(L,x),B=Split(A.se,x+1);
      PutAns(1LL*(B.fi->x-1)*m+B.fi->y); putchar('\n');
      B.fi->key=n; if(B.se) B.se->tag-=1,B.se->key-=1;
      L=Merge(A.fi,Merge(B.se,B.fi));
    }
    else{
      Dnode A=Split(ro[x],y),B=Split(A.se,y+1);
      node *cur;
      if(B.fi) cur=B.fi; else cur=lst(A.fi);
      PutAns(1LL*(cur->x-1)*m+cur->y+y-cur->key); putchar('\n');

      if((!B.se || fst(B.se)->key!=y+1) && y+1!=m)
    B.se=Merge(crt(cur->x,cur->y+y-cur->key+1,y+1,rand()),B.se);
      if(B.se) B.se->tag-=1,B.se->key-=1;
      ro[x]=Merge(A.fi,B.se);
      Dnode C=Split(L,x),D=Split(C.se,x+1);
      D.fi->key=m-1; if(D.se) D.se->tag-=1,D.se->key-=1;
      ro[x]=Merge(ro[x],D.fi);
      if(cur->key==y)
    cur->key=n,L=Merge(C.fi,Merge(D.se,cur));
      else
    L=Merge(C.fi,Merge(D.se,crt(cur->x,cur->y+y-cur->key,n,rand())));
    }
  }
  return 0;
}

分数跟去年一样490
D1T3+D2T2+D2T3=30+10+30=70
因为自己傻逼少了70分

怕是要退役了

LowestJN
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
NOIP2016滚粗记
commonc的博客
11-21 1044
DAY0 早上来之后听说老师让我们颓了 ljss今天直接没来 但是我还有一个TC的hard写了一半，决定写完了再颓.... 结果一写就写了一上午.... 中午有个高一学弟请我们吃了披萨..吃完ahcisy就回家了 然后下午一开始还是不想颓，决定打一个250分题再玩 然后就基本搞了一个下午，边看视频边写.... 晚上回家拿着手机颓了好半天，之后就睡觉了 DAY1 以为要求
NOIP2017游记
WorldWide_D的博客
11-18 1120
day0今年跑到了郊区考场…晚饭吃得比较早，所以到那里时又叫了夜宵。 晚上也没怎么腐，和平时差不多的时间就睡了。day1坐进了考场后，今年连试机的时间也没有。。坐在那干等了15min开场。待到8点半比赛开始，解压了文件后，发现今年只有example_data和题目。 没有注意事项、没有要求写什么自我介绍txt不过没有太在意。首先看完所有题。第一题竟然不会，第二题感觉就是个模拟，第三题大概是拓扑序
NOIP2017爆零记
薇小薇
11-20 1848
跟着我左手右手切割磁感线
NOIP2016 酱油记
lych的博客
11-19 1446
Day1：        赛前被连续毒奶，感觉要gg。        结果真的gg了，大家再见。。T2 lca后的O(N)怎么做啊，QAQ。 Day2：        翻盘失败，OI再见。
【游记】记NOIP2017滚粗之旅
MEET YOU FIRST TIME
11-15 725
高一滚粗记
NOIP2017官方测试数据
11-01
noip2017 、官方测试数据、包含普及组和提高组两个组别的测试数据。
NOIP2017普及组参考答案
05-02
NOIP2017普及组参考答案
1.rar（noip2017测试数据）
07-06
noip测试数据 noip2017提高组复赛
NOIP2017_day1.pdf
11-11
最新出炉NOIP2017提高组复赛day1试题，有志者速下载！
解题报告——noip 2017提高组
11-14
解题报告——noip 2017提高组,初赛提高组，pascal及其解析
NOIP2016水分记
masterwater的博客
11-23 388
Day1   拿到题目前非常紧张啊毕竟蒟蒻，貌似之前模拟题连第一题都很难水出来，赶紧奶了自己几口。     T1       果不其然第一题还是良心的水题，暴力模拟跟去年差不多。纸上手胡一下搞了半个小时就弄出来了。mogician什么的我才不认识呢一定是出题人英语不好     T2       画风骤变一脸懵逼，说好的第一天前两题是让蒟蒻安心水过拿个2=回家的呢。按自己的水平貌似只能暴力
NOIP2017 滚粗记
袋鼠与蒟蒻
11-12 479
Day 0     下午去看了看考场。。。我分配到了最菜的一个机房（ccf评测机都有2.8GHz，我们那个居然只有1.8GHz），键盘的空格键许多都有问题（跷跷板）。不过总的来说还是不错，给椅子点个赞。        晚上突然发慌（Flag），于是上洛谷疯狂水模板题，一晚上做了20多道。。。结果最后打算做一道模拟题找信心，结果从9:20调到10:30。我当时说NOIP模拟不可能这样（Flag*
NOIP2016很久后的总结和吐槽
照月之白
01-10 2316
NOIP过去了真的太久了..现在也没什么感觉了..
信息学竞赛中AC、WA、RE、CE、TLE、MLE、PE、OLE、分别什么意思
热门推荐
Annie_bing的博客
01-03 1万+
信息学竞赛中AC、WA、RE、CE、TLE、MLE、PE、OLE、分别什么意思 ‌AC Accepted 答案正确/通过 ‌WA Wrong Answer 答案错误 RE Runtime Error 运行时错误 这表明你的程序在运行过程中因为出锅而崩溃了，通常可能是访问非法内存等问题，出现这个提示但你还能过样例的话，大概率是数组没开够，仔细检查一下。 CE Complie Error 编译错误 这表明你的程序没有通过编译。如果在本地编译可以通过的话，检查一下提交语言是不是选对了或是有没有引用一些不该引用的东
ATP的NOIP2016爆零记
ADP+Pi==ATP
11-20 2192
被圆锥曲线搞到心力交瘁的ATP。。。
noip2017滚粗记
mengbi_er的博客
11-12 744
warning:这是一篇noi胸牌滚粗，noip预计不超过400的蒟蒻的游记…day -1喜闻乐见的颓废日前期…（明明就是颓废日好吗），然而我还是敲了几个模板233day 0喜闻乐见的颓废日…，望着一群会打红警的dalao，本蒟蒻瑟瑟发抖…有点虚，于是敲了敲模板…（alone_wolf:敲模板有用？），水了几道cf的…div2A…好似还有人与我讨论离散化和动态开点线段树那个快来着…day 1以为家里
NOIP2016 提高组总结
Jacky35
11-21 7885
题解每天的第二题和第三题题解已经在我的博客中放出 链接： NOIP2016提高组 第一天第二题 天天爱跑步running 题解 NOIP2016提高组 第一天第三题 换教室 classroom 题解 NOIP2016 提高组 第二天第二题 蚯蚓earthworm 题解 NOIP2016提高组 第二天第三题 愤怒的小鸟angrybirds 题解 洛谷题目链接第二次参加提高组，原来的目标是5
NOIP 2017 完挂记
დ♂Hany01's Blog Space♂ღ
11-12 465
Day0Day0Day 0 上午搞了个什么鬼动员大会，下午晚上在网上瞎逛，本想看一看splay什么的，但觉得联赛中这种东西绝对不考（立下flag），所以只是看了些联赛原题和题解，早早地就睡了。 Day1Day1Day 1 早上晕车了，感觉今天要完，刚进机房，坐到位置上，啥？没鼠标？？(差评*1)然后让监考老师从后面弄了个鼠标来，却发现电脑完全开不了（差评*2），只能换到最后一排，前面是_...
Noip2017丢人记
praying
11-12 732
...
[NOIP2017]时间复杂度
最新发布
07-27
对于一个算法的时间复杂度，可以用大O符号来表示。大O表示法描述了算法在最坏情况下的运行时间增长率。 常见的时间复杂度有： - O(1)：常数时间复杂度，表示算法的执行时间不随输入规模增加而增加。...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值