LOJ #508. 「LibreOJ NOI Round #1」失控的未来交通工具 数论 并查集

题目链接：「LibreOJ NOI Round #1」失控的未来交通工具

题目简述

一个带边权无向图，有两种操作：加边以及询问在 $x,x+b,...,x+(c−1)b$ 这些数中，有多少存在一条与之模 $m$ 同余的从 $u$ 到 $v$ 的路径（可以不是简单路径）。

数据范围

Subtask #	分值	$n, q$ 的限制 $m$ 的限制	$c$ 的限制 附加限制 1 $11$ $1 \leq n, q \leq 100$ $1 < m \leq 100$ $1 \leq c \leq 5$ 无 2 $21$ $1 \leq n, q \leq 2\times 10^5$ $m=2$ 3 $13$ $m$ 是质数
4	$7$ $m$ 是奇数	图中任何时刻都不会出现简单环
5	$10$ 无 无 6 $5$	无	$1 \leq c \leq 5$ 图中任何时刻不会有度数大于 $2$ 的点
7	$13$ 无 8 $20$ $1 \leq n, q \leq 10^6$ 无 题解 Subtask# 1 对于一开始这个数据 $n,m,q\le 100$，显然是个暴力。 但是这是个（丧）题，所以说模拟是肯定拿不到分的嘛对吧，于是就考虑怎么优化这个模拟啊，我们发现由于 $n,m$ 都很小，所以我们可以考虑每次把跑的记 $f_{i,j}$ 表示从 $u$ 开始，走到 $i$，此时 $\mod m$ 同余 $j$。这样一来所有的状态也就只有 $n*m$ 种，非常地科学。 查询直接 $\text{for}$ 一遍就可以了反正 $c$ 非常小。 所以复杂度就是 $O(cnq(n+q))$（边数也是 $q$），代码如下： （除最后全代码以外，都省略 $n,mod,q$ 的读入） namespace subtask_1{ bool f[111][111]; struct Edge{ int to,w; Edge *next; }*head[1111]; inline void add(R int u,R int v,R int w){ static Edge E[1111],*e=E; *e=(Edge){v,w,head[u]}; head[u]=e++; } int cnt; void dfs(R int u,R int p){ if(f[u][p])return; f[u][p]=1; for(R Edge *i=head[u];i;i=i->next) dfs(i->to,(p+i->w)%mod); } int main(){ for(R int opt,u,v,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); add(u,v,w);add(v,u,w); }else{ read(u);read(v);read(x);read(b);read(c); memset(f,0,sizeof f); ++cnt; dfs(u,0); ans=0; for(int i=1;i<=c;++i){ ans+=(f[v][x]); x=(x+b)%mod; } printf("%d\n",ans); } } return 0; } } 规律 然后我们要发现一个规律，就是一个非简单路径，可以分成一个简单路径和一堆环。 Subtask# 2 当 $m=2$ 时，只需要考虑奇数偶数的情况。我们还是把路径拆成一个简单路径和一些环。 对于这个简单路径上，若是有一个点连出去了一个奇环（就是可以从这个点出发，走奇数距离回到这个节点），记为点 $X$ 吧。那么从 $X$ 到 $v$ 的路径若是偶数就在这个环上绕一圈，若是奇数就在环上绕两圈凑个偶数加起来就奇数了。 如果没有这样的 $X$ 那么奇偶性就是两个点连起来的简单路径的奇偶性了。 然后用并查集维护一下联通块内是否存在奇环以及每个点到并查集根的路径长度就行了。 查询还是暴力，复杂度 $O(cq+q(\alpha n+\log(w))$ namespace subtask_2{ int id[maxn][2],fa[maxn],idcnt; int find(R int x){return x!=fa[x]?fa[x]=find(fa[x]):x;} inline void merge(int x,int y){ R int xx=find(x),yy=find(y); if(xx>yy)swap(xx,yy); if(xx!=yy)fa[xx]=yy; } int main(){ for(R int i=1;i<=n;++i){ id[i][0]=++idcnt,fa[idcnt]=idcnt; id[i][1]=++idcnt,fa[idcnt]=idcnt; } for(R int opt,u,v,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); merge(id[u][0],id[v][w&1]); merge(id[v][0],id[u][w&1]); merge(id[u][1],id[v][(w+1)&1]); merge(id[v][1],id[u][(w+1)&1]); }else{ read(u);read(v);read(x);read(b);read(c); ans=0; for(;c--;x=(x+b)&1){ ans+=(find(id[u][0])==find(id[v][x])); } printf("%d\n",ans); } } return 0; } } Subtask# 3 $m$ 是个奇质数，就可以考虑在一条边上不断绕圈。 假如这条边为 $w$，则绕一圈走了 $2w$。$m$ 又是奇质数，这样 $\gcd(2w,m)=1$。然后我们就可以用费马小定理构造出 $\equiv 1\ (\text{mod}\ m)$ 的长度，自然也就可以构造出任意长度的路径了。 那么如果 $w\equiv 0\ (\text{mod}\ m)$ 的话，显然只可以走出 $\equiv 0\ (\text{mod}\ m)$ 的路径。 所以我们还是开一个并查集维护这个并查集之中是否只存在 $\equiv 0\ (\text{mod}\ m)$ 的边就行了。 复杂度 $O(q(\alpha n+\log(w))$。代码： namespace subtask_3{ int fa[maxn]; char f[maxn]; int find(R int x){return x!=fa[x]?fa[x]=find(fa[x]):x;} inline void merge(R int x,R int y,R int w){ R int xx=find(x),yy=find(y); if(xx>yy)swap(xx,yy); if(xx!=yy)fa[yy]=xx,f[xx]|=f[yy]; if(w%mod)f[xx]=f[yy]=1; } int main(){ for(R int i=1;i<=n;++i)fa[i]=i,f[i]=0; for(R int opt,u,v,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); merge(u,v,w); }else{ read(u);read(v);read(x);read(b);read(c); ans=0; u=find(u);v=find(v); if(u!=v){ puts("0"); continue; } for(;c--;x=(x+b)%mod){ ans+=(bool)((!x)||(f[u])); } printf("%d\n",ans); } } return 0; } } Subtask# 4 当 $m$ 是一个奇数的时候，如果是质数就按照 Subtask# 3 的做法，否则的话我们考虑当前有一条边长度是 $w$，设 $\gcd(m,2w)=\gcd(m,w)=g$。 显然我们无论怎么走或者绕圈，对答案的贡献都一定是 $g$ 的倍数。推而广之，对于所有的边，显然路径的长度一定是 $G=\gcd(m,w_1,w_2,\cdots,w_{siz})$ 的倍数。 又有对于每一条边，我们可以绕出 $g$ 的任意倍数，那么就肯定可以走出 $G$ 的任意倍数。 那如何单独让一条边贡献呢？假设这条边为 $(u',v')$ 那么我们就跑 $m$ 遍路径 $(u,u')$ 即可（m是一个奇数，所以就是跑一堆来回加一次）。 因此在序列中所有是 $G$ 倍数的都满足条件。 然后同样开一个并查集维护联通块边的 $\gcd$ 就行了。 复杂度与 Subtask# 3 相同，代码由于与 Subtask# 3 以及 Subtask# 5 十分类似，请读者自行理解，就不再赘述了。 Subtask# 5 这个子任务与上一个的不同之处在于 $c$ 变得很大了，但是我们要找的东西是不变的，就是序列中存在多少个数是 $G$ 的倍数。 观察题目中给出的序列生成方式，这是一个等差数列（废话），我们记这个数列为 $\{a_n\}$。这样一来我们如果找到了第一个 $a_i$ 满足 $a_i \equiv 0\ (\text{mod}\ G)$，那么 $a_{i+1} \equiv b\ (\text{mod}\ G)$， $a_{i+j} \equiv b*j\ (\text{mod}\ G)$，则 $a_{i+G} \equiv 0\ (\text{mod}\ G)$。即从 $a_i$ 到 $a_{i+G-1}$ 构成了一个模 $G$ 意义下的完全剩余系，就是这模 $G$ 的值是存在循环的，循环节为 $G$。 这样一来我们就将题目化为了如何找到第一个 $i$ 使得 $a_i \equiv 0\ (\text{mod}\ G)$。由于 $a_i=x+b*i$，那么就是要求最小的 $i$ 使得 $x+b*i \equiv 0\ (\text{mod}\ G)$，即是求方程 $b*i\equiv -x\ (\text{mod}\ G)$ 的最小正整数解 $i$。 这东西可以化为 $b*i+G*y=G-x$，用扩展欧几里德算法求解即可。 复杂度的话 $O(q(\alpha n+\log(w)))$ 代码如下：（略去欧几里德算法以及扩展欧几里德算法的定义） namespace subtask_4_and_5{ int fa[maxn],g[maxn]; int find(R int x){return x!=fa[x]?fa[x]=find(fa[x]):x;} inline void merge(R int x,R int y,R int w){ R int xx=find(x),yy=find(y); if(xx>yy)swap(xx,yy); if(xx!=yy){ fa[yy]=xx; g[xx]=_Math_orzyrt_::gcd( _Math_orzyrt_::gcd(g[xx],g[yy]) ,w); }else{ g[xx]=_Math_orzyrt_::gcd(g[xx],w); } } int main(){ for(R int i=1;i<=n;++i)fa[i]=i,g[i]=mod; R int aa,bb,cc,dd,xx,yy; for(int opt,u,v,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); merge(u,v,w%mod); }else{ read(u);read(v);read(x);read(b);read(c); ans=0; u=find(u);v=find(v); if(u!=v){ puts("0"); continue; } aa=b%g[u];bb=g[u];cc=g[u]-x%g[u]; dd=_Math_orzyrt_::gcd(aa,bb); if(cc%dd){ puts("0"); continue; } aa/=dd;bb/=dd;cc/=dd; _Math_orzyrt_::exgcd(aa,bb,xx,yy); xx=((LL)xx*cc%bb+bb)%bb; if(xx<c)ans=(c-1-xx)/bb+1; printf("%d\n",ans); } } return 0; } } 然后Subtask# 6&7 我也不懂单独的部分分怎么做，就跳过了，因为我的分析只与 $m$ 有关。 正解 看我这个标题突然变大就知道是要讲正解了对吧。对于最后的这些部分分，$m$ 都是没有限制的，就没办法表示成某个数的整数倍了，怎么办呢？ 按照套路，我们把路径拆成一个简单路径和一堆环，显然对于每一个环，假如其环长为 $W$，那么就可以在环上走出 $g=\gcd(W,m)$ 的任意倍数的贡献。 怎么加入一个环呢？直接从开头到环上一点跑 $m$ 个来回就行了。 于是我们就可以仿照 $m$ 为奇数时的构造方法，使得最终答案为 $G=\gcd(m,W_1,W_2,\cdots,W_{siz})$ 的倍数加上一条简单路径长度。 于是用并查集维护连通块内的 $G$ 以及每个点到并查集的根路径长，回答询问用扩展欧几里得算法。 那么怎么用并查集维护联通块内环的 $\gcd$ 呢？在加入连通块间的边时肯定是直接合并的咯没啥问题；加入连通块内的边（设其为 $(u',v’)$，长为$w$）时，设原来环的 $\gcd$ 为 $G$，$u'$ 到 $v'$ 某路径长为 $L$，新加入的边，则新的 $G'=\gcd(G,2w,w+L)$。 复杂度 $O(q(\alpha n+\log(w)))$。 本题完整代码如下： #include <bits/stdc++.h> #define LL long long #define R register template<class TT>inline TT Max(R TT a,R TT b){return a<b?b:a;} template<class TT>inline TT Min(R TT a,R TT b){return a<b?a:b;} using namespace std; char __B[1<<18],*__S=__B,*__T=__B; #define getchar() (__S==__T&&(__T=(__S=__B)+fread(__B,1,1<<15,stdin),__S==__T)?EOF:*__S++) template<class TT>inline void read(R TT &x){ x=0;R bool f=0;R char ch=getchar(); for(;ch<48||ch>57;ch=getchar())f|=(ch=='-'); for(;ch>47&&ch<58;ch=getchar()) x=(x<<1)+(x<<3)+(ch^48); (f)&&(x=-x); } template<class TT> char write(R TT x){ x<0&&(putchar('-'),x=-x); x>9&&(write(x/10)); putchar(x%10+48); } int n,q,mod; namespace _Math_orzyrt_{ template<class TT> inline bool isprime(R TT &x){ for(R TT i=2;i*i<=x;++i) if(x%i==0)return 0; return 1; } template<class TT> TT gcd(R TT a,R TT b){return b?gcd(b,a%b):a;} template<class TT> void exgcd(R TT a,R TT b,R TT &x,R TT &y){ b?(exgcd(b,a%b,y,x),y-=a/b*x):(x=1,y=0); } } namespace subtask_1{ bool f[111][111]; struct Edge{ int to,w; Edge *next; }*head[1111]; inline void add(R int u,R int v,R int w){ static Edge E[1111],*e=E; *e=(Edge){v,w,head[u]}; head[u]=e++; } int cnt; void dfs(R int u,R int p){ if(f[u][p])return; f[u][p]=1; for(R Edge *i=head[u];i;i=i->next) dfs(i->to,(p+i->w)%mod); } int main(){ for(R int opt,u,v,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); add(u,v,w);add(v,u,w); }else{ read(u);read(v);read(x);read(b);read(c); memset(f,0,sizeof f); ++cnt; dfs(u,0); ans=0; for(int i=1;i<=c;++i){ ans+=(f[v][x]); x=(x+b)%mod; } printf("%d\n",ans); } } return 0; } } const int maxn = 1e6; namespace subtask_2{ int id[maxn][2],fa[maxn],idcnt; int find(R int x){return x!=fa[x]?fa[x]=find(fa[x]):x;} inline void merge(int x,int y){ R int xx=find(x),yy=find(y); if(xx>yy)swap(xx,yy); if(xx!=yy)fa[xx]=yy; } int main(){ for(R int i=1;i<=n;++i){ id[i][0]=++idcnt,fa[idcnt]=idcnt; id[i][1]=++idcnt,fa[idcnt]=idcnt; } for(R int opt,u,v,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); merge(id[u][0],id[v][w&1]); merge(id[v][0],id[u][w&1]); merge(id[u][1],id[v][(w+1)&1]); merge(id[v][1],id[u][(w+1)&1]); }else{ read(u);read(v);read(x);read(b);read(c); ans=0; for(;c--;x=(x+b)&1){ ans+=(find(id[u][0])==find(id[v][x])); } printf("%d\n",ans); } } return 0; } } namespace subtask_3{ int fa[maxn]; char f[maxn]; int find(R int x){return x!=fa[x]?fa[x]=find(fa[x]):x;} inline void merge(R int x,R int y,R int w){ R int xx=find(x),yy=find(y); if(xx>yy)swap(xx,yy); if(xx!=yy)fa[yy]=xx,f[xx]|=f[yy]; if(w%mod)f[xx]=f[yy]=1; } int main(){ for(R int i=1;i<=n;++i)fa[i]=i,f[i]=0; for(R int opt,u,v,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); merge(u,v,w); }else{ read(u);read(v);read(x);read(b);read(c); ans=0; u=find(u);v=find(v); if(u!=v){ puts("0"); continue; } for(;c--;x=(x+b)%mod){ ans+=(bool)((!x)||(f[u])); } printf("%d\n",ans); } } return 0; } } namespace subtask_4_and_5{ int fa[maxn],g[maxn]; int find(R int x){return x!=fa[x]?fa[x]=find(fa[x]):x;} inline void merge(R int x,R int y,R int w){ R int xx=find(x),yy=find(y); if(xx>yy)swap(xx,yy); if(xx!=yy){ fa[yy]=xx; g[xx]=_Math_orzyrt_::gcd( _Math_orzyrt_::gcd(g[xx],g[yy]) ,w); }else{ g[xx]=_Math_orzyrt_::gcd(g[xx],w); } } int main(){ for(R int i=1;i<=n;++i)fa[i]=i,g[i]=mod; R int aa,bb,cc,dd,xx,yy; for(int opt,u,v,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); merge(u,v,w%mod); }else{ read(u);read(v);read(x);read(b);read(c); ans=0; u=find(u);v=find(v); if(u!=v){ puts("0"); continue; } aa=b%g[u];bb=g[u];cc=g[u]-x%g[u]; dd=_Math_orzyrt_::gcd(aa,bb); if(cc%dd){ puts("0"); continue; } aa/=dd;bb/=dd;cc/=dd; _Math_orzyrt_::exgcd(aa,bb,xx,yy); xx=((LL)xx*cc%bb+bb)%bb; if(xx<c)ans=(c-1-xx)/bb+1; printf("%d\n",ans); } } return 0; } } namespace final{ int fa[maxn],g[maxn],d[maxn]; int find(R int x){ if(x!=fa[x]){ R int f=fa[x]; fa[x]=find(fa[x]); d[x]=(d[x]+d[f])%mod; } return fa[x]; } inline void merge(R int x,R int y,R int w){ R int xx=find(x),yy=find(y); if(xx>yy)swap(xx,yy); if(xx!=yy){ fa[yy]=xx; d[yy]=((unsigned)d[x]+d[y]+w)%mod; g[xx]=_Math_orzyrt_::gcd( _Math_orzyrt_::gcd(g[xx],g[yy]) ,w<<1); }else{ g[xx]=_Math_orzyrt_::gcd( _Math_orzyrt_::gcd((unsigned)w<<1,(unsigned)d[x]+d[y]+w) ,(unsigned)g[xx]); } } int main(){ for(R int i=1;i<=n;++i)fa[i]=i,g[i]=mod,d[i]=0; R int aa,bb,cc,dd,xx,yy; for(R int opt,u,v,uu,vv,w,x,b,c,ans;q--;){ read(opt); if(opt==1){ read(u);read(v);read(w); merge(u,v,w%mod); }else{ read(u);read(v);read(x);read(b);read(c); ans=0; uu=find(u);vv=find(v); if(uu!=vv){ puts("0"); continue; } aa=b%g[uu];bb=g[uu];cc=d[u]+d[v]-x; dd=_Math_orzyrt_::gcd(aa,bb); if(cc%dd){ puts("0"); continue; } aa/=dd;bb/=dd;cc/=dd; _Math_orzyrt_::exgcd(aa,bb,xx,yy); xx=((LL)xx*cc%bb+bb)%bb; if(xx<c)ans=(c-1-xx)/bb+1; write(ans); putchar('\n'); } } return 0; } } int main(){ read(n);read(mod);read(q); if(mod==2)return subtask_2::main(); if(n<111&&111>q)return subtask_1::main(); if(mod==2)return subtask_2::main(); if(_Math_orzyrt_::isprime(mod))return subtask_3::main(); if(mod&1)return subtask_4_and_5::main(); return final::main(); } 确定要放弃本次机会？ 福利倒计时 : : 立减 ¥ 普通VIP年卡可用 立即使用 Night2002 关注 关注 2 点赞 踩 0 收藏 觉得还不错? 一键收藏 知道了 0 评论 复制链接 分享到 QQ 分享到新浪微博 扫一扫 专栏目录 up_717085_google-api-php-client_71loj.rar 12-22 1. **OAuth2身份验证**：Google API通常需要通过OAuth2协议进行身份验证，确保应用程序能安全地访问用户的资源。开发者需要理解如何获取和使用OAuth2的客户端ID和秘密。 2. **Google API Client Library for PHP的... LOJ#508「LibreOJ NOI Round #1」失控的未来交通工具 Log_x's Blog 02-27 340 Address LOJ#508 Algorithm 1 设 fx,if_{x,i}fx,i​ 表示是否存在终点在点 xxx 且权值和模 mmm 的值为 iii 的路径，对每组询问记忆化搜索即可，时间复杂度 O(nmq)\mathcal O(nmq)O(nmq)，期望得分 11pts\text{11pts}11pts。 Algorithm 2 针对 m=2m = 2m=2 的数据。 询问时 uu... 参与评论 您还未登录，请先 登录 后发表或查看评论 LibreOJ NOI Round #1 Day 1 B. 失控的未来交通工具 weixin_30294295的博客 07-08 110 瞬间移动 官方题解 题意：一个带边权无向图，加边以及询问在x,x+b,...,x+(c−1)bx,x+b,...,x+(c-1)bx,x+b,...,x+(c−1)b这些数中，有多少存在一条与之模m同余的从u到v的路径（可以不是简单路径）。 考场上读错题系列，以为边是有向的，然后就完全不可做了对不对…… 由于是无向边，而且路径可以不是简单路径，那就意味着我们可以... LOJ508失控的未来交通工具(图论绕环gcd) xiaxiaoguang_的博客 07-07 152 模数为奇数的时候,我们可以在某个环上绕圈 此时单独一个边是可以被记入方案里的,构造方法就是我们在一开始在u想要到达v,然后中间有(x,y)(x,y)(x,y)一条边,那么我们只需要u,xu,xu,x绕m次,最后一定停在x,然后走一步,然后y,vy,vy,v绕m次走到v即可 那么我们考虑此时最优的一定是这些边,不需要绕任何一个环,我们所有边方案组合就能得到最优 而想要所有边我们只需要找到任意一个经过所有边若干次的路径就好了,然后再mod  m\mod mmodm意义下假设我们有一个G=gcd⁡(x1..xm, #508. 「LibreOJ NOI Round #1」失控的未来交通工具（带权并查集+exgcd） Freopen的博客 11-01 297 Portal 不会写带权并查集+1 tutorial 额， 首先可以发现像树上最短异或路径一样， 对于每个联通块的求一个生成树， 对于每个环我们可以在走(u,v)(u,v)(u,v)的树上路径的同时逃到这个环上走任意圈然后走回来。 而中间的路我们只需要重复走mmm趟来回就可以看作没有贡献。 然后在(modm)\pmod m(modm)的情况下，有一个长度为kkk的环，我们就可以达到:树上路径长+p... [数论] LOJ #508. 「LibreOJ NOI Round #1」失控的未来交通工具 雯舞 07-23 1138 这种非简单路 一般转化成任意一条路加上若干环 这里大概是任意一条路加上若干环长的gcd 任意一条路 可以弄出任意一颗生成树？ 但是这里实际上只需要一个带权并查集 详见官方题解#include<cstdio> #include<cstdlib> #include<algorithm> using namespace std; typedef long long ll; typedef pair< Loj X:创建逻辑电路-开源 07-04 Loj X 是一个开源项目，专注于提供一个平台，让人们能够创建和设计逻辑电路，从而实现基本的计算功能，甚至构建小型的完整处理器。这个工具集成了各种元件，如逻辑门、存储器和程序读取器，使得用户可以构建出复杂的... log4j配置，jar包 06-01 总结来说，Log4j是Java开发中不可或缺的工具，通过配置`log4j.properties`文件，我们可以精确地控制日志的输出，便于调试和监控系统运行状态。而`log4j-1.2.11.jar`则包含了实现这些功能所需的类库。 scicloj.github.com：Scicloj网站 02-08 scicloj.github.com：Scicloj网站 动态规划_方泓杰.pdf 09-14 这个概念由美国数学家Richard Bellman在20世纪50年代提出，其核心在于将一个复杂问题分解为若干个子问题，并通过解决这些子问题来找到原问题的最优解。 动态规划的基本原理是利用重叠子问题和最优子结构两个关键... 【LOJ NOI Round#2 Day1 T1】单枪匹马【矩阵】 gigo_64 sink in my world 10-11 146 传送门 支持后面加数，区间进行奇怪的要求操作算值。 你看这个操作是个线性的迭代，这种时候就应该想到线性代数，那就想到矩阵。 你就强行用矩阵表示，然后手玩一下结果矩阵来反推原来的矩阵。 手玩出来每个数的初始矩阵是 0 1 1 a[i] 其中右下角是分子，右上角是分母。 正确性？找两个矩阵乘起来就正确了。 那么一个区间可以视作1到r的矩阵乘去掉1到l-1。 这个去掉操作联想乘逆元，... [反演] LOJ #509. 「LibreOJ NOI Round #1」动态几何问题 雯舞 07-23 1251 μ2(n)=∑d2|nμ(d)\mu^2(n)=\sum_{d^2|n} \mu(d) 然后就是xjb推 反正退役了 我也就弃坑了 95分代码 复杂度分析及优化详见官方题解#include<cstdio> #include<cstdlib> #include<cmath> #include<algorithm> using namespace std; typedef long long l [数论] LOJ #510. 「LibreOJ NOI Round #1」北校门外的回忆 雯舞 07-11 1031 这个题跟树状数组没有半毛钱关系 首先这是一个最低位翻倍的过程，如果这一位最终会变成 00 ，那么步数是 O(logn)O(\log n)的？ 要是这一位不能变成 00 ，也就是在环上跑了，似乎跑到环上的步数也是 O(logn)O(\log n)的？ 然后就变成了，修改要是跑到 00 就暴力跑，跑到环上，询问照样 设某个修改最早到环上的起点是w+z×kmw+z\times k^m，那么会对 w 507 LOJ 「LibreOJ NOI Round #1」接竹竿 weixin_30566063的博客 02-10 203 思路：这个题是在一个n2的复杂度之上进行了dp方程的优化变形，最后变成线性dp，但还是有一些不懂（菜是本质，%oi爷们） 代码： #include <bits/stdc++.h> using namespace std; const int maxn=1e6+5; typedef long long LL; int n,k; int col[maxn]; LL dp... 【LOJ3059】「HNOI2019」序列 cz_xuyixuan的博客 04-24 1124 【题目链接】 点击打开链接 【思路要点】 没有修改的做法在《IOI2018中国国家候选队论文集——浅谈保序回归问题》中有所介绍。 具体做法如下： (1)(1)(1) 、注意到若所有 BiB_iBi​ 均相等，最小化 ∑i=1N(Ai−B)2=∑i=1NAi2−2AiB+B2\sum_{i=1}^{N}(A_i-B)^2=\sum_{i=1}^{N}A_i^2-2A_iB+B^2∑i=1... LOJ #539. 「LibreOJ NOIP Round #1」旅游路线 倍增floyd + 思维 EM LGH 08-20 149 考试的时候是这么想的： 求出每一个点花掉 $i$ 的花费向其他点尽可能走的最长距离，然后二分这个花费，找到第一个大于 $d$ 的就输出$.$然而，我这个记忆化搜索 $TLE$ 的很惨$.$这里讲一下正解： 上面的大题思路是正确的，但是记忆化搜索太慢，考虑倍增 $floyd.$令 $f[i][j]$ 表示 $i$ 号点花费 $j$ 能走的最远距离$.$令 $go[i][j][k]$ 表示 ... loj 1002(spfa变形） weixin_33806300的博客 10-11 75 题目链接：http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=25828 题意：求所有点到给定的目标顶点的路径上的权值的最大值的最小值。 思路：spfa的应用，更新的时候判断max(dist[u],w(u,v))&lt;dist[v]即可。 1 #include&lt;iostream&gt; 2 #inclu... LOJ P507 「LibreOJ NOI Round #1」接竹竿【动态规划】 Yucohny 06-15 674 题目分析： 这仍然是一道显然的DP题。 首先由于取出一段区间的要求是这段区间的两个端点是一样的，所以我们可以考虑用类似于链式前向星的方法储存相同端点的到达。 由数据范围，考虑定一个一维的状态，一维的状态始终是十分好定的： DP[I]表示前I张扑克处理后能够得到的最大得分 接下来考虑一下转移： 对于第I张照片，有可以取出计分与不去出计分两种情况。 若第I张照片不取出计分，则DP[I]... LOJ10172涂抹果酱 zhangjingyanzjyer的博客 10-18 387 https://loj.ac/problem/10172 题目大意就是给了你n行m列的矩阵，让你给这个矩阵的每个格子染色，有1,2,3三种颜色可以选择。同时，必须满足任意两个相邻的格子的颜色不能相同。并且第k行已经完成了染色，无法被更改。求合法的染色方案数。 这道题n的范围很大，但m的范围很小，所以很容易想到可以使用状压dp来求解。由于每一个点的状态有三种，所以要用三进制数来表示每一行的状态。... LOJ6354 & 洛谷4366：[Code+#4]最短路——题解 最新发布 05-17 算法的基本思想是：从源点开始，每次选择当前最短路径的节点，并更新其周围节点的距离。 具体步骤如下： 1. 初始化距离数组 dist，将源点到自己的距离设为 0，其他点到源点的距离设为无穷大。 2. 创建一个集合 ... Night2002 CSDN认证博客专家 CSDN认证企业博客 码龄7年 暂无认证 12 原创 42万+ 周排名 65万+ 总排名 5862 访问 等级 207 积分 7 粉丝 12 获赞 3 评论 6 收藏 私信 关注 热门文章 Night的数据结构杂谈-虚树 799 LOJ #508. 「LibreOJ NOI Round #1」失控的未来交通工具 数论 并查集 666 关于一类积性函数求和的优化——杜教筛 621 Night关于数学的杂谈-插值法 573 Night的算法杂谈-2SAT问题 558 分类专栏 比赛总结 1篇 算法与数据结构 6篇 题解们 5篇 最新评论 Night的算法杂谈-2SAT问题 Steaunk: 您这个逻辑符号写反了吧？应该是与不是或吧？ 记2018.2.27 NOIP班-数据结构(9) Night2002 回复 zxy_good: 看oj啊 Night的数据结构杂谈-虚树 zxy_good: 学习了，谢谢！！！ 大家在看 【人工智能时代】- 机器学习开源框架大总结 Cocos2d-x：Cocos2d-x性能优化技巧 python获取同花顺涨停连板跌停数据代码可用 531 问：Java中抽象类和接口的差异，真心了解不？ Postman：Postman API测试：Postman中创建第一个API请求 1455 最新文章 关于一类积性函数求和的优化——杜教筛 Night的数学杂谈——模数不互质的CRT [bzoj3622]已经没有什么好害怕的了 动态规划 2018年12篇 目录 目录 分类专栏 比赛总结 1篇 算法与数据结构 6篇 题解们 5篇 目录 评论 被折叠的  条评论 为什么被折叠? 到【灌水乐园】发言 查看更多评论 添加红包 祝福语 请填写红包祝福语或标题 红包数量 个 红包个数最小为10个 红包总金额 元 红包金额最低5元 余额支付 当前余额3.43元 前往充值 > 需支付：10.00元 取消 确定 下一步 知道了 成就一亿技术人! 领取后你会自动成为博主和红包主的粉丝 规则 hope_wisdom 发出的红包 实付元 使用余额支付 点击重新获取 扫码支付 钱包余额 0 抵扣说明： 1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。 2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。 余额充值