IOI2020集训队作业-6 (CF571E, CF696F, ARC093E)

最新推荐文章于 2021-02-23 18:54:43 发布
最新推荐文章于 2021-02-23 18:54:43 发布
阅读量553 收藏 1
点赞数 1
文章标签： 集训队作业
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45313881/article/details/102905117
版权

A - CF571E Geometric Progressions

题意

给出 n n n个数列，第 i i i个数列是： a i , a i b i , a i b i 2 , a i b i 3 ⋯ a_i,a_ib_i,a_ib_i^2 ,a_ib_i^3 \cdots ai​,ai​bi​,ai​bi2​,ai​bi3​⋯

问最小的被每一个数列都包含了的数对 1 0 9 + 7 10^9+7 109+7取模的值。无解输出 − 1 -1 −1。

1 ≤ n ≤ 100 , 1 ≤ a i , b i ≤ 1 0 9 1\le n\le 100,1\le a_i,b_i\le 10^9 1≤n≤100,1≤ai​,bi​≤109

Sol

两个等比数列的交，要么也是一个等比数列（公比可能为 1 1 1），要么为空集。

设第一个等比数列的第 x 1 x_1 x1​项和第二个等比数列的第 x 2 x_2 x2​项相同（ x 1 , x 2 ≥ 0 x_1,x_2\ge 0 x1​,x2​≥0），则有

a 1 b 1 x 1 = a 2 b 2 x 2 a_1 b_1^{x_1} = a_2b_2^{x_2} a1​b1x1​​=a2​b2x2​​

对每一个质因子进行考虑，设质因子 p p p在 m m m中的次数为 e p ( m ) e_p(m) ep​(m)，则：

e p ( a 1 ) + e p ( b 1 ) x 1 = e p ( a 2 ) + e p ( b 2 ) x 2 e_{p}(a_1) + e_p(b_1) x_1 = e_p(a_2)+e_p(b_2)x_2 ep​(a1​)+ep​(b1​)x1​=ep​(a2​)+ep​(b2​)x2​

这样我们将得到若干个关于 x 1 x_1 x1​和 x 2 x_2 x2​的二元一次方程。

这个方程组要么无解，要么有无数组解，要么有唯一解。对应的就是等比数列的交集为空，交集为公比不为 1 1 1的等比数列，交集为公比为 1 1 1的等比数列。

每次合并两个等比数列，直到最终只剩下一个等比数列，输出剩下的等比数列的首项就可以了。

Code

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#define ll long long
using namespace std;
template <class T>
inline void rd(T &x) {
	x=0; char c=getchar(); int f=1;
	while(!isdigit(c)) { if(c=='-') f=-1; c=getchar(); }
	while(isdigit(c)) x=x*10-'0'+c,c=getchar(); x*=f;
}
inline ll Abs(ll x) { return x>0?x:-x; }
ll gcd(ll a,ll b) { return b?gcd(b,a%b):a; }
struct item {
	int p; ll k;
	item(int p=0,ll k=0): p(p),k(k) {}
};
struct Num {
	item p[810];
	int num;
	void init(int x) {
		for(int i=2;i*(ll)i<=x;++i) if(x%i==0) {
			int c=0;
			while(x%i==0) x/=i,c++;
			p[num++]=item(i,c);
		}
		if(x>1) p[num++]=item(x,1);
	}
	item query(int x) {
		for(int i=0;i<num;++i) if(p[i].p==x) return p[i];
		return item(x,0);
	}
}A[110],B[110];
int n;
void FAIL() { puts("-1"); exit(0); }

