因为一些翻译问题铂金只过了 B 。之后看只剩余大约 1h 于是弃疗了。
但是某沙雕同学提供的假翻译似乎也挺能做而且还算有趣,以后可能会把它出出来?也算不错了,膜拜带翻译家
不过并不能改变最后排名垫底了。
后来想想感觉 A 还是很能做的
C 是真的不会,考后找某集训队学长请教,结果只学会了部分分的做法。。。
以上是废话
Bronze
三个普及组题。按题意模拟即可
代码懒得放了。。。
Silver
A
求第 (k) 个既不被 (3) 整除,也不被 (5) 整除的数字
(kle 10^9)
发现 (15) 个形成一个循环,边角直接加上去就行了
#include <bits/stdc++.h>
using namespace std;
int l[]={0,1,2,4,7,8,11,13,14};
int main (){
freopen ("moobuzz.in","r",stdin);
freopen ("moobuzz.out","w",stdout);
int n;scanf ("%d",&n);
int m=(n-1)/8;
printf ("%d",m*15+l[(n-1)%8+1]);
return 0;
}
B
长度为 (L) 的数轴上有 (n) 个点,每个点有权值和左/右的方向,每个时刻一个单位的移动。若两个点碰撞则交换方向(并不交换权值)后继续运动,若一个点到了 (0) 或 (L) 则停止运动。问停止运动的权值之和恰好大于 (frac{sum w_i}{2}) 的时候发生了多少次碰撞。
(Lle 10^9,nle 5 imes 10^4,w_ile 10^3)
挺有意思的一个题,这一类题一般都是看似交换,实则当成不交换处理。一开始以为这题也可以直接无视交换方向,后来发现权值这种东西并不能直接这样搞。。。
考虑一个向右的和一串向左的连续碰撞之后的结果,应该是等价于把向右的那个反向,一串向左的最后一个反向,剩余的不受影响。
推广一下这个结论,假设不考虑碰撞,有 (x) 个点可以到达 (0) ,那么这 (x) 个点的权值之和一定是从左往右的前 (x) 个之和,证明比较显然。到 (L) 的也同理。
所以可以二分答案找出这个时刻,后面求碰撞次数就可以无视交换权值方向直接做了,具体实现可以用二分。
不知道为什么出题人放 (n^2) 过了
#include <bits/stdc++.h>
using namespace std;
const int N=5e4+5;
int n,L,sum=0;
int w[N],x[N],d[N];
struct Node{
int pos,w;
}a[N];
bool check(int md){
int S=0;int c1=0,c2=0;
for (int i=1;i<=n;i++)
if (x[i]+d[i]*md<=0) c1++;
else if (x[i]+d[i]*md>=L) c2++;
for (int i=1;i<=c1;i++) S+=a[i].w;
for (int i=1;i<=c2;i++) S+=a[n-i+1].w;
return S*2>=sum;
}
bool cmp(Node a,Node b){
return a.pos<b.pos;
}
int p1[N],p2[N];
int main (){
freopen ("meetings.in","r",stdin);
freopen ("meetings.out","w",stdout);
scanf ("%d%d",&n,&L);
int c1=0,c2=0;
for (int i=1;i<=n;i++){
scanf ("%d%d%d",&w[i],&x[i],&d[i]);sum+=w[i];
if (d[i]==1) p1[++c1]=x[i];
else p2[++c2]=x[i];
a[i]=(Node){x[i],w[i]};
}
sort(a+1,a+n+1,cmp);
sort(p1+1,p1+c1+1);
sort(p2+1,p2+c2+1);
int l=0,r=L,T;
while (l<=r){
int mid=(l+r)>>1;
if (check(mid)) T=mid,r=mid-1;
else l=mid+1;
}
long long ans=0;
for (int i=1;i<=c1;i++){
l=lower_bound(p2+1,p2+c2+1,p1[i])-p2;
r=upper_bound(p2+1,p2+c2+1,p1[i]+2*T)-p2-1;
ans+=r-l+1;
}
printf ("%lld",ans);
return 0;
}
C
给一棵树点权为 (0/1),(q) 次询问路径 ((u,v)) 上有没有 (0/1)
(nle 10^5,qle 10^5)
一个挺无聊的题,感觉已经被出烂了。可以直接倍增维护是否存在,记向上 (2^k) 步有没有 (0/1) ,求 LCA 的时候顺便搞出来就行了
#include <bits/stdc++.h>
using namespace std;
const int N=1e5+5;
int fa[20][N];
bool G[20][N],H[20][N];
char s[N];
int Head[N],Next[N<<1],Adj[N<<1],tot=0;
void addedge(int u,int v){
Next[++tot]=Head[u],Head[u]=tot,Adj[tot]=v;
Next[++tot]=Head[v],Head[v]=tot,Adj[tot]=u;
}
int deep[N];
void dfs(int x,int f){
fa[0][x]=f;
for (int e=Head[x];e;e=Next[e])
if (Adj[e]!=f) deep[Adj[e]]=deep[x]+1,dfs(Adj[e],x);
}
int main (){
freopen ("milkvisits.in","r",stdin);
freopen ("milkvisits.out","w",stdout);
int n,m;cin>>n>>m;
cin>>(s+1);
for (int i=1;i<=n;i++)
if (s[i]=='H') H[0][i]=1;
else G[0][i]=1;
for (int i=1,u,v;i<n;i++) {
cin>>u>>v;
addedge(u,v);
}
deep[1]=1;
dfs(1,0);
for (int j=1;j<=18;j++)
for (int i=1;i<=n;i++){
fa[j][i]=fa[j-1][fa[j-1][i]];
G[j][i]=G[j-1][i]|G[j-1][fa[j-1][i]];
H[j][i]=H[j-1][i]|H[j-1][fa[j-1][i]];
}
while (m--){
int u,v;char c;cin>>u>>v>>c;
if (deep[u]<deep[v]) swap(u,v);
if (c=='G') {
bool ans=false;
for (int j=18;j>=0;j--)
if (deep[fa[j][u]]>=deep[v])
ans|=G[j][u],u=fa[j][u];
ans|=G[0][u];
if (u!=v){
for (int j=18;j>=0;j--)
if (fa[j][u]!=fa[j][v])
ans|=G[j][u]|G[j][v],u=fa[j][u],v=fa[j][v];
ans|=G[0][fa[0][u]],ans|=G[0][u],ans|=G[0][v];
}
printf ("%d",ans);
}else{
bool ans=false;
for (int j=18;j>=0;j--)
if (deep[fa[j][u]]>=deep[v])
ans|=H[j][u],u=fa[j][u];
ans|=H[0][u];
if (u!=v){
for (int j=18;j>=0;j--)
if (fa[j][u]!=fa[j][v])
ans|=H[j][u]|H[j][v],u=fa[j][u],v=fa[j][v];
ans|=H[0][fa[0][u]],ans|=H[0][u],ans|=H[0][v];
}
printf ("%d",ans);
}
}
return 0;
}
Gold
A
给一张图,每条边有流量限制和边权,选一条 (1) 到 (N) 路径,使得(frac{min{f_i}}{sum w_i}) 最大,输出这个值
(nle 10^3,mle 10^3,f_ile 10^3)
没想太多。。看到流量很小,可以枚举流量,然后在原图上跑 dijkstra,复杂度 (O(fmlog m)) 。不知道有没有高论
#include <bits/stdc++.h>
using namespace std;
const int N=1005,E=2005;
struct Edge{
int u,v,f,w;
}e[E];
int n,m;
int dis[N];bool vis[N];
int Head[N],Next[E],Adj[E],Weight[E],tot=0;
void addedge(int u,int v,int w){
Next[++tot]=Head[u],Head[u]=tot,Adj[tot]=v,Weight[tot]=w;
Next[++tot]=Head[v],Head[v]=tot,Adj[tot]=u,Weight[tot]=w;
}
priority_queue < pair <int,int> > Q;
int dijkstra(int lim){
memset (vis,0,sizeof(vis));
memset (dis,0x3f,sizeof(dis));
memset (Head,0,sizeof(Head));
memset (Next,0,sizeof(Next));
memset (Adj,0,sizeof(Adj));
memset (Weight,0,sizeof(Weight));tot=0;
Q.push(make_pair(0,1));
dis[1]=0;
for (int i=1;i<=m;i++)
if (e[i].f>=lim) addedge(e[i].u,e[i].v,e[i].w);
while (!Q.empty()){
int x=Q.top().second;Q.pop();
if (vis[x]) continue;vis[x]=true;
for (int e=Head[x];e;e=Next[e])
if (dis[x]+Weight[e]<dis[Adj[e]]){
dis[Adj[e]]=dis[x]+Weight[e];
Q.push(make_pair(-dis[Adj[e]],Adj[e]));
}
}
return dis[n];
}
int main (){
freopen ("pump.in","r",stdin);
freopen ("pump.out","w",stdout);
scanf ("%d%d",&n,&m);
for (int i=1;i<=m;i++) scanf ("%d%d%d%d",&e[i].u,&e[i].v,&e[i].w,&e[i].f);
double ans=0;
for (int i=1;i<=1000;i++)
ans=max(ans,(double)i/(double)dijkstra(i));
printf ("%lld",(long long)(ans*1000000));
return 0;
}
B
银组 t3 加强版,把两种属性变成了 (10^5) 种
(nle 10^5,mle 10^5)
比赛时写了个很蠢的 (nlog^2n) 的做法,因为懒得重写直接把以前代码拷来改改就行了,用 (c) 棵动态开点线段树维护所有颜色,然后直接树剖在对应线段树上看有没有点即可
另一个很好写的单 (log) 做法:树上差分,用主席树维护每个节点到根每种颜色出现次数,查询可以用两个点的数量减去 LCA 和 LCA父亲 处的数量,判断是否为 (0) 即可
因为没有强制在线,或许能有更低的复杂度?
#include <bits/stdc++.h>
using namespace std;
const int N=100005;
int Head[N],Next[N<<1],Adj[N<<1],tot=0;
void addedge(int u,int v){
Next[++tot]=Head[u],Head[u]=tot,Adj[tot]=v;
Next[++tot]=Head[v],Head[v]=tot,Adj[tot]=u;
}
int fa[N],deep[N],son[N],size[N],top[N];
void dfs(int x,int f){
size[x]=1;
for (int e=Head[x];e;e=Next[e])
if (Adj[e]!=f){
fa[Adj[e]]=x;
deep[Adj[e]]=deep[x]+1;
dfs(Adj[e],x);
size[x]+=size[Adj[e]];
son[x]=(size[son[x]]<size[Adj[e]]?Adj[e]:son[x]);
}
}
int dfn[N],to[N],Time=0;
void dfs2(int x,int tp){
top[x]=tp;dfn[x]=++Time,to[Time]=x;
if (!son[x]) return;dfs2(son[x],tp);
for (int e=Head[x];e;e=Next[e])
if (Adj[e]!=fa[x]&&Adj[e]!=son[x])
dfs2(Adj[e],Adj[e]);
}
int rt[N];
int ls[N*30],rs[N*30];
int cnt=0;
#define mid ((l+r)>>1)
void insert(int&root,int l,int r,int x){
if (!root) root=++cnt;
if (l==r) return;
if (x<=mid) insert(ls[root],l,mid,x);
else insert(rs[root],mid+1,r,x);
}
bool query(int root,int l,int r,int L,int R){
if ((!root)||r<L||l>R) return false;
if (L<=l&&r<=R) return true;
return query(ls[root],l,mid,L,R)||query(rs[root],mid+1,r,L,R);
}
int s[N];
int main (){
freopen ("milkvisits.in","r",stdin);
freopen ("milkvisits.out","w",stdout);
int n,m;scanf ("%d%d",&n,&m);
for (int i=1;i<=n;i++) scanf ("%d",&s[i]);
for (int i=1;i<n;i++){
int u,v;scanf ("%d%d",&u,&v);
addedge(u,v);
}
dfs(1,0),dfs2(1,1);
for (int i=1;i<=n;i++) insert(rt[s[i]],1,n,dfn[i]);
while (m--){
int u,v,w;scanf ("%d%d%d",&u,&v,&w);
bool ans=false;
while (top[u]!=top[v]){
if (deep[top[u]]<deep[top[v]]) swap(u,v);
ans|=query(rt[w],1,n,dfn[top[u]],dfn[u]);
u=fa[top[u]];
}
if (dfn[u]>dfn[v]) swap(u,v);
ans|=query(rt[w],1,n,dfn[u],dfn[v]);
printf ("%d",ans);
}
return 0;
}
C
给定一个字符串和一个 (M imes M) 的矩阵 (a),表示可以消耗 (a_{i,j}) 的代价把字符 (i) 改成字符 (j),求把原字符串改成每个连续段长度都不小于 (K) 的最小代价。
(n,Kle 10^5,Mle 26)
首先直接用原矩阵的交换方法不一定最优,拿到手先跑一遍 floyd。
然后考虑 dp ,(f_i) 表示前 (i) 个中每个连续段都 (ge k) 的最小代价,显然有转移
其中 (cost(i,j,c)) 表示把 ([i,j]) 全都改成字符 (c) 的最小代价,这个可以前缀和预处理。
观察转移,发现可以在枚举到 (i) 时再把 (i-k) 的决策加进去,这样维护一个前缀 (min) 即可做到 (O(1)) 转移。
复杂度 (O(n imes M))
#include <bits/stdc++.h>
using namespace std;
const int N=100005;
int c[30][30],f[N];
char s[N];
int sum[30][N];
int query(int x,int l,int r){
return sum[x][r]-sum[x][l-1];
}
int mn[30];
int main (){
freopen ("cowmbat.in","r",stdin);
freopen ("cowmbat.out","w",stdout);
int n,m,k;scanf ("%d%d%d",&n,&m,&k);
scanf ("%s",s+1);
for (int i=0;i<m;i++)
for (int j=0;j<m;j++)
scanf ("%d",&c[i][j]);
for (int l=0;l<m;l++)
for (int i=0;i<m;i++)
for (int j=0;j<m;j++)
if (i!=j&&j!=l&&i!=l)
c[i][j]=min(c[i][j],c[i][l]+c[l][j]);
for (int i=0;i<m;i++)
for (int j=1;j<=n;j++)
sum[i][j]=c[s[j]-'a'][i]+sum[i][j-1];
memset (mn,0x3f,sizeof(mn));
memset (f,0x3f,sizeof(f));f[0]=0;
for (int i=k;i<=n;i++)
for (int j=0;j<m;j++)
mn[j]=min(mn[j]+c[s[i]-'a'][j],f[i-k]+query(j,i-k+1,i)),f[i]=min(f[i],mn[j]);
printf ("%d",f[n]);
return 0;
}
Platinum
B
给定一棵以 (1) 为根节点的树,支持以下两种操作:
- 给定 (x,c) ,给 (x) 的子树中每个点的属性集合里插入一种属性 (c)
- 给定 (x) ,求 (x) 的子树内所有点各有多少种不同的属性,需要输出属性个数之和
(nle 10^5,qle 10^5)
dfs 序上子树表示为一段连续的区间,因此两个子树要么包含,要么不相交
考虑修改一个点,影响的范围是 (x) 的子树,我们给这个范围内所有点的权值都 +1
但是这样会产生冲突,如果一个属性已经有过了,那么就会有重复,为了消除重复下面分两种情况:
- 考虑当前操作的区间被某个区间包含过,那么这是个废区间,直接不做
- 如果这个区间包含之前的某些区间,那么撤销那些区间的影响,再加上这个区间
实现方面,可以用 set 维护情况 2 并支持撤销,操作 1 我写了个很蠢的线段树,现在想想大概是不必要的
这样,查询操作就可以做一次线段树区间求和
时间复杂度 (O(qlog n))
代码里有很多因为假题意留下的无用的东西,并不推荐食用
#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int N=100005,mx=100000;
int Head[N],Next[N<<1],Adj[N<<1],tot=0;
void addedge(int u,int v){
Next[++tot]=Head[u],Head[u]=tot,Adj[tot]=v;
Next[++tot]=Head[v],Head[v]=tot,Adj[tot]=u;
}
set < pair<int,int> > S[N];
set <int> P[N];
int dfn[N],to[N],Time=0;
int deep[N],size[N],son[N],top[N],fa[N];
void dfs(int x,int f){
size[x]=1,fa[x]=f;
for (int e=Head[x];e;e=Next[e])
if (Adj[e]!=f){
deep[Adj[e]]=deep[x]+1;
dfs(Adj[e],x);
size[x]+=size[Adj[e]];
son[x]=(size[son[x]]>size[Adj[e]]?son[x]:Adj[e]);
}
}
void dfs2(int x,int tp){
top[x]=tp;dfn[x]=++Time,to[Time]=x;
if (!son[x]) return;
dfs2(son[x],tp);
for (int e=Head[x];e;e=Next[e])
if (Adj[e]!=fa[x]&&Adj[e]!=son[x])
dfs2(Adj[e],Adj[e]);
}
inline int LCA(int u,int v){
if ((!u)||(!v)) return 0;
for (;top[u]!=top[v];deep[top[u]]<deep[top[v]]?v=fa[top[v]]:u=fa[top[u]]);
return deep[u]<deep[v]?u:v;
}
struct SegTree{
ll val[N<<2];
int tag[N<<2];
#define mid ((l+r)>>1)
void pushdown(int root,int l,int r){
if (tag[root]){
tag[root<<1]+=tag[root];
tag[(root<<1)|1]+=tag[root];
val[root<<1]+=1ll*tag[root]*(mid-l+1);
val[(root<<1)|1]+=1ll*tag[root]*(r-mid);
tag[root]=0;
}
}
ll query(int root,int l,int r,int L,int R){
if (r<L||l>R) return 0;
if (L<=l&&r<=R) return val[root];
pushdown(root,l,r);
return query(root<<1,l,mid,L,R)+query((root<<1)|1,mid+1,r,L,R);
}
void update(int root,int l,int r,int L,int R,int v){
if (r<L||l>R) return;
if (L<=l&&r<=R){
val[root]+=1ll*(r-l+1)*v;
tag[root]+=v;
return;
}
pushdown(root,l,r);
update(root<<1,l,mid,L,R,v);
update((root<<1)|1,mid+1,r,L,R,v);
val[root]=val[root<<1]+val[(root<<1)|1];
}
#undef mid
}T;
int rt[N],cnt=0;
struct Seg2{
#define mid ((l+r)>>1)
int ls[N*30],rs[N*30];
void insert(int&root,int l,int r,int x){
if (!root) root=++cnt;
if (l==r) return;
if (x<=mid) insert(ls[root],l,mid,x);
else insert(rs[root],mid+1,r,x);
}
bool used(int root,int l,int r,int L,int R){
if ((!root)||r<L||l>R) return false;
if (L<=l&&r<=R) return true;
return used(ls[root],l,mid,L,R)||used(rs[root],mid+1,r,L,R);
}
#undef mid
}Tree;
int n,q;
bool vis(int c,int x){
bool f=false;
while (x){
f|=Tree.used(rt[c],1,n,dfn[top[x]],dfn[x]);
x=fa[top[x]];
}
return f;
}
void ins(int c,int x){
Tree.insert(rt[c],1,n,dfn[x]);
}
int main (){
freopen ("snowcow.in","r",stdin);
freopen ("snowcow.out","w",stdout);
scanf ("%d%d",&n,&q);
for (int i=1;i<n;i++){
int u,v;scanf ("%d%d",&u,&v);
addedge(u,v);
}
deep[1]=1;dfs(1,0),dfs2(1,1);
while (q--){
int opt;scanf ("%d",&opt);
if (opt==1){
int x,c;scanf ("%d%d",&x,&c);
// fprintf (stderr,"%d %d
",x,c);
//Case1
if (vis(c,x)) continue;ins(c,x);
//Case2
int l=dfn[x],r=dfn[x]+size[x]-1;
set < pair<int,int> > :: iterator it;
while (1){
it=S[c].lower_bound(make_pair(l,r));
if (it==S[c].end()) break;
if ((*it).first>r) break;
int u=to[(*it).first];
T.update(1,1,n,dfn[u],dfn[u]+size[u]-1,-1);
S[c].erase(it);
}
S[c].insert(make_pair(l,r));
T.update(1,1,n,dfn[x],dfn[x]+size[x]-1,1);
}else{
int x;scanf ("%d",&x);
printf ("%lld
",T.query(1,1,n,dfn[x],dfn[x]+size[x]-1));
}
}
return 0;
}