HDU - 6982 J - Road Discount wqs二分 + 模型转换 + 优化

传送门

文章目录

  • 题意:
  • 思路:

题意:

给你一个nnn个点mmm条边的图,每个边有一个代价以及折扣价,你需要输出nnn行,第iii行代表你可以选i−1i-1i1条边使其变成优惠价,问每次的最小生成树的代价是多少。
n≤1e3,m≤2e5,ci,di≤1e3n\le 1e3,m\le2e5,c_i,d_i\le 1e3n1e3,m2e5,ci,di1e3

思路:

直接考虑折扣价不是很好想,所以考虑能不能把这些边单独拿出来。
下面我们假定原边是白边,折扣边是黑边,那么对于每次要输出的,问题就转换成了选kkk条黑边的最小生成树的代价是多少。
显然这是一个wqswqswqs二分的一个经典问题,我们设这个函数是f(k)f(k)f(k),这是一个凸函数,我们二分一个值midmidmid,之后将所有黑边的权值都加上midmidmid,让后跑最小生成树,假设选择了cntcntcnt条黑边,且总代价是sumsumsum,那么如果cnt>=kcnt>=kcnt>=k的话,显然可以更新ans=sum−k∗midans=sum-k*midans=sumkmid ,让后调整一下左右边界即可。
直接跑的话复杂度是O(nmlognlogn)O(nmlognlogn)O(nmlognlogn)的,虽然第二个logloglog是最小生成树的,常数很小,但仍是过不了,考虑优化。
考虑每次都有很多无用边,即非树边是无用的,所以直接去掉非树边即可,将边缩小到O(n)O(n)O(n)级别的,但是n2lognlognn^2lognlognn2lognlogn想过这个有101010个测试点的题还是不可能的。
继续优化,考虑对于每次二分,他的信息是可以复用的,且边权在[1,1000][1,1000][1,1000]范围内,所以可以预处理出来,之后每次查询二分的话直接使用已有信息即可。
复杂度O(n2a+nlogn)O(n^2a+nlogn)O(n2a+nlogn),其中n2an^2an2aaaa是并查集的常数,很小。

// Problem: J - Road Discount
// Contest: Virtual Judge - 2021多校第三场补题
// URL: https://vjudge.net/contest/449636#problem/J
// Memory Limit: 524 MB
// Time Limit: 6000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
int p[N],ans,cnt;
PII f[N];
struct Node {int x,y,w,add,op;bool operator < (const Node &W) const {return w<W.w;}
}edge1[N],edge2[N],edge[N];int find(int x) {return x==p[x]? x:p[x]=find(p[x]);
}PII check(int mid) {for(int i=1;i<=m;i++) edge2[i].w+=mid;int tot=0;edge1[m+1].w=INF; edge2[m+1].w=INF;for(int i=1,j=1;i<=m||j<=m;) {if(edge1[i].w<edge2[j].w) edge[++tot]=edge1[i++];else edge[++tot]=edge2[j++];}cnt=0; ans=0;for(int i=1;i<=n;i++) p[i]=i;for(int i=1;i<=tot;i++) {int a=edge[i].x,b=edge[i].y,w=edge[i].w,op=edge[i].op;int pa=find(a),pb=find(b);if(pa==pb) continue;p[pa]=pb; cnt+=op==0;ans+=w;}for(int i=1;i<=m;i++) edge2[i].w-=mid;return {cnt,ans};
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);int _; scanf("%d",&_);while(_--) {scanf("%d%d",&n,&m);for(int i=1;i<=m;i++) scanf("%d%d%d%d",&edge1[i].x,&edge1[i].y,&edge1[i].w,&edge1[i].add),edge1[i].op=1;for(int i=1;i<=m;i++) {edge2[i]=edge1[i],edge2[i].w=edge1[i].add;edge2[i].op=0;}sort(edge1+1,edge1+1+m); sort(edge2+1,edge2+1+m);int tot=0;for(int i=1;i<=n;i++) p[i]=i;for(int i=1;i<=m;i++) {int a=edge1[i].x,b=edge1[i].y,w=edge1[i].w,op=edge1[i].op;int pa=find(a),pb=find(b);if(pa==pb) continue;p[pa]=pb; edge1[++tot]=edge1[i];}tot=0;for(int i=1;i<=n;i++) p[i]=i;for(int i=1;i<=m;i++) {int a=edge2[i].x,b=edge2[i].y,w=edge2[i].w,op=edge2[i].op;int pa=find(a),pb=find(b);if(pa==pb) continue;p[pa]=pb; edge2[++tot]=edge2[i];}m=n-1;for(int i=-1010;i<=1010;i++) f[i+1010]=check(i);for(int k=0;k<n;k++) {int l=-1010,r=1010,res;while(l<=r) {int mid=(l+r)/2;if(f[mid+1010].X>=k) res=f[mid+1010].Y-mid*k,l=mid+1;else r=mid-1;}printf("%d\n",res);}}return 0;
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314753.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