item A1[3210],B1[3210],A2[3210],B2[3210];
int pri[3210],num;
ll exgcd(ll a,ll b,ll &x,ll &y) {
	if(!b) { x=1,y=0; return a; }
	ll g=exgcd(b,a%b,y,x); y-=a/b*x; return g;
}
void Inter(ll A,ll B,ll C,ll a,ll b,ll c,ll &x1,ll &x2) {
	while(a) {
		ll t=A/a;
		A-=t*a,B-=t*b,C-=t*c;
		swap(A,a),swap(B,b),swap(C,c);
	}
	if(c%b) FAIL();
	x2=-c/b;
	if((C+B*x2)%A) FAIL();
	x1=(-C-B*x2)/A;
}
void UN(Num &a1,Num &b1,Num &a2,Num &b2) {
	num=0;
	for(int i=0;i<a1.num;++i) pri[num++]=a1.p[i].p;
	for(int i=0;i<a2.num;++i) pri[num++]=a2.p[i].p;
	for(int i=0;i<b1.num;++i) pri[num++]=b1.p[i].p;
	for(int i=0;i<b2.num;++i) pri[num++]=b2.p[i].p;
	sort(pri,pri+num);
	num=unique(pri,pri+num)-pri;
	for(int i=0;i<num;++i) A1[i]=a1.query(pri[i]);
	for(int i=0;i<num;++i) A2[i]=a2.query(pri[i]);
	for(int i=0;i<num;++i) B1[i]=b1.query(pri[i]);
	for(int i=0;i<num;++i) B2[i]=b2.query(pri[i]);
	
	ll A=0,B=0,C=0;
	int flg1=0; ll x1,x2;
	for(int i=0;i<num;++i) {
		ll a=B1[i].k,b=-B2[i].k,c=A1[i].k-A2[i].k;
		if(a==0&&b==0) {
			if(c) FAIL();
			continue;
		}
		ll g=gcd(a,gcd(-b,Abs(c)));
		a/=g,b/=g,c/=g;
		if(!b) {
			if(c%a) FAIL();
			if(-c/a<0) FAIL();
		}
		if(!A&&!B) { A=a,B=b,C=c; continue; }
		if(!B) {
			if(b) { Inter(A,B,C,a,b,c,x1,x2),flg1=1; break; }
			if(C/A!=c/a) FAIL();
			continue;
		}
		if(A*b==a*B) {
			if(c*A==C*a) continue;
			FAIL();
		}
		Inter(A,B,C,a,b,c,x1,x2),flg1=1;
		break;
	}
	if(flg1) {
		for(int i=0;i<num;++i) {
			ll a=B1[i].k,b=-B2[i].k,c=A1[i].k-A2[i].k;
			if(a*x1+b*x2+c) FAIL();
		}
		for(int i=0;i<num;++i) {
			A1[i].k=A1[i].k+B1[i].k*x1;
			B1[i].k=0;
		}
		for(int i=0;i<num;++i)
			a1.p[i]=A1[i],b1.p[i]=B1[i];
		a1.num=b1.num=num;
		return;
	}
	ll g=exgcd(A,B,x1,x2);
	
	if(C%g) FAIL();
	x1*=-C/g,x2*=-C/g;
	
	ll tx=Abs(-B/g),ty=Abs(A/g);
	if(C>0||!ty) {
		x1=(x1%tx+tx)%tx;
		if(B) x2=-(A*x1+C)/B;
		else x2=0;
	}
	else {
		x2=(x2%ty+ty)%ty;
		if(A) x1=(-B*x2-C)/A;
		else x1=0;
	}
	for(int i=0;i<num;++i) {
		A1[i].k=A1[i].k+B1[i].k*x1;
		B1[i].k=tx*B1[i].k;
	}
	for(int i=0;i<num;++i)
		a1.p[i]=A1[i],b1.p[i]=B1[i];
	a1.num=b1.num=num;
}
const int mod=1e9+7;
int Pow(int x,ll y) {
	int res=1;
	for(y%=(mod-1);y;x=x*(ll)x%mod,y>>=1)
		if(y&1) res=res*(ll)x%mod;
	return res;
}
int main() {
	rd(n);
	for(int i=1,a,b;i<=n;++i) rd(a),rd(b),A[i].init(a),B[i].init(b);
	for(int i=2;i<=n;++i) UN(A[1],B[1],A[i],B[i]);
	long long ans=1;
	for(int i=0;i<A[1].num;++i)
		ans=ans*(ll)Pow(A[1].p[i].p,A[1].p[i].k)%mod;
	printf("%lld",ans);
	return 0;
}

B - CF696F …Dary!

题意

有一个包含 n n n个顶点的凸多边形。

你需要在凸多边形内选两个点，以这两个点为圆心、以某个实数 r r r为半径画两个圆，使得凸多边形的每一条边所在直线至少与一个圆有交。

你需要最小化 r r r。

n ≤ 300 n\le 300 n≤300，坐标绝对值不超过 1 0 4 10^4 104，答案要求精确到 1 0 − 6 10^{-6} 10−6。

Sol

首先二分 r r r。

