20200519 hz【对称中心,前缀和优化DP,LCT最小生成树(边权k±x)】_adia前缀的含义-程序员宅基地

技术标签: 最小生成树  DP优化  LCT  

T1:小B的棋盘

在这里插入图片描述
K个棋子可以不用完,对称指旋转180°后完全重合。
在这里插入图片描述
而前 k + 1 k+1 k+1与后 k + 1 k+1 k+1个点中一定有一对点匹配,这样就可以确定对称中心,枚举后check即可, O ( k 2 n ) O(k^2n) O(k2n)

Code:

#include<bits/stdc++.h>
#define maxn 100005
using namespace std;
int n,K;
struct node{
    
	int x,y;
	bool operator < (const node &p)const{
    return x==p.x?y<p.y:x<p.x;}
	node operator + (const node &p)const{
    return (node){
    x+p.x,y+p.y};}
	bool operator == (const node &p)const{
    return x==p.x&&y==p.y;}
}a[maxn];
map<node,bool>vis;
bool check(node o){
    
	if(vis[o]) return 0;
	vis[o]=1;
	int cnt=0;
	for(int l=1,r=n;l<=r;){
    
		if(a[l]+a[r]==o) l++,r--;
		else{
    
			if(++cnt>K) return 0;
			if((a[l]+a[r])<o) l++;
			else r--;
		}
	}
	return 1;
}
int main()
{
    
	scanf("%d%d",&n,&K);
	if(n<=K) return puts("-1"),0;
	for(int i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y);
	sort(a+1,a+1+n);
	int ans=0;
	for(int i=1;i<=K+1;i++) 
		for(int j=n-K+i-1;j<=n;j++) 
			ans+=check(a[i]+a[j]);
	printf("%d\n",ans);
}
T2:小B的夏令营

在这里插入图片描述
在这里插入图片描述
改为直接对 r r r计算 l l l前缀和是优化的关键。

Code:

#include<bits/stdc++.h>
#define maxn 1505
using namespace std;
const int mod = 1e9+7;
int n,m,p,k,g[maxn],sg[maxn],sum[maxn],psum[maxn],f[maxn],fac[100005],inv[100005];
int C(int n,int m){
    return 1ll*fac[n]*inv[m]%mod*inv[n-m]%mod;}
int Pow(int a,int b){
    int s=1;for(;b;b>>=1,a=1ll*a*a%mod) if(b&1) s=1ll*s*a%mod; return s;}
void Pre(){
    
	for(int i=1;i<=m;i++) sum[i]=(sum[i-1]+f[i])%mod,psum[i]=(psum[i-1]+1ll*g[i]*sum[i])%mod;
}
int main()
{
    
	int a,b;
	scanf("%d%d%d%d%d",&n,&m,&a,&b,&k),p=1ll*a*Pow(b,mod-2)%mod;
	fac[0]=inv[0]=1; for(int i=1;i<=k;i++) fac[i]=1ll*fac[i-1]*i%mod;
	inv[k]=Pow(fac[k],mod-2); for(int i=k-1;i>=1;i--) inv[i]=1ll*inv[i+1]*(i+1)%mod;
	for(int i=0,lim=min(k,m);i<=lim;i++) g[i]=1ll*C(k,i)*Pow(p,i)%mod*Pow(1-p,k-i)%mod;
	sg[0]=g[0]; for(int i=1;i<=m;i++) sg[i]=(sg[i-1]+g[i])%mod;
	sum[m]=1;
	for(int t=1;t<=n;t++,Pre()){
    
		for(int i=1;i<=m;i++) f[i]=1ll*g[m-i]*(1ll*sg[i-1]*(sum[m]-sum[m-i])%mod-psum[i-1])%mod;
	}
	printf("%d\n",(sum[m]+mod)%mod);
}
T3:小B的图

在这里插入图片描述
在这里插入图片描述
也就是说 k k k决定了负边的优先度,而不是发生替换的时间(不相交的不影响,相交的 k k k大的能替换 k k k小的也能替换)

Code:

#include<bits/stdc++.h>
#define maxn 300005
#define LL long long
using namespace std;
char cb[1<<20],*cs,*ct;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<20,stdin),cs==ct)?0:*cs++)
inline void read(int &a){
    
	char c;bool f=0;while(!isdigit(c=getc())) c=='-'&&(f=1);
	for(a=c-'0';isdigit(c=getc());a=a*10+c-'0'); f&&(a=-a);
}
const int inf = 1e9+5;
int n,m,A,B,cnt;
pair<int,LL>ans[maxn];
struct node{
    
	int x,y,w;
	bool operator < (const node &p)const{
    return w<p.w;}
}a[maxn],b[maxn];
namespace LCT{
    
