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

智能推荐

艾美捷Epigentek DNA样品的超声能量处理方案-程序员宅基地

文章浏览阅读15次。空化气泡的大小和相应的空化能量可以通过调整完全标度的振幅水平来操纵和数字控制。通过强调超声技术中的更高通量处理和防止样品污染,Epigentek EpiSonic超声仪可以轻松集成到现有的实验室工作流程中,并且特别适合与表观遗传学和下一代应用的兼容性。Epigentek的EpiSonic已成为一种有效的剪切设备,用于在染色质免疫沉淀技术中制备染色质样品,以及用于下一代测序平台的DNA文库制备。该装置的经济性及其多重样品的能力使其成为每个实验室拥有的经济高效的工具,而不仅仅是核心设施。

11、合宙Air模块Luat开发:通过http协议获取天气信息_合宙获取天气-程序员宅基地

文章浏览阅读4.2k次,点赞3次,收藏14次。目录点击这里查看所有博文  本系列博客,理论上适用于合宙的Air202、Air268、Air720x、Air720S以及最近发布的Air720U(我还没拿到样机,应该也能支持)。  先不管支不支持,如果你用的是合宙的模块,那都不妨一试,也许会有意外收获。  我使用的是Air720SL模块,如果在其他模块上不能用,那就是底层core固件暂时还没有支持,这里的代码是没有问题的。例程仅供参考!..._合宙获取天气

EasyMesh和802.11s对比-程序员宅基地

文章浏览阅读7.7k次,点赞2次,收藏41次。1 关于meshMesh的意思是网状物,以前读书的时候,在自动化领域有传感器自组网,zigbee、蓝牙等无线方式实现各个网络节点消息通信,通过各种算法,保证整个网络中所有节点信息能经过多跳最终传递到目的地,用于数据采集。十多年过去了,在无线路由器领域又把这个mesh概念翻炒了一下,各大品牌都推出了mesh路由器,大多数是3个为一组,实现在面积较大的住宅里,增强wifi覆盖范围,智能在多热点之间切换,提升上网体验。因为节点基本上在3个以内,所以mesh的算法不必太复杂,组网形式比较简单。各厂家都自定义了组_802.11s

线程的几种状态_线程状态-程序员宅基地

文章浏览阅读5.2k次,点赞8次,收藏21次。线程的几种状态_线程状态

stack的常见用法详解_stack函数用法-程序员宅基地

文章浏览阅读4.2w次,点赞124次,收藏688次。stack翻译为栈,是STL中实现的一个后进先出的容器。要使用 stack,应先添加头文件include<stack>,并在头文件下面加上“ using namespacestd;"1. stack的定义其定义的写法和其他STL容器相同, typename可以任意基本数据类型或容器:stack<typename> name;2. stack容器内元素的访问..._stack函数用法

2018.11.16javascript课上随笔(DOM)-程序员宅基地

文章浏览阅读71次。<li> <a href = "“#”>-</a></li><li>子节点:文本节点(回车),元素节点,文本节点。不同节点树:  节点(各种类型节点)childNodes:返回子节点的所有子节点的集合,包含任何类型、元素节点(元素类型节点):child。node.getAttribute(at...

随便推点

layui.extend的一点知识 第三方模块base 路径_layui extend-程序员宅基地

文章浏览阅读3.4k次。//config的设置是全局的layui.config({ base: '/res/js/' //假设这是你存放拓展模块的根目录}).extend({ //设定模块别名 mymod: 'mymod' //如果 mymod.js 是在根目录,也可以不用设定别名 ,mod1: 'admin/mod1' //相对于上述 base 目录的子目录}); //你也可以忽略 base 设定的根目录,直接在 extend 指定路径(主要:该功能为 layui 2.2.0 新增)layui.exten_layui extend

5G云计算:5G网络的分层思想_5g分层结构-程序员宅基地

文章浏览阅读3.2k次,点赞6次,收藏13次。分层思想分层思想分层思想-1分层思想-2分层思想-2OSI七层参考模型物理层和数据链路层物理层数据链路层网络层传输层会话层表示层应用层OSI七层模型的分层结构TCP/IP协议族的组成数据封装过程数据解封装过程PDU设备与层的对应关系各层通信分层思想分层思想-1在现实生活种,我们在喝牛奶时,未必了解他的生产过程,我们所接触的或许只是从超时购买牛奶。分层思想-2平时我们在网络时也未必知道数据的传输过程我们的所考虑的就是可以传就可以,不用管他时怎么传输的分层思想-2将复杂的流程分解为几个功能_5g分层结构

基于二值化图像转GCode的单向扫描实现-程序员宅基地

文章浏览阅读191次。在激光雕刻中,单向扫描(Unidirectional Scanning)是一种雕刻技术,其中激光头只在一个方向上移动,而不是来回移动。这种移动方式主要应用于通过激光逐行扫描图像表面的过程。具体而言,单向扫描的过程通常包括以下步骤:横向移动(X轴): 激光头沿X轴方向移动到图像的一侧。纵向移动(Y轴): 激光头沿Y轴方向开始逐行移动,刻蚀图像表面。这一过程是单向的,即在每一行上激光头只在一个方向上移动。返回横向移动: 一旦一行完成,激光头返回到图像的一侧,准备进行下一行的刻蚀。

算法随笔:强连通分量-程序员宅基地

文章浏览阅读577次。强连通:在有向图G中,如果两个点u和v是互相可达的,即从u出发可以到达v,从v出发也可以到达u,则成u和v是强连通的。强连通分量:如果一个有向图G不是强连通图,那么可以把它分成躲个子图,其中每个子图的内部是强连通的,而且这些子图已经扩展到最大,不能与子图外的任一点强连通,成这样的一个“极大连通”子图是G的一个强连通分量(SCC)。强连通分量的一些性质:(1)一个点必须有出度和入度,才会与其他点强连通。(2)把一个SCC从图中挖掉,不影响其他点的强连通性。_强连通分量

Django(2)|templates模板+静态资源目录static_django templates-程序员宅基地

文章浏览阅读3.9k次,点赞5次,收藏18次。在做web开发,要给用户提供一个页面,页面包括静态页面+数据,两者结合起来就是完整的可视化的页面,django的模板系统支持这种功能,首先需要写一个静态页面,然后通过python的模板语法将数据渲染上去。1.创建一个templates目录2.配置。_django templates

linux下的GPU测试软件,Ubuntu等Linux系统显卡性能测试软件 Unigine 3D-程序员宅基地

文章浏览阅读1.7k次。Ubuntu等Linux系统显卡性能测试软件 Unigine 3DUbuntu Intel显卡驱动安装,请参考:ATI和NVIDIA显卡请在软件和更新中的附加驱动中安装。 这里推荐: 运行后,F9就可评分,已测试显卡有K2000 2GB 900+分,GT330m 1GB 340+ 分,GT620 1GB 340+ 分,四代i5核显340+ 分,还有写博客的小盒子100+ 分。relaybot@re...

推荐文章

热门文章

相关标签