如果圆与某条边所在直线有交，那么这个圆的圆心一定在这条边向多边形的中心平移 r r r之后得到的直线所划分出的、异于平移方向的那个半平面内。

所以，可以在多边形找出两个满足条件的圆心，等价于能够把这样的 n n n个半平面分成至多两组，每组和多边形的交不为空集。

一种暴力的做法是：圆心一定可以是某两个半平面的交点，我们不妨去枚举这个交点 p p p，然后删掉所有与它有交的半平面，检查剩下的半平面的交与多边形的交是否为空集。这样的复杂度是 O ( n 3 ) O(n^3) O(n3)的。

可以证明，我们的决策一定是将原多边形上连续的一段边对应的半平面分在一组。

考虑 l l l这条边和 r r r这条边所对应的半平面的交点 p p p，如果在 l l l到 r r r之间的存在一条边对应的半平面与 p p p没有交，并且在 r r r到 l l l之间也存在一条边对应的半平面与 p p p没有交，那么 p p p一定不能作为圆心，这是因为那两个半平面与多边形的交一定是空集。

所以可以考虑对于每条边求出一个最大的 l l l，满足它和它之后的 l l l条边对应的半平面的与多边形的交不为空集。 l l l显然具有单调性，可以双指针。

时间复杂度 O ( n 2 log ⁡ V ) O(n^2\log V) O(n2logV)。

Code

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#define ll long long
#define db long double
using namespace std;
template <class T>
inline void rd(T &x) {
	x=0; char c=getchar(); int f=1;
	while(!isdigit(c)) { if(c=='-') f=-1; c=getchar(); }
	while(isdigit(c)) x=x*10-'0'+c,c=getchar(); x*=f;
}
const db eps=1e-9;
const int N=310;
int sgn(db x) { return x>-eps?(x>eps):-1; }
struct Point {
	db x,y;
	Point(db x=0,db y=0): x(x),y(y) {}
	friend Point operator +(Point A,Point B) { return Point(A.x+B.x,A.y+B.y); }
	friend Point operator -(Point A,Point B) { return Point(A.x-B.x,A.y-B.y); }
	friend Point operator *(Point A,db B) { return Point(A.x*B,A.y*B); }
	friend Point operator /(Point A,db B) { return Point(A.x/B,A.y/B); }
};
db Cross(Point A,Point B) { return A.x*B.y-A.y*B.x; }
db len(Point A) { return sqrt(A.x*A.x+A.y*A.y); }
struct Line {
	Point u,v;
	Line(Point u=Point(0,0),Point v=Point(0,0)): u(u),v(v) {}
};
Point Inter(Line A,Line B) {
    db s1=Cross(B.u-A.u,A.v-A.u);
    db s2=Cross(A.v-A.u,B.v-A.u);
    return (B.u*s2+B.v*s1)/(s1+s2);
}
int n;
Line p[N<<1];
int rnk[N<<1];
int pos[N<<1];
int vis[N<<1];
db ag[N<<1];
bool cmpag(int x,int y) { return ag[x]<ag[y]; }
bool HPI(int l,int r,int flg) {
	static Line q[N<<1];
	static Point s[N<<1];
	static int id[N<<1];
	for(int i=l;i!=r;++i) vis[i%n+n]=1;
	int hd=1,tl=0;
	for(int i=0;i<n*2;++i) if(vis[pos[i]]) {
		Line c=p[pos[i]];
		while(hd<tl&&sgn(Cross(c.v-c.u,s[tl]-c.u))<0) tl--;
		while(hd<tl&&sgn(Cross(c.v-c.u,s[hd+1]-c.u))<0) hd++;
		if(hd<=tl&&sgn(Cross(c.v-c.u,q[tl].v-q[tl].u))==0) {
			if(sgn(Cross(c.v-c.u,q[tl].v-c.u))<0) {
				if(sgn(Cross(q[tl].v-q[tl].u,c.v-q[tl].u))>=0)
					q[tl]=c,id[tl]=pos[i];
				else {
					for(int i=l;i!=r;++i) vis[i%n+n]=0;
					return 0;
				}
			}
		}		
		else q[++tl]=c,id[tl]=pos[i];
		if(hd<tl) s[tl]=Inter(q[tl],q[tl-1]);
	}
	while(hd<tl&&sgn(Cross(q[hd].v-q[hd].u,s[tl]-q[hd].u))<0) tl--;
	for(int i=l;i!=r;++i) vis[i%n+n]=0;
	if(tl-hd+1<3) return 0;
	s[hd]=Inter(q[hd],q[tl]);
	if(flg) printf("%.10Lf %.10Lf\n",s[hd].x,s[hd].y);
	return 1;
}
Point P[N];
bool check(db d) {
	for(int i=0;i<n;++i) {
		Point t=p[i].v-p[i].u;
		Point r=Point(-t.y,t.x);
		r=r*(d/len(p[i].v-p[i].u));
		p[i+n].u=p[i].v+r;
		p[i+n].v=p[i].u+r;
	}
	int r=1;
	for(int i=0;i<n;++i) {
		r=max(r,i+1);
		while(r+1<i+n&&HPI(i,r+1,0)) r++;
		if(HPI(r,i+n,0)) return 1;
	}
	return 0;
}
void sol(db d) {
	printf("%.10Lf\n",d);
	for(int i=0;i<n;++i) {
		Point t=p[i].v-p[i].u;
		Point r=Point(-t.y,t.x);
		r=r*(d/len(p[i].v-p[i].u));
		p[i+n].u=p[i].v+r;
		p[i+n].v=p[i].u+r;
	}
	int r=1;
	for(int i=0;i<n;++i) {
		r=max(r,i+1);
		while(r+1<i+n&&HPI(i,r+1,0)) r++;
		if(HPI(r,i+n,0)) {
			HPI(i,r,1);
			HPI(r,i+n,1);
			return;
		}
	}
}
int main() {
	rd(n);
	for(int i=0;i<n;++i) scanf("%Lf%Lf",&P[i].x,&P[i].y);
	p[0]=Line(P[n-1],P[0]);
	for(int i=1;i<n;++i) p[i]=Line(P[i-1],P[i]);
	for(int i=0;i<n;++i) p[i+n]=p[i],swap(p[i+n].u,p[i+n].v);
	for(int i=0;i<2*n;++i) ag[i]=atan2(p[i].v.y-p[i].u.y,p[i].v.x-p[i].u.x),pos[i]=i;
	sort(pos,pos+2*n,cmpag);
	for(int i=0;i<2*n;++i) rnk[pos[i]]=i;
	for(int i=0;i<n;++i) vis[i]=1;
	db l=0,r=1e5;
	while(r-l>1e-7) {
		db mid=(l+r)/2;
		if(check(mid)) r=mid;
		else l=mid;
	}
	sol(r);
	return 0;
}