	#define pa fa[x]
	int fa[maxn],ch[maxn][2],mxp[maxn],val[maxn];
	bool rev[maxn];
	bool isr(int x){
    return ch[pa][0]!=x&&ch[pa][1]!=x;}
	bool isc(int x){
    return ch[pa][1]==x;}
	void upd(int x){
    
		mxp[x]=x;
		if(val[mxp[x]]<val[mxp[ch[x][0]]]) mxp[x]=mxp[ch[x][0]];
		if(val[mxp[x]]<val[mxp[ch[x][1]]]) mxp[x]=mxp[ch[x][1]];
	}
	void rot(int x){
    
		int y=fa[x],z=fa[y],c=isc(x);
		if(!isr(y)) ch[z][isc(y)]=x;
		fa[ch[y][c]=ch[x][!c]]=y;
		fa[ch[x][!c]=y]=x,fa[x]=z;
		upd(y);
	}
	void pd(int x){
    if(rev[x]) swap(ch[x][0],ch[x][1]),rev[ch[x][0]]^=1,rev[ch[x][1]]^=1,rev[x]=0;}
	void pdpath(int x){
    if(!isr(x)) pdpath(pa); pd(x);}
	void splay(int x){
    
		for(pdpath(x);!isr(x);rot(x)) if(!isr(pa)) rot(isc(pa)==isc(x)?pa:x); upd(x);
	}
	void access(int x){
    
		for(int y=0;x;x=fa[y=x]) splay(x),ch[x][1]=y,upd(x);
	}
	void bert(int x){
    access(x),splay(x),rev[x]^=1;}
	int sert(int x){
    access(x),splay(x);for(;ch[x][0];x=ch[x][0]);return x;}
	int split(int x,int y){
    bert(x),access(y),splay(y);return y;}
	void link(int x,int y){
    bert(x),fa[x]=y;}
	void cut(int x,int y){
    split(x,y),fa[x]=ch[y][0]=0,upd(y);}
}
using namespace LCT;
int main()
{
    
	read(n),read(A),read(B),read(m);
	for(int i=1;i<=A;i++) read(a[i].x),read(a[i].y),read(a[i].w);
	for(int i=1;i<=B;i++) read(b[i].x),read(b[i].y),read(b[i].w);
	sort(a+1,a+1+A),sort(b+1,b+1+B);
	for(int i=1;i<=n;i++) val[i]=-inf,mxp[i]=i; val[0]=-inf;
	ans[0].first=-inf;
	for(int i=1;i<=A;i++) if(bert(a[i].x),sert(a[i].y)!=a[i].x){
    
		val[n+i]=a[i].w,mxp[n+i]=n+i;
		link(n+i,a[i].x),link(n+i,a[i].y);
		ans[0].second+=a[i].w;
	}
	for(int i=1;i<=B;i++){
    
		int t=mxp[split(b[i].x,b[i].y)];
		if(val[mxp[t]]==-inf) continue;
		int tmp=b[i].w-a[t-n].w;
		ans[++cnt]=make_pair(tmp/2+(tmp>0&&(tmp&1)),tmp);
		cut(t,a[t-n].x),cut(t,a[t-n].y);
		link(b[i].x,b[i].y);
	}
	sort(ans+1,ans+1+cnt);
	for(int i=1;i<=cnt;i++) ans[i].second+=ans[i-1].second;
	for(int v,l,r,mid;m--;){
    
		read(v),l=0,r=cnt;
		while(l<r) mid=(l+r+1)>>1,ans[mid].first>v?(r=mid-1):(l=mid);
		printf("%lld\n",ans[l].second+1ll*(n-1-2*l)*v);
	}
}
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/C20181220_xiang_m_y/article/details/106220354

智能推荐

计算机毕业设计Java疫情防控医用品管理(系统+源码+mysql数据库+Lw文档)_疫情防护用品销售管理系统 论文-程序员宅基地

文章浏览阅读467次。计算机毕业设计Java疫情防控医用品管理(系统+源码+mysql数据库+Lw文档)springboot基于SpringBoot的婚庆策划系统的设计与实现。JSP健身俱乐部网站设计与实现sqlserver和mysql。JSP网上测试系统的研究与设计sqlserver。ssm基于SpringMvC的流浪狗领养系统。ssm基于Vue.js的音乐播放器设计与实现。ssm校园流浪猫图鉴管理系统的设计与实现。_疫情防护用品销售管理系统 论文

android插件化开发打包,Android项目开发如何设计整体架构-程序员宅基地

