• USACO2019 DEC 部分题解


    因为一些翻译问题铂金只过了 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) 的最小代价,显然有转移

    [f_i=min{ f_j +cost(i,j,c)} (jle i-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) 为根节点的树,支持以下两种操作:

    1. 给定 (x,c) ,给 (x) 的子树中每个点的属性集合里插入一种属性 (c)
    2. 给定 (x) ,求 (x) 的子树内所有点各有多少种不同的属性,需要输出属性个数之和

    (nle 10^5,qle 10^5)

    dfs 序上子树表示为一段连续的区间,因此两个子树要么包含,要么不相交

    考虑修改一个点,影响的范围是 (x) 的子树,我们给这个范围内所有点的权值都 +1

    但是这样会产生冲突,如果一个属性已经有过了,那么就会有重复,为了消除重复下面分两种情况:

    1. 考虑当前操作的区间被某个区间包含过,那么这是个废区间,直接不做
    2. 如果这个区间包含之前的某些区间,那么撤销那些区间的影响,再加上这个区间

    实现方面,可以用 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;
    }
    
  • 相关阅读:
    企业库相关资源下载Enterprise Library Download[2007.5.29更新]
    企业库EntLib初识(3)配置管理应用程序块ConfigurationApplicationBlock
    我喜欢博客园
    今天是个值得纪念的日子
    Three collaboration teamwork Tools
    How to do ‘undo’ in TortoiseSVN
    继承windows service的安装类出现的多余服务被安装的问题
    List Sort 使用的注意几点
    下面两段代码,你认为那一段写的更好?
    asp.net 2.0教程 c#中的接口、抽象、委托
  • 原文地址:https://www.cnblogs.com/crazyzh/p/12070114.html
Copyright © 2020-2023  润新知