C - ARC093E Bichrome Spanning Tree

题意

有 n n n个点、 m m m条边的带权无向图。问有多少种把边黑白染色的方式，使得这张图的同时包含了两种颜色的边的生成树，边权和最小为 X X X。答案对 1 0 9 + 7 10^9+7 109+7取模。

n ≤ 1000 , m ≤ 2000 , 1 ≤ w ≤ 1 0 9 , X ≤ 1 0 12 n\le 1000,m\le 2000,1\le w \le 10^9,X\le 10^{12} n≤1000,m≤2000,1≤w≤109,X≤1012（ w w w是边权）

Sol

首先求出最小生成树的边权 W W W，设 Δ = X − W \Delta = X - W Δ=X−W。

如果 Δ < 0 \Delta < 0 Δ<0则无解。

考虑如果最小生成树中所有边都是同一种颜色，我们肯定是从图中另选一条与树边颜色不同的边，然后断掉这条边连接的点在树上的路径之间边权最大的边。

对于每条不在最小生成树上的边求出（这条边的边权 - 这条边连接的点在树上的路径之间的边的最大边权），设 x 0 , x 1 , x 2 x_0,x_1,x_2 x0​,x1​,x2​分别表示这个值小于 Δ \Delta Δ，等于 Δ \Delta Δ，大于 Δ \Delta Δ的边的数量。

如果 Δ = 0 \Delta = 0 Δ=0，那么要么原来的最小生成树就合法，要么最小生成树不合法、但 x 1 x_1 x1​中有边与最小生成树上的边的颜色不一样。答案就是 ( 2 n − 1 − 2 ) ⋅ 2 m − ( n − 1 ) + 2 ⋅ ( 2 x 1 − 1 ) ⋅ 2 x 2 (2^{n-1}-2)\cdot 2^{m-(n-1)} + 2\cdot (2^{x_1}-1)\cdot 2^{x_2} (2n−1−2)⋅2m−(n−1)+2⋅(2x1​−1)⋅2x2​。

