为啥代码在CF里AC,在luogu,过不去?
查看原帖
为啥代码在CF里AC,在luogu,过不去?
609170
xin_fu楼主2022/8/27 18:52
#include<bits/stdc++.h>
#define ll long long
using namespace std;

const int N=2e5+10;

struct edge{
	ll x,y,w;
	ll c,d;
}g[N];

ll n,m,mn,S,sum;
int head[N],to[N<<1],nxt[N<<1],e[N<<1],id[N<<1],tot;
int t[N];int f[N][21],h[N][21],d[N],id1[N][21];
int fa[N];

void add(int x,int y,int w,int i)
{
	to[++tot]=y;
	id[tot]=i;
	nxt[tot]=head[x];
	head[x]=tot;
	e[tot]=w;
}

void dfs(int x)
{
	for(int i=1;i<20;i++)
	{
		if(h[x][i-1]>h[f[x][i-1]][i-1])
		{
			h[x][i]=h[x][i-1];
			id1[x][i]=id1[x][i-1];
		}
		else
		{
			h[x][i]=h[f[x][i-1]][i-1];
			id1[x][i]=id1[f[x][i-1]][i-1];
		}
		f[x][i]=f[f[x][i-1]][i-1];
	}
	for(int i=head[x];i;i=nxt[i])
	{
		int y=to[i],w=e[i];
		if(y==f[x][0])continue;
		f[y][0]=x;
		h[y][0]=w;
		id1[y][0]=id[i];
		d[y]=d[x]+1;
		sum+=w;
		dfs(y);
	}
}

int lca(int x,int y)
{
	int ans;
	if(d[x]>d[y])swap(x,y);
	for(int i=19;i>=0;i--)
	{
		if(d[f[y][i]]>=d[x])
		{
			if(!ans)ans=id1[y][i];
			else if(h[y][i]>g[ans].w)
			{
				ans=id1[y][i];
			}
			y=f[y][i];
		}
	}
	if(x==y)
	return ans;
	for(int i=19;i>=0;i--)
	{
		if(f[x][i]!=f[y][i])
		{
			if(!ans)ans=id1[x][i];
			if(h[x][i]>g[ans].w)
				ans=id1[x][i];
			if(h[y][i]>g[ans].w)
				ans=id1[y][i];
			x=f[x][i],y=f[y][i];
		}
	}
	if(!ans)ans=id1[x][0];
	if(h[x][0]>g[ans].w)
		ans=id1[x][0];
	if(h[y][0]>g[ans].w)
		ans=id1[y][0];
	return ans;
}

bool cmp(edge a,edge b)
{
	return a.w<b.w;
}

int fd(int x)
{
	if(x!=fa[x])fa[x]=fd(fa[x]);
	return fa[x];
}

signed main()
{
	cin>>n>>m;
	for(int i=1;i<=n;i++)fa[i]=i;
	for(int i=1;i<=m;i++)
		cin>>g[i].w;
	for(int i=1;i<=m;i++)
		cin>>g[i].c;
	for(int i=1;i<=m;i++)
	{
		cin>>g[i].x>>g[i].y;
		g[i].d=i;
	}
	cin>>S;
	sort(g+1,g+1+m,cmp);
	int cnt=0;
	for(int i=1;i<=m;i++)
	{
		int x=fd(g[i].x),y=fd(g[i].y),w=g[i].w;
		if(x!=y)
		{
			add(g[i].x,g[i].y,g[i].w,i);
			add(g[i].y,g[i].x,g[i].w,i);
			fa[y]=x;
			if(g[mn].c>g[i].c)
				mn=i;
			t[i]=1;
			cnt++;
			if(mn==0)
				mn=i;
//			cout<<g[i].x<<" "<<g[i].y<<endl;
		}
		if(cnt>=n-1)
			break;
	}
	dfs(1);
	ll tep=S/g[mn].c;
	int t1,t2,t3=0;
	t1=mn,t2=0;
	for(int i=1;i<=m;i++)
	{
		if(t[i] || g[i].c>g[mn].c)continue;
		int l=lca(g[i].x,g[i].y);
		ll res=(g[i].w-g[l].w)*g[i].c;
		if(res>=S)continue;
		res=S-res;
		if((res/g[i].c)>tep)
		{
			tep=res/g[i].c;
			t1=i,t2=l,t3=g[i].w-g[l].w;
		}
	}
	cout<<sum-tep<<endl;
	for(int i=1;i<=m;i++)
	{
		if(t[i] && i!=t2 && i!=t1)
		{
			cout<<g[i].d<<" "<<g[i].w<<endl;
		}
		else if(i==t1)
		{
			cout<<g[i].d<<" "<<g[i].w-tep-t3<<endl;
		}
	}
	return 0;
}
2022/8/27 18:52
加载中...