你可能不知道的Docker资源限制

本篇内容涉及Docker的内存与CPU限制&#xff0c;可以用于在实际开发中为指定容器设置限制最大使用的资源量&#xff0c;预计阅读时间为5分钟。01—What is 资源限制&#xff1f;默认情况下&#xff0c;容器是没有资源限制的&#xff0c;它会尽可能地使用宿主机能够分配给它的资…

2021牛客暑期多校训练营4 E - Tree Xor 线段树 + 拆分区间

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一棵树&#xff0c;每个点原本都有一个权值wiw_iwi​&#xff0c;但是你只知道相邻两个点之间的wu⊕wvw_u\oplus w_vwu​⊕wv​&#xff0c;问你有多少种w1,2,...,nw_{1,2,...,n}w1,2,...,n​ n≤1e5,wi…

东南亚的IT公司,我劝你善良!

来源公众号&#xff1a;半佛仙人&#xff08;ID&#xff1a;banfoSB&#xff09;“真的救我一条‘狗命’&#xff0c;可以吗&#xff0c;谢谢您了&#xff0c;我真的撑不住了。”5月3日&#xff0c;小巴在朋友圈看到这条信息&#xff0c;附着一张长图&#xff0c;定位&#xff…

2021牛客暑期多校训练营4 B - Sample Game 期望dp\生成函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个生成器&#xff0c;每次生成1−n1-n1−n其中的某个数的概率为pip_ipi​&#xff0c;生成的规则如下&#xff1a; (1)(1)(1)随机生成一个数加入集合。 (2)(2)(2)判断生成的数是否是集合中的最大值&am…

小白开学Asp.Net Core 《九》

小白开学Asp.Net Core 《九》 — — 前端篇&#xff08;不务正业&#xff09;在《小白开学Asp.Net Core 三》中使用了X-admin 2.x 和 Layui将管理后端的界面重新布局了&#xff0c;里面简单的介绍了layui table 的使用以及页面table所需的数据做了简单的封装。今天扩展学习下。…

P4173 残缺的字符串 FFT匹配含有通配符的字符串

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个长度为m,nm,nm,n的串a,ba,ba,b&#xff0c;问你bbb串中每个长度为mmm的连续字串能否与aaa完全匹配&#xff0c;其中含有通配符∗*∗&#xff0c;输出每个位置的开头。 n,m≤3e5n,m\le3e5n,m≤3e5 思…

.net持续集成sonarqube篇之sonarqube安装与基本配置

Sonarqube下载与安装Sonarqube下载地址是:https://www.sonarqube.org/downloads/下载版本有两个,一个是长期支持版,另一个是最新版,此处安装的是最新版,目前版本是7.3,下载的时候点击醒目的蓝色按钮即可(此时下载的是社区版),下面有三个无底色按钮下载链接,分别对应的是开发者版…

Codeforces Round #296 (Div. 1) D. Fuzzy Search FFT匹配字符串

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n,m,k≤2e5n,m,k\le2e5n,m,k≤2e5 思路&#xff1a; 直接考虑fftfftfft来匹配字符串。 由于kkk是给定的&#xff0c;所以难度低了很多&#xff0c;普通的字符串匹配不能处理这种可以范围相等的情况&#xf…

迁移 WinForm 应用从 dotnet framework 到 dotnetcore3.0

迁移 WinForm 应用从 dotnet framework 到 dotnetcore3.0Intro微软从 dotnetcore3.0 开始已经开始支持 wpf 以及 winform 的开发&#xff0c;dotnet core 3.0 preview7 已经发布&#xff0c;官方称已经可以准备上生产了 Production Ready迁移这篇WPF的迁移还是比较不错的&#…