文章浏览阅读988次,点赞28次,收藏28次。最后小编想说:不论以后选择什么方向发展,目前重要的是把Android方面的技术学好,毕竟其实对于程序员来说,要学习的知识内容、技术有太多太多,要想不被环境淘汰就只有不断提升自己,从来都是我们去适应环境,而不是环境来适应我们!这里附上我整理的几十套腾讯、字节跳动,京东,小米,头条、阿里、美团等公司19年的Android面试题。把技术点整理成了视频和PDF(实际上比预期多花了不少精力),包含知识脉络 + 诸多细节。由于篇幅有限,这里以图片的形式给大家展示一小部分。

基于单片机数码管秒表控制系统设计-程序员宅基地

文章浏览阅读600次,点赞11次,收藏6次。*单片机设计介绍,基于单片机数码管秒表控制系统设计。

Python小程序之验证码图片生成_小程序图片验证码后端生成-程序员宅基地

文章浏览阅读235次。python小程序之验证码图片的生成定义随机字母的生成函数定义随机颜色生成函数,采用RGB格式,生成一个元组调用Image,生成画布,填充底色为白色调用画笔函数Draw,传入画布对象填充画布的每一个色块,作为背景在画布上控制间距,填上每一个字在最后的图上进行模糊操作代码# 生成一个随机的二维码小程序from PIL import Image,ImageDraw,ImageF..._小程序图片验证码后端生成

思科自防御网络安全方案典型配置_思科设备怎么ranga)服务器区域独立防护;-程序员宅基地

文章浏览阅读2.2k次。 1. 用户需求分析客户规模:客户有一个总部,具有一定规模的园区网络; 一个分支机构,约有20-50名员工; 用户有很多移动办公用户 客户需求:组建安全可靠的总部和分支LAN和WAN; 总部和分支的终端需要提供安全防护,并实现网络准入控制,未来实现对VPN用户的网络准入检查; 需要提供IPSEC/SSLVPN接入; 在内部各主要部门间,及内外网络间进_思科设备怎么ranga)服务器区域独立防护;

苹果账号迁移流程_apple 账号迁移-程序员宅基地

文章浏览阅读445次。4、转移账号生成的 p8 文件(证书文件)1、转移苹果账号的 teamID。2、接受苹果账号的 teamID。5、接受账号生成的 p8 文件。3、转移应用的 AppID。_apple 账号迁移

随便推点

深度学习中优化方法之动量——momentum、Nesterov Momentum、AdaGrad、Adadelta、RMSprop、Adam_momentum seg-程序员宅基地

文章浏览阅读1k次。https://blog.csdn.net/u012328159/article/details/80311892_momentum seg

动态数据生成静态html页_监听数据变更自动生成静态html-程序员宅基地

文章浏览阅读816次。主要的原理就是替换模板里的特殊字符。 1、静态模板页面 template.html,主要是定义了一些特殊字符,用来被替换。 HTML code DOCTYPE HT_监听数据变更自动生成静态html

预防按钮的多次点击 恶意刷新-程序员宅基地

文章浏览阅读494次。 今日在做一个新闻系统的评论时. 想到了预防"提交"按钮的多次点击的问提 (prevent multiple clicks of a submit button in ASP.NET). 以前碰到此类问提总是用重定位页面来解决. 这次我想找到一个一劳永逸的办法. 通过查讯Google,找到了一些代码,挑选一些较好的修改了一下。public void pa

sokcs5软件dante配置指南_dante 代理 配置pam用户名密码 模式-程序员宅基地

文章浏览阅读4.7k次。近来公司业务有需要做socks5代理的需求,研究了一下,主要的开源实现有2个:dante http://www.inet.no/dante/ss5 http://ss5.sourceforge.net/比较了一下,还是比较倾向于dante,因为看到有人这样评价ss5:Project has an incredibly poor source code quality. Th_dante 代理 配置pam用户名密码 模式

Excel vba 求助。_vba countifs 源码-程序员宅基地

文章浏览阅读809次。在excel vba 中用到countifs 函数,但用来统计带有特殊符号* 时总是统计chu_vba countifs 源码

web前端基础——实现动画效果_web前端实现图片动画效果-程序员宅基地

文章浏览阅读2.6k次。当两个效果之间变换时,可以使用transition过渡属性,但是有多个效果来回变换时,就需要使用动画效果,且动画过程可控(重复播放,画面暂停,最终画面等)文章目录1、简介2、实现步骤3、复合属性animation4、动画属性1、简介动画的本质是快速切换大量图片在人脑中形成的具有连续性的画面构成动画的最小单元:帧或者动画帧2、实现步骤定义动画@keyframes 动画名称{ from{} to{}}@keyframes 动画名称{ 0%{} 10%{} 20%{} 50._web前端实现图片动画效果

推荐文章

热门文章

相关标签