如果 Δ > 0 \Delta > 0 Δ>0，那么就要求原来的最小生成树不合法，并且所有 x 0 x_0 x0​的边都与最小生成树上的边的颜色相同，并且存在 x 1 x_1 x1​的边与最小生成树上的边颜色不同。答案就是 2 ⋅ ( 2 x 1 − 1 ) 2 x 2 2\cdot (2^{x_1}-1)2^{x_2} 2⋅(2x1​−1)2x2​。

Code

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#define ll long long
using namespace std;
template <class T>
inline void rd(T &x) {
	x=0; char c=getchar(); int f=1;
	while(!isdigit(c)) { if(c=='-') f=-1; c=getchar(); }
	while(isdigit(c)) x=x*10-'0'+c,c=getchar(); x*=f;
}
const int mod=1e9+7;
int Pow(int x,int y) {
	int res=1;
	while(y) {
		if(y&1) res=res*(ll)x%mod;
		x=x*(ll)x%mod,y>>=1;
	}
	return res;
}
const int N=1010,M=2010;
int head[N],ecnt;
struct ed { int to,next,w; }e[N<<1];
void ad(int x,int y,int w) {
	e[++ecnt]=(ed){y,head[x],w}; head[x]=ecnt;
	e[++ecnt]=(ed){x,head[y],w}; head[y]=ecnt;
}
int gl,qmx;
void dfs(int u,int last,int mx) {
	if(u==gl) { qmx=mx; return; }
	for(int k=head[u];k;k=e[k].next) {
		int v=e[k].to; if(v==last) continue;
		dfs(v,u,max(mx,e[k].w));
	}
}
int Query(int x,int y) {
	gl=y;
	dfs(x,0,0);
	return qmx;
}
struct edge {
	int x,y,w;
	friend bool operator < (edge A,edge B) { return A.w<B.w; }
}E[M];
int used[M];
int fa[N];
int find(int x) { return fa[x]==x?x:fa[x]=find(fa[x]); }
bool un(int x,int y) {
	if(find(x)==find(y)) return 0;
	fa[find(x)]=find(y); return 1;
}
ll X;
int n,m;
int main() {
	rd(n),rd(m),rd(X);
	for(int i=1;i<=m;++i) rd(E[i].x),rd(E[i].y),rd(E[i].w);
	sort(E+1,E+m+1);
	for(int i=1;i<=n;++i) fa[i]=i;
	ll ans=0;
	for(int i=1;i<=m;++i)
		if(un(E[i].x,E[i].y)) ad(E[i].x,E[i].y,E[i].w),ans+=E[i].w,used[i]=1;
	int cnt0=0,cnt1=0,cnt2=0;
	ll del=X-ans;
	if(del<0) {
		printf("0"); return 0;
	}
	for(int i=1;i<=m;++i) if(!used[i]) {
		int t=E[i].w-Query(E[i].x,E[i].y);
		if(t>del) cnt2++;
		else if(t==del) cnt1++;
		else cnt0++;
	}
	if(del==0)
		printf("%d\n",((Pow(2,n-1)-2)*(ll)Pow(2,m-(n-1))%mod+2ll*(Pow(2,cnt1)-1)*Pow(2,cnt2)%mod+mod)%mod);
	else
		printf("%d\n",(2ll*Pow(2,cnt2)*(ll)(Pow(2,cnt1)-1)%mod+mod)%mod);
	return 0;
}
E[i].w,used[i]=1;
	int cnt0=0,cnt1=0,cnt2=0;
	ll del=X-ans;
	if(del<0) {
		printf("0"); return 0;
	}
	for(int i=1;i<=m;++i) if(!used[i]) {
		int t=E[i].w-Query(E[i].x,E[i].y);
		if(t>del) cnt2++;
		else if(t==del) cnt1++;
		else cnt0++;
	}
	if(del==0)
		printf("%d\n",((Pow(2,n-1)-2)*(ll)Pow(2,m-(n-1))%mod+2ll*(Pow(2,cnt1)-1)*Pow(2,cnt2)%mod+mod)%mod);
	else
		printf("%d\n",(2ll*Pow(2,cnt2)*(ll)(Pow(2,cnt1)-1)%mod+mod)%mod);
	return 0;
}
zhongyuwei20031224
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
IOI2021集训队作业
mayaohua
12-21 971
（即将）退役选手诈尸啦！ 今年的集训队作业来自于 21 场 ICPC 竞赛，所有竞赛均在 Codeforces 的 GYM 中可以找到（下面列表中结尾的数字代表其在 GYM 中的场次编号，例如 101221 对应 http://codeforces.com/gym/101221）。 A 2014 ACM-ICPC World Finals，101221。 B 2015 ACM-ICPC World Finals，101239。 C 2016 ACM-ICPC World Finals，101242。 D 2
Atcoder arc093E
mayaohua
11-20 450
考虑先不管染色，任意求出原图的一棵MST，设为TTT。这个可以用经典的kruskal算法解决。 我们称生成树TTT的权重和为w(T)w(T)w(T)。 若w(T)>Xw(T)>Xw(T)>X，显然一定无解。接下来我们假设w(T)≤Xw(T)\leq Xw(T)≤X。 那么对于一个染色方案，若TTT中同时有两种颜色的边，那它本身就是符合题意的生成树中权重和最小的之一。显然当且仅当w...
ioi2020集训队作业_IOI2020集训队作业
weixin_36150766的博客
02-23 327
跟着做一下，不是自己独立思考出来的会打个$\Xi$。(然后全是$\Xi$AGC023E不会，再说。AGC023F感觉怎么写个假贪心就过了啊。感性理解一下，肯定$0$越靠前越好，但可能出现一个子树有$0$又有$1$，那怎么办呢？考虑从单个节点开始，合并出一个又一个连通块，最终合并成一棵树。记录这个连通块里的$0$的个数与$1$的个数，那么我们可以在每次合并的时候算出这次合并所产生的逆序对数。那合并的...
ioi2018集训队作业1
12-26
希望喜欢算法的童鞋一起学习大佬们的思维。IOI是国际信息学奥林匹克竞赛，此文档是中国IOI2018集训队的作业题。
集训队作业
wxh010910
11-04 5104
非常抱歉，这篇文章鸽了。
IOI国家集训队论文集
11-17
包括IOI中国国家集训队从1999年到2017年的论文 部分目录： 国家集训队1999论文集 陈宏：《数据结构的选择与算法效率——从IOI98试题PICTURE谈起》 来煜坤：《把握本质，灵活运用——动态规划的深入探讨》 齐鑫：...
IOI国家集训队论文1999-2021(缺2020)-2021.08.03.rar
08-03
IOI国家集训队论文1999-2021(缺2020)-2021.08.03.rar
IOI国家集训队论文集1999-2019
09-20
# 国家集训队论文列表（1999-2019） ___点击目录快速跳转：___ - _国家集训队论文列表（1999-2019）_ * [_1999_](#1999) * [_2000_](#2000) * [_2001_](#2001) * [_2002_](#2002) * [_2003_](#2003) * [_...
IOI2020中国国家候选队论文集 非正式版1
08-03
【标题】：“IOI2020中国国家候选队论文集 非正式版1” 【描述】：本文主要探讨了分散层叠算法在优化经典分块问题中的应用，这是一种加速多序列二分查找的技术。 【部分内容】：分散层叠算法，又称为Fractional ...
arc 093 E Bichrome Spanning Tree
Rising_shit的博客
06-19 647
         题意：你有一个n个点m条边的有向图，你需要给每条边染黑色或者白色，要求，原图必须存在一个由黑边和白边组成的生成树，且在那些由黑边与白边构成的生成树中，最小的生成树的边权为x，问有多少种染色方案。         首先我们可以求一棵原图的最小生成树，令它的权值为 t，显然如果x大于t是无解的，接着我们来讨论t=x 与t&lt;x 的情况。         当t=x的时候，只要某条边...
【图论】【生成树】(AtCoder Regular Contest 093 E) Bichrome Spanning Tree
氧化钠的博客
03-28 568
题意： 给出一个N个点M条边无向连通图，现在需要给每条边染色（黑/白），染色完毕后，算出一颗至少包含一条黑边与一条白边的最小生成树，要求这个生成树的权值和为X，求染色的方案数。 N≤1000,M≤2000N≤1000,M≤2000N\leq 1000,M\leq 2000 分析： Atcoder的题的确很锻炼思维。 首先，一个众所周知的最小生成树算法：破圈算法。这就是本题解法的基础。...
IOI2020集训队作业-19 (CF571D, CF708E, ARC093F)
这个博客已经搬到了：zhongyuwei.github.io
12-26 457
A - CF571D Campus Sol 首先对大学和警卫室分别建个Kruskal重构树并求出dfs序。询问等价于：询问某个宿舍在最后一次被清空了之后，所有对它的加操作的和。可以差分一下转化成最后一次清空和询问这两个时间点之前的加操作的和，然后对时间做扫描线，用线段树维护加操作即可。 Code #include <cstdio> #include <iostream> #...
CF--Geometric Progression--思维
lanshan1111的博客
08-14 192
C. Geometric Progression time limit per test 1 second memory limit per test 256 megabytes input standard input output standard output Polycarp loves geometric progressions very much. Since he...
CF 567C. Geometric Progression
便纵有千种风情
09-24 528
题目链接题意：给出一个序列,求公比为k，且长度为3的子等比序列个数。(1 ≤ n, k ≤ 2e5)( - 1e9 ≤ ai ≤ 1e9) #解法: dp[0][x]表示当前x的数量。 dp[1][x]表示当前以值x结尾，长度为2且公比为K的等比数列数量。代码：#include<cstring> #include<iostream> #include<map> using namespace
CF 696E ...Wait for it... 链剖练手题
YxuanwKeith的博客
08-18 810
题目大意给定一颗NN个点的树，现在有MM个物品，每个物品一开始都在节点AiA_i上，权值为ii。现在有两种操作，先读入OrdOrd： Ord=1Ord = 1：再读入三个整数u,v,Limu,v,Lim表示从uu到vv这条路径上从小到大把最多LimLim个物品从树上取出，并且从小到大输出。 Ord=2Ord = 2：再读入两个数u,Valu, Val表示给在uu这棵子树的所有物品加上权值ValV
IOI2020第一轮选拔模拟4（线性基、点分治）
chiyankuan的博客
12-09 170
今天的题也很难，目前只改出了两题。 T1：这题要用到一个叫线性基的工具。 首先我们发现3^b给予我们一个贪心的思路：一定是要优先满足b较大的列。 紧接着就有一个做法了：从大到小枚举b，判断计算当前b的最大贡献。如果b有+1而没有-1，那么就判断一下3^b是否可行。有-1而没有+1同理。注意有+1也有-1时两列都为0和两列都为1的情况是同理的，在这种情况下我们可以定义限制为两列异或起来的...
UOJ#419. 【集训队作业2018】圆形（格林公式）
weixin_34245169的博客
04-22 219
题面 传送门 题解 首先您得会用格林公式计算圆的面积并 这里只需要动态维护一下圆弧就可以了 时间复杂度\(O(n^2\log n)\) //minamoto #include<bits/stdc++.h> #define R register #define inline __inline__ __attribute__((always_inline)) #define fp(i,a,...
arc 093 F Dark Horse
Rising_shit的博客
06-19 538
         题意： 我们现在有一个有2的n次方的人的淘汰赛，每次相邻的两个人打，赢得人晋级，以此类推到只剩一个人为止。给定一个集合a[i], 规则如下：1号和a[i]的人打a[i]胜，1号和其他人打 1号胜利，其他号和其他的打，字典序小的胜。问1号胜利的方案数有多少种。        我们先假装第一个人就在第一个位置，由于都是从一个2^1的人中取胜，到从2^n的人中取胜，所以一号人在哪方案数...
【集训队作业】IOI 2020 集训队作业 试题泛做 13
cz_xuyixuan的博客
01-18 873
Codeforces 679E Bear and Bad Powers of 42 不难发现数列中的元素不可能达到很大，我们只需要考虑 424242 的前若干个幂。 考虑没有赋值操作的做法，则可用线段树维护区间中最接近 424242 的下一个幂的数与这个幂的差值，在区间加时，只需要在线段树上 DFS 找到差值变负的位置更新其与 424242 下一个幂的差值即可。 由于一个数导致进行 DFS 的次数...
[IOI2000] 邮局
最新发布
08-17
- *1* *2* [P4767 [IOI2000]邮局（四边形不等式）](https://blog.csdn.net/gyp0205/article/details/123991528)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值