漫长的补结题报告之路 poj3622
题意(引用):题意:有很多种草,有两个属性:价格和新鲜度;有很多牛,它们都会各自需求一种草,要求是其价格最低为ai,新鲜度最低为bi,且这些牛不希望自己的草和别人的一样。问要满足所有需求的最小花费是多少?
一开始想的都是各种匹配,结果正解是贪心……
应该来说想不到好方法时,不是二分答案就是贪心了吧?
先按新鲜度为第一关键字,价格为第二关键字排序
从最挑剔(新鲜度要求越高)的牛开始考虑,每次选择应当是价格最小并且能能满足当前牛的牧草
可以这样想:当前处理的牛使没选过的牛中最挑剔的,新鲜度能满足当前牛的一定也能满足之前牛
再来考虑价格,因为价格是最低要求,这次价格选择最低限度的,显然会为之后的选择打开更多空间
所以每次每次选择应当是价格最小并且能能满足当前牛的牧草是最优的
朴素模拟复杂度O(nm) ,所以需要各种数据结构优化
考虑到好久没有写平衡树(其实是一开始想到的就是平衡树)
于是就写了splay
var cv,cw,v,w,d,fa:array[..] of longint;
son:array[..,..] of longint;
root,i,j,p,n,m,t,h:longint;
ans:int64;
function succ(x:longint):longint; //找后继
var p:longint;
begin
p:=son[x,];
while son[p,]<> do p:=son[p,];
exit(p);
end; procedure rotate(x,w:longint);
var y:longint;
begin
y:=fa[x];
if fa[y]<> then
begin
if son[fa[y],]=y then son[fa[y],]:=x
else son[fa[y],]:=x;
end;
fa[x]:=fa[y];
son[y,-w]:=son[x,w];
if son[x,w]<> then fa[son[x,w]]:=y;
son[x,w]:=y;
fa[y]:=x;
end; procedure splay(x:longint); //拍习惯了就快了,还是那句话在哪里就往反向转
var y:longint;
begin
while fa[x]<> do
begin
y:=fa[x];
if fa[y]= then
begin
if son[y,]=x then rotate(x,)
else rotate(x,);
end
else begin
if son[fa[y],]=y then
begin
if son[y,]=x then
begin
rotate(y,);
rotate(x,);
end
else begin
rotate(x,);
rotate(x,);
end;
end
else begin
if son[y,]=x then
begin
rotate(x,);
rotate(x,);
end
else begin
rotate(y,);
rotate(x,);
end;
end;
end;
end;
root:=x;
end; procedure delete(x:longint); //删除写得比较丑陋
var p,y,q,u:longint;
begin
y:=fa[x];
if y= then q:=
else if son[y,]=x then q:=
else if son[y,]=x then q:=;
if (son[x,]<>) and (son[x,]<>) then
begin
p:=succ(x);
if (son[p,]<>) and (fa[p]<>x) then
begin
son[fa[p],]:=son[p,];
fa[son[p,]]:=fa[p];
end
else if fa[p]<>x then son[fa[p],]:=;
son[p,]:=son[x,];
fa[son[x,]]:=p;
if fa[p]<>x then
begin
son[p,]:=son[x,];
fa[son[x,]]:=p;
end;
if y<> then son[y,q]:=p;
fa[p]:=y;
splay(p);
end
else begin
if y<> then
begin
if son[x,]<> then
begin
son[y,q]:=son[x,];
fa[son[x,]]:=y;
end;
if son[x,]<> then
begin
son[y,q]:=son[x,];
fa[son[x,]]:=y;
end;
end;
if son[x,]<> then u:=son[x,] else u:=son[x,];
if y= then
begin
root:=u;
fa[u]:=;
end
else splay(y);
end;
dec(h);
fa[x]:=;
son[x,]:=;
son[x,]:=;
d[x]:=;
end; procedure insert(x:longint);
var p:longint;
begin
inc(t);
inc(h);
d[t]:=x;
if h= then
begin
root:=;
fa[t]:=;
end
else begin
p:=root;
repeat
if d[p]>=x then
begin
if son[p,]= then break;
p:=son[p,];
end
else begin
if son[p,]= then break;
p:=son[p,];
end;
until false;
fa[t]:=p;
if d[p]>=x then son[p,]:=t else son[p,]:=t;
splay(t);
end;
end; procedure swap(var a,b:longint);
var c:longint;
begin
c:=a;
a:=b;
b:=c;
end; procedure sortc(l,r:longint);
var i,j,x,y: longint;
begin
i:=l;
j:=r;
x:=cv[(l+r) shr ];
y:=cw[(l+r) shr ];
repeat
while (cv[i]<x) or ((cv[i]=x) and (cw[i]<y)) do inc(i);
while (x<cv[j]) or ((cv[j]=x) and (cw[j]>y)) do dec(j);
if not(i>j) then
begin
swap(cv[i],cv[j]);
swap(cw[i],cw[j]);
inc(i);
j:=j-;
end;
until i>j;
if l<j then sortc(l,j);
if i<r then sortc(i,r);
end; procedure sort(l,r:longint);
var i,j,x,y: longint;
begin
i:=l;
j:=r;
x:=v[(l+r) shr ];
y:=w[(l+r) shr ];
repeat
while (v[i]<x) or ((v[i]=x) and (w[i]<y)) do inc(i);
while (x<v[j]) or ((v[j]=x) and (w[j]>y)) do dec(j);
if not(i>j) then
begin
swap(w[i],w[j]);
swap(v[i],v[j]);
inc(i);
j:=j-;
end;
until i>j;
if l<j then sort(l,j);
if i<r then sort(i,r);
end;
begin
readln(n,m);
for i:= to n do
readln(cw[i],cv[i]);
for i:= to m do
readln(w[i],v[i]);
if m<n then
begin
writeln(-);
halt;
end;
sortc(,n);
sort(,m);
j:=m;
t:=;
root:=;
fillchar(son,sizeof(son),);
fillchar(fa,sizeof(fa),);
for i:=n downto do
begin
while v[j]>=cv[i] do
begin
insert(w[j]);
dec(j);
end;
insert(cw[i]);
p:=succ(t);
if p<> then
begin
ans:=ans+d[p];
delete(t);
delete(p);
end
else begin
ans:=-;
break;
end;
end;
writeln(ans);
end.
话说splay终于写对了还是1Y,好高兴
漫长的补结题报告之路 poj3622的更多相关文章
- 《基于Arm实验箱的国密算法应用》课程设计 结题报告
<基于Arm实验箱的国密算法应用>课程设计 结题报告 小组成员姓名:20155206赵飞 20155220吴思其 20155234昝昕明 指导教师:娄嘉鹏 设计方案 题目要求:基于Arm实 ...
- 《基于Cortex-M4的ucOS-III的应用》课程设计 结题报告
<基于Cortex-M4的ucOS-III的应用>课程设计 结题报告 小组成员姓名:20155211 解雪莹 20155217 杨笛 20155227 辜彦霖 指导教师:娄嘉鹏 一.设计方 ...
- [置顶] 白话最小边覆盖总结--附加 hdu1151结题报告
刚开始看到这个题目的时候就觉得想法很明了,就是不知道如何去匹配... 去网上看了不少人的解题报告,但是对于刚接触“最小边覆盖”的我来说....还是很困难滴....于是自己又开始一如以往学习“最大独立集 ...
- hdu1281结题报告
哎哎...自己刚刚一看到这个题目居然.....什么都想不到...看了一下别人的解题报告说最大匹配...于是就自己开始构思啦... 对于这个棋盘,有K个可以放棋子的位置....那么 首先我们开始可以求出 ...
- 2013山东省ICPC结题报告
A.Rescue The Princess 已知一个等边三角形的两个顶点A.B,求第三个顶点C,A.B.C成逆时针方向. 常规的解题思路就是用已知的两个点列出x,y方程,但这样求出方程的解的表达式比较 ...
- uva401 - Palindromes结题报告
题目地址 : http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...
- 有向图强连通分支的Tarjan算法讲解 + HDU 1269 连通图 Tarjan 结题报告
题目很简单就拿着这道题简单说说 有向图强连通分支的Tarjan算法 有向图强连通分支的Tarjan算法伪代码如下:void Tarjan(u) {dfn[u]=low[u]=++index//进行DF ...
- 2016noipday1t1玩具迷题结题报告
经常读这个代码有益于比赛时想起一些思路.... day1t1,洛谷dalao称之为水题...??然后我去年还是没拿到分,就这个,我还就写了40%的数据,AC到40,然而这不是关键,注释了freopen ...
- 2017 五一 清北学堂 Day1模拟考试结题报告
预计分数:100+50+50 实际分数:5+50+100 =.= 多重背包 (backpack.cpp/c/pas) (1s/256M) 题目描述 提供一个背包,它最多能负载重量为W的物品. 现在给出 ...
随机推荐
- Node.js之【express 安装问题】
经常在全局安装express后,在cmd里面会找不到express命令, 本地模式安装express:'express' 不是内部或外部命令,也不是可运行的程序或批处理文件. 1.先全局安装expre ...
- web页面开发笔记(不断更新)
布局: 1.左右分列2端:使用float:left+float:right,如果一边有多列的话,另一列会对齐最下面那列.所以一般把左右各设一列,一列中再细分多行. 2.div不跟随:使用clear:b ...
- 解决Execwb 导致 ado崩溃的问题。
http://qc.embarcadero.com/wc/qcmain.aspx?d=61255
- Sharepoint delegate control
<?xml version="1.0" encoding="utf-8"?> <Elements xmlns="http://sch ...
- APP中数据加载的6种方式-b
我们看到的APP,往往有着华丽的启动界面,然后就是漫长的数据加载等待,甚至在无网络的时候,整个处于不可用状态.那么我们怎么处理好界面交互中的加载设计,保证体验无缝衔接,保证用户没有漫长的等待感,而可以 ...
- 5.5 准备创建bean
5.5 准备创建bean 我们不可能指望在一个函数中完成一个复杂的逻辑,而且我们跟踪了这么多Spring代码,经历了这么多函数,或多或少也发现了一些规律:一个真正干活的函数其实是以do开头的,比如d ...
- 闭包小demo
var a = (function(){ var c= 0; return function(){ return ++c; } }()); var g = a(); console.log(g); v ...
- 团体程序设计天梯赛-练习集L1-003. 个位数统计
L1-003. 个位数统计 时间限制 400 ms 内存限制 65536 kB 代码长度限制 8000 B 判题程序 Standard 作者 陈越 给定一个k位整数N = dk-1*10k-1 + . ...
- c++实现文本中英文单词和汉字字符的统计
源代码下载:http://download.csdn.net/detail/nuptboyzhb/4987141 1.统计文本中汉字的频数,为后续的文本分类做基础.对于汉字的统计,需要判断读取的是否为 ...
- 扩展DJANGO的LISTVIEW
不用MODEL,不用QUERYSET,而用get_queryset方法来扩展LISTVIEW, 从而实现特定过滤或搜索功能. class DVListView(ListView): template_ ...