【codeforces 733F】Drivers Dissatisfaction
【题目链接】:http://codeforces.com/problemset/problem/733/F
【题意】
给你n个点m条边;
让你从中选出n-1条边;
形成一个生成树;
(即让n个点都联通);
然后,你有S的预算;
每次可以选择一条边i,然后花费ci的预算,把这条边的权值递减1;
(边一开始的权值为wi);
问你最后的最小生成树是多少;
【题解】
/*
肯定是找某一条边一直减(ci最小的那一个,因为代价最小,又都是减少1);
把m条边按照w升序排;
做个最小生成树;
把最小生成树里面c最小的那条边一直减就好,这个作为ans1;
然后在这个最小生成树上写个dfs,
搞出来每个点上面的第2^i个节点是谁,以及这个点到这个点上面的第2^i个点之间,最大的w是哪条边.
然后枚举所有的非树边,这里的非树边,它的c的值一定要小于最小生成树里面的树边的最小的c值;
这样它才有可能成为那条一直减小的边;然后最后比最小生成树里面的某条边边权来得小;
假如非树边的两个点是u和v;
则如果你要把u-v这条边加到MST里面的话,则必然要在从u到v的路径上选一条边删掉;
删掉的边应是最大的那条边.(记录这条边是什么,删掉之后下次如果还要删的话,得还原)
具体方式就是找到u和v的LCA;
然后从u到lca的路径中找最大w的边;
然后从v到lca的路径总找最大w的边;
取两个边的w的较大者;
然后把这条边删掉;
然后把这条枚举的边加进去;一直减;然后更新答案;
*/
【Number Of WA】
0
【完整代码】
#include <bits/stdc++.h>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define LL long long
#define rep1(i,a,b) for (int i = a;i <= b;i++)
#define rep2(i,a,b) for (int i = a;i >= b;i--)
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define ms(x,y) memset(x,y,sizeof x)
typedef pair<int,int> pii;
typedef pair<LL,LL> pll;
const int dx[9] = {0,1,-1,0,0,-1,-1,1,1};
const int dy[9] = {0,0,0,-1,1,-1,1,-1,1};
const double pi = acos(-1.0);
const int M = 2e5+100;
const int INF = 21e8;
struct abc
{
int w,c,u,v,id;
};
int n,m,f[M],b[M],fa[M][22],d[M][22],deep[M],C,num,bin[22],bb[M],pri[M];
LL tot,ans,S;
abc a[M];
vector <pii> g[M];
int ff(int x)
{
if (f[x]==x)
return x;
else
return f[x] = ff(f[x]);
}
int maxw(int x,int y)
{
if (a[x].w>a[y].w)
return x;
else
return y;
}
void dfs(int x)
{
rep1(i,1,20)
fa[x][i] = fa[fa[x][i-1]][i-1];
rep1(i,1,20)
d[x][i] = maxw(d[x][i-1],d[fa[x][i-1]][i-1]);
for (pii temp:g[x])
{
int y = temp.fi,w = temp.se;
if (y==fa[x][0]) continue;
fa[y][0] = x,d[y][0] = w;
deep[y] = deep[x]+1;
dfs(y);
}
}
int lcq(int x,int y)
{
if (deep[x]<deep[y])
swap(x,y);
//deep[x]>=deep[y];
int temp = deep[x]-deep[y];
rep1(i,0,20)
if (temp&bin[i])
x = fa[x][i];
//deep[x]==deep[y];
rep2(i,20,0)
if (fa[x][i]!=fa[y][i])
x = fa[x][i],y = fa[y][i];
return x==y?x:fa[x][0];
}
int query(int x,int y)
{
//deep[x]>=deep[y]
int temp = deep[x]-deep[y];
int ret = 0;
rep1(i,0,20)
if (temp&bin[i])
{
ret = maxw(ret,d[x][i]);
x = fa[x][i];
}
return ret;
}
int main()
{
//freopen("F:\\rush.txt","r",stdin);
ios::sync_with_stdio(false),cin.tie(0);//scanf,puts,printf not use
//init??????
bin[0] = 1;
rep1(i,1,20) bin[i] = bin[i-1]<<1;
cin >> n >> m;
rep1(i,1,m) cin >> a[i].w;
rep1(i,1,m) cin >> a[i].c;
rep1(i,1,m)
{
cin >> a[i].u >> a[i].v;
a[i].id = i;
}
cin >>S;
sort(a+1,a+1+m,[&](abc a,abc b) {return a.w<b.w;});
C = INF;
rep1(i,1,n) f[i] = i;
rep1(i,1,m)
{
int x = a[i].u,y = a[i].v;
int r1 = ff(x),r2 = ff(y);
if (r1!=r2)
{
f[r1] = r2;
g[x].pb(mp(y,i));
g[y].pb(mp(x,i));
b[i] = 1;
if (a[i].c<C)
{
C = a[i].c;
num = i;
}
tot+=a[i].w;
}
}
ans = tot-S/C;
dfs(1);
int dl = 0,prei = 0;
rep1(i,1,m)
if (a[i].c<C)
{
int x = a[i].u,y = a[i].v;
int z = lcq(x,y);
int tmp = maxw(query(x,z),query(y,z));
LL temp = tot-a[tmp].w+a[i].w;
temp-=S/a[i].c;
if (temp<ans)
{
b[dl] = 1;b[dl=tmp] = 0;
b[prei] = 0;prei = i;
num = i;b[i] = 1;
ans = temp;
}
}
cout << ans << endl;
rep1(i,1,m)
if (b[i])
{
bb[a[i].id] = 1;
if (num==i)
pri[a[i].id] = a[i].w-S/a[i].c;
else
pri[a[i].id] = a[i].w;
}
rep1(i,1,m)
if (bb[i])
cout << i <<' '<<pri[i]<<endl;
return 0;
}
【codeforces 733F】Drivers Dissatisfaction的更多相关文章
- 【codeforces 733F】 Drivers Dissatisfaction
http://codeforces.com/problemset/problem/733/F (题目链接) 题意 给出一张n个点的无向图,每一条变有两个特征值:${w,c}$:分别表示这条边的权值为$ ...
- 【codeforces 415D】Mashmokh and ACM(普通dp)
[codeforces 415D]Mashmokh and ACM 题意:美丽数列定义:对于数列中的每一个i都满足:arr[i+1]%arr[i]==0 输入n,k(1<=n,k<=200 ...
- 【codeforces 707E】Garlands
[题目链接]:http://codeforces.com/contest/707/problem/E [题意] 给你一个n*m的方阵; 里面有k个联通块; 这k个联通块,每个连通块里面都是灯; 给你q ...
- 【codeforces 707C】Pythagorean Triples
[题目链接]:http://codeforces.com/contest/707/problem/C [题意] 给你一个数字n; 问你这个数字是不是某个三角形的一条边; 如果是让你输出另外两条边的大小 ...
- 【codeforces 709D】Recover the String
[题目链接]:http://codeforces.com/problemset/problem/709/D [题意] 给你一个序列; 给出01子列和10子列和00子列以及11子列的个数; 然后让你输出 ...
- 【codeforces 709B】Checkpoints
[题目链接]:http://codeforces.com/contest/709/problem/B [题意] 让你从起点开始走过n-1个点(至少n-1个) 问你最少走多远; [题解] 肯定不多走啊; ...
- 【codeforces 709C】Letters Cyclic Shift
[题目链接]:http://codeforces.com/contest/709/problem/C [题意] 让你改变一个字符串的子集(连续的一段); ->这一段的每个字符的字母都变成之前的一 ...
- 【Codeforces 429D】 Tricky Function
[题目链接] http://codeforces.com/problemset/problem/429/D [算法] 令Si = A1 + A2 + ... + Ai(A的前缀和) 则g(i,j) = ...
- 【Codeforces 670C】 Cinema
[题目链接] http://codeforces.com/contest/670/problem/C [算法] 离散化 [代码] #include<bits/stdc++.h> using ...
随机推荐
- @value 注解获取属性文件中的值
一.属性文件 db.properties name=jack 二.配置文件 applicationContext.xml <!-- 加载配置文件,该节点只能存在一个,所以用 * ,加载所有属性文 ...
- SSM知识巩固2
数据回显 1.springmvc默认对pojo数据进行回显. pojo数据传入controller方法后,springmvc自动将pojo数据放到request域,key等于pojo类型(首字母小写) ...
- Apach配置本地域名
分三步: 一,开启虚拟域名设置,在路径中找到(apache\apache2.4.9\conf)httpd.conf 配置文件, 去掉 " Include conf/extra/httpd-v ...
- 浙大PAT考试1013~1016(最伤的一次。。)
我能说我1016WA了几天都不得最后还是拿别人代码交的么. .. 真心找不到那个神数据.. . 自己把整个程序的流程都画出来了.细致推敲是木有问题的啊... 题目地址:点击打开链接 先从1013開始介 ...
- maven 创建web项目出错
错误的信息: .m2/repository/org/apache/maven/archetypes/maven-archetype-webapp/maven-metadata-central.xml ...
- Nagios监控nginx服务具体过程
1在nginx 服务器上安装nrpe客户端: Nginx的服务须要监控起来.不然万一down了而不及时修复,会影响web应用.例如以下web应用上面启动的nginx后台进程[root@lb-net-2 ...
- 具体图解 Flume介绍、安装配置
写在前面一: 本文总结"Hadoop生态系统"中的当中一员--Apache Flume 写在前面二: 所用软件说明: 一.什么是Apache Flume 官网:Flume is a ...
- #定位系统性能瓶颈# perf
perf是一个基于Linux 2.6+的调优工具,在liunx性能測量抽象出一套适应于各种不同CPU硬件的通用測量方法,其数据来源于比較新的linux内核提供的 perf_event 接口 系统事件: ...
- inconsistent line endings 解决方法
I'm using Unity 3D in combination with Visual Studio 2008 on a Windows 7 64 bit system. When savi ...
- VS2013找不到SDKDDKVer.h
今天在升级vs2010 的project的时候遇到了一个这种问题.提示:找不到SDKDKVer.h 通过查找资料发现,原来是vs版本号之间Windows SDK的路径宏定义不同,有些坑. 网上有人说能 ...