HDU - 6992 Lawn of the Dead 线段树 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一张n∗mn*mn∗m的图&#xff0c;其中有kkk个点不能走&#xff0c;你只能向下和向右走&#xff0c;问你能到达多少点。 n,m,k≤1e5n,m,k\le1e5n,m,k≤1e5 思路&#xff1a; 可以发现每个点如果其左边和…

EZNEW.NET开发框架100%重磅开源

EZNEW.NET是一套基于领域驱动开发(DDD)为指导思想的企业级项目通用开发框架。通过将当前主流的开发技术和最佳的开发实践相结合&#xff0c;将开发中常见且严重影响开发效率的繁琐技术细节进行了模块化的封装&#xff0c;让开发人员能将更多更多精力聚焦在系统业务的分析中去&a…

黑暗爆炸OJ 3028. 食物 生成函数

传送门 文章目录题意&#xff1a;思考题意&#xff1a; 思考 考虑将每个条件转换成生成函数&#xff1a; (1)f1(x)1x2...11−x2(1)f_1(x)1x^2...\frac{1}{1-x^2}(1)f1​(x)1x2...1−x21​ (2)f2(x)1x1−x21−x(2)f_2(x)1x\frac{1-x^2}{1-x}(2)f2​(x)1x1−x1−x2​ (3)f3(x)1x…

「Azure」数据分析师有理由爱Azure之一-Azure能带给我们什么?

前面我们以相同的方式从数据分析师的视角介绍了Sqlserver&#xff0c;本系列亦同样地延续下去&#xff0c;同样是挖掘数据分析师值得使用的Azure云平台的功能。因云平台功能太多&#xff0c;笔者所接触的面也十分有限&#xff0c;有更专业的读者欢迎补充。对云服务的一点点小认…

POJ - 3734 Blocks 指数生成函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 一段长度为nnn的序列&#xff0c;你有红黄蓝绿四种颜色的砖块&#xff0c;问你铺砖的方案数&#xff0c;每块砖长度为111&#xff0c;其中红黄颜色个数必须为偶数。 思路&#xff1a; 考虑多重集合排列数&…

理想化的DevOps团队里只需要有Dev就够了?

&#xff08;图片来源于网络&#xff09;几天前&#xff0c;本公众号发布的一篇译文列举了9种DevOps团队结构适用类型与7种反型&#xff08;点击查看原文&#xff09;。文章转发到朋友圈之后&#xff0c;很多DevOps同行留言&#xff08;吐槽&#xff09;了自己团队的现状&#…

hdu 1521 排列组合 多重集排列 + 指数生成函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 显然是多重集排列数&#xff0c;我们考虑构造指数生成函数&#xff0c;让后模拟一下多项式乘法即可啦。 由于存在分数&#xff0c;所以直接用doubledoubledouble即可。 //#pragma GCC opt…

架构杂谈《七》

容器VS虚拟机一、什么是虚拟机虚拟机&#xff08;Virtual Machine&#xff09;指通过软件模拟的具有完整硬件系统功能的、运行在一个完全隔离环境中的完整计算机系统。虚拟系统通过生成现有操作系统的全新虚拟镜像&#xff0c;它具有真实windows系统完全一样的功能&#xff0c;…

#3027. [Ceoi2004]Sweet 生成函数 + 组合数学

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 这个题显然可以容斥来写&#xff0c;刚学生成函数就来水一下。 对于每一堆iii我们写出其生成函数Fi(x)∑k0mi(1xx2...xmi)1−x1mi1−xF_i(x)\sum_{k0}^{m_i}(1xx^2...x^{m_i})\frac{1-x^{1…

坚持一个好习惯该有多难?

引子这段时间以来我都在坚持写博客&#xff0c;但是更新得比较少&#xff0c;每周才能一更&#xff0c;不过好歹也勉勉强强也能坚持下去了。我们社区还组建了一个博客互勉群&#xff0c;不过除了个别人能经常打卡外&#xff0c;大部分人都是混个脸熟&#xff0c;凑个热闹。在这…

#3771. Triple 生成函数 + FFT + 容斥

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 注意到这个题是求若干个数的组合数&#xff0c;(a,b),(b,a)(a,b),(b,a)(a,b),(b,a)视为一种方案&#xff0c;所以我们考虑生成一个普通型生成函数。 考虑到每个数只能选一次&#xff0c;但…