题目源自>https://contest.xinyoudui.com/contest/236/problem/1277
题目
出题屡屡令人不爽,尤其是上传数据传不上、不让传时。
小 M 出了一道题,试图上传数据,然而不幸的是数据可能会爆上传空间限制,于是就先挑了几个必须上传的测试点。
简单来说,数据分为 n 个 Subtask,第 i 个 Subtask 共造了p𝑖个测试点。为了保证数据强度,必须上传至少wi个测试点。然而,由于上传数据的总空间限制,我们总计上传的测试点数不能超过m 个。
换句话说,假设第 i 个 Subtask 上传了qi个数据,那么 qi至少为wi,最多为pi
并且qi 的总和不能超过m。
在这基础上,小 M 希望最小的一个qi/pi尽可能大。你可以帮帮他吗?
输入格式
第一行两个整数 n,m。
接下来 n 行,每行两个整数 pi,wi
输出格式
若无方案满足,请输出 Zip file size too big after extraction.Poor M!。
否则第一行以最简分数形式输出最大化的 qi/pi 最小值(形如 x/y,使得 gcd{𝑥,𝑦}=1;特别地,
1为 1/1,0 为 0/1)。第二行输出你构造的数列 q1 ~ qn。
注意第一行行末不可有多余空格。
样例
input1
10 100
1 1
11 1
45 4
14 5
19 1
26 4
8 1
17 9
19 8
1 0
output1
10/17
1 7 27 9 12 16 5 10 12 1
input2
10 3
2 0
2 0
2 0
2 0
2 0
2 0
2 0
2 0
2 0
2 0
output2
10 3
2 0
2 0
2 0
2 0
2 0
2 0
2 0
2 0
2 0
2 0
数据范围
题解
二分答案算法题解(本题还可用贪心算法完成)
因为题中关键词有最小值的最大,优先考虑二分答案。
在区间
[
0
,
1
)
[0,1)
[0,1)中最小p/q的值进行二分。
整体思路:
1.输入数据:保存数据的同时记录:输入的同时统计p总数存到变量sumP和w总数存到变量sumW
1.特判特殊情况:
(1)sumW>m:无方案
(2)sumP<=m:可以全部上传:q==p,比值为1/1
2.进行二分答案,逼近最小值的最大值
3.处理q/p的化简问题
二分答案部分大致思路:
将p/q的最小值进行二分枚举,那么就是已知最小的p/q,想办法让所有的点上传后其q/p的值大于等于这个最小值。如果按照这个比值进行上传,最后总上传数仍然小于m,则让这个最小值再大一些,再做尝试。如果按照这个比值不能保证所有点上传后,所有的q/p都大于等于该值,则让这个值再小一些
double MINN=1e-10
double l=0,r=1;
while(r-l>MINN){
double mid =(l+r)/2;
for(int i=0;i<n;i++){
对每一个数据考察:让q在[w,p]增长,在总数不超过m的前提下,q/p能否>=mid值:
如果能超过,则考察下一个数据点;否则,mid值过大,退出循环。
}
如果所有点都考察完了,总数量仍小于m,则mid值过小,下次在[mid,r)区间搜索。
如果是考察中途退出,则表示mid值过大,下次在[l,mid)区间搜索。
}
最后考察边界值:0
正解代码:
#include<iostream>
#include<algorithm>
using namespace std;
typedef struct data{
int p;
int q;
int w;
int order;//原本考虑按比值排序后再二分,以提升效率。后来发现不排也可以通过测试。这个变量可以删除不用
}D;
typedef struct ans{
int q;
int order;
}A;
D ds[1005];
A as[1005];
int sumP,sumW;
int gcd(int x,int y);
double MINN=1e-10;
int main(){
int n,m;
cin>>n>>m;
int minP=1,minQ=2;
for(int i=0;i<n;i++){
cin>>ds[i].p>>ds[i].w;
ds[i].q=ds[i].w;
ds[i].order=i;
sumP+=ds[i].p;
sumW+=ds[i].w;
as[i].q=ds[i].q;
as[i].order=i;
if(ds[i].q*minP<ds[i].p*minQ){
minP=ds[i].p;
minQ=ds[i].q;
}
}
if(sumW>m){//如果最低要求都大于m,则没有方案
cout<<"Zip file size too big after extraction.Poor M!";
return 0;
}
if(sumP<=m){//如果pi全部上传都不超过m
cout<<"1/1"<<endl;
for(int i=0;i<n;i++){
cout<<ds[i].p<<" ";
}
return 0;
}
if(sumW==m){//这个特判不写,下面的二分也能分出结果。感觉写了更快。
int t=gcd(minQ,minP);
cout<<minQ/t<<"/"<<minP/t<<endl;
for(int i=0;i<n;i++){
cout<<as[i].q<<" ";
}
return 0;
}
double l=0,r=1.0;
int last=m-sumW;
while(r-l>MINN){
double mid=(l+r)/2;
last=m-sumW;
bool canGo=true;
for(int i=0;i<n;i++){
ds[i].q=ds[i].w;
}
for(int i=0;i<n;i++){
if(ds[i].q>=ds[i].p*mid)continue;
while(ds[i].q<ds[i].p*mid){
if(last>0){
ds[i].q++;
last--;
}
else break;
}
if(ds[i].q<ds[i].p*mid){
r=mid;
canGo=false;
break;
}
}
if(canGo){
//更新minQ,minP
//更新as[]
minQ=2,minP=1;
for(int i=0;i<n;i++){
as[i].q=ds[i].q;
if(ds[i].q*minP<ds[i].p*minQ){//因为浮点数存在精度问题,能不用除法就不用。
//所以将不等式q/p<minQ/minP转化
minP=ds[i].p;
minQ=ds[i].q;
}
}
if(last>0)l=mid;
if(last==0)break;
}
}
int t=gcd(minQ,minP);
cout<<minQ/t<<"/"<<minP/t<<endl;//最简分数
for(int i=0;i<n;i++){
cout<<as[i].q<<" ";
}
return 0;
}
int gcd(int x,int y){
if(x==0)return y;
while(y%x!=0){
int t=y%x;
y=x;
x=t;
}
return x;
}
贪心题解提示
本题优先队列的用法:
#include<queue>
struct Data{
int p;
int q;
int w;
int order;
double fen;
bool operator < (const Data x)const{//重载运算符<,使得优先队列比较依据为变量fen;默认的大根优先队列变成小根优先队列。即top为fen值最小的数据点
return fen>x.fen;
}
};
int main(){
priority_queue<Data> qu;
//常用函数:
//qu.top() 返回队列中的最小值
//qu.pop() 删除队列中的最小值
//qu.push()插入一个元素
/*
大致思路:
将所有的data push进队列
last 保存可添加数据点的数量
while(last>0){
Data d=qu.top();//取出top,即最小的fen
qu.pop();//弹出top元素
d.q++;
last--;
qu.push(d);//处理好的数据点重新push进队
}
*/
}
```