题意:开始有一个空序列s,一个变量c=0,接着从左往右依次将数组a中的数字放入s的尾部,每放一个数字就检测一次混乱度K,当混乱度k大于M时就清空序列并让c=c+1

   K = Bi * Vi(1<=i<=k(序列总长度)的总和),Bi表示序列中第i小的数字,Vi是给定的非递减的数,输出每次加入序列后的变量c

题解:首先发现当没有清空的时候每次向后增加K都不会减小,所以对于不清空来说K一定是非递减的

   定义左端点L二分寻找每个右端点R,保证现在R是[L,R]第一个大于M的位置,寻找内部是直接排序暴力

   但是会超时,因为可能每个R都离L不远,这样每次减去的数字就很少

   我们考虑另一种二分,首先根据L暴力找到p,使[L,2^P]是第一个大于M的位置,这儿最多是暴力logn次,而且当2^p不是很大时,每次暴力内部都很快([L,2^p]数字少)

   接着二分寻找[L,2^(p-1)]与[L,2^p]中最小的R满足[L,R]第一个大于M,这儿减去的数字少的话也很快,减去的多等等循环的次数就少了

   关键点:循环的次数多则减去的数字多,减去的数字少则循环的次数少;这样可以节约时间

#include<set>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<string>
#include<cstdio>
#include<cstring>
#include<iomanip>
#include<stdlib.h>
#include<iostream>
#include<algorithm>
using namespace std;
#define eps 1E-8
/*注意可能会有输出-0.000*/
#define sgn(x) (x<-eps? -1 :x<eps? 0:1)//x为两个浮点数差的比较,注意返回整型
#define cvs(x) (x > 0.0 ? x+eps : x-eps)//浮点数转化
#define zero(x) (((x)>0?(x):-(x))<eps)//判断是否等于0
#define mul(a,b) (a<<b)
#define dir(a,b) (a>>b)
typedef long long ll;
typedef unsigned long long ull;
const int Inf=<<;
const ll INF=1LL<<;
const double Pi=acos(-1.0);
const int Mod=1e9+;
const int Max=;
ll a[Max],v[Max],temp[Max];
int Pow[],ans[Max];
void Init()
{
Pow[]=;
for(int i=; i<; ++i)
Pow[i]=Pow[i-]*;
return;
}
int Judge(int l,int r,ll m)//判断[l,r]是否大于M
{
ll sum=0LL;
int coun=;
for(int j=l; j<=r; ++j)
{
temp[coun++]=a[j];
}
sort(temp,temp+coun);//此段排序
for(int j=; j<coun; ++j) //暴力计算
{
sum+=temp[j]*v[j];
}
if(sum>m)
return ;
return ;
}
int Dic(int l,int r1,int r2,ll m)//关键:二分范围找到最小的位置满足刚好大于m
{
while(r1<r2)
{
int midr=(r1+r2>>);
if(Judge(l,midr,m))
{
r2=midr;
}
else
{
r1=midr+;
}
}
return r1;
}
int check(int l,int n,ll m)
{
int r1=l,r2=l;//第一个大于M的一定在[r1,r2]之间
for(int i=; ; ++i,r2+=Pow[i]) //关键:枚举logn次,找到最小的i,满足[l,l+^i]刚好大于m
{
if(r2>n)//注意大于n的情况
r2=n;
if(Judge(l,r2,m))
break;
r1=r2+;
if(r2==n)
break;
}
return Dic(l,r1,r2,m);//二分[r1,r2]
}
void Solve(int n,ll m)
{
for(int i=; i<n;) //左端点
{
int j=check(i,n,m);
ans[j]=ans[i==?:i-]+;//第j个位置刚好第一次大于M,则先加值
i=j+;//再清空
}
for(int i=; i<n; ++i)
{
ans[i]=max(ans[i],ans[i-]);//上面只是跟新了拐点的答案
}
return;
}
int main()
{
Init();
int n;
ll m;
while(~scanf("%d %lld",&n,&m))
{
ll u;
for(int i=; i<n; ++i)
{
scanf("%lld",&a[i]);
}
for(int i=; i<n; ++i)
{
scanf("%lld",&v[i]);
}
memset(ans,,sizeof(ans));
Solve(n,m);
for(int i=; i<n; ++i)
printf("%d%c",ans[i],i==n-?'\n':' ');
}
return ;
}

“玲珑杯”ACM比赛 Round #13 B -- 我也不是B(二分排序)的更多相关文章

  1. “玲珑杯”ACM比赛 Round #13 题解&源码

    A 题目链接:http://www.ifrog.cc/acm/problem/1111 分析:容易发现本题就是排序不等式, 将A数组与B数组分别排序之后, 答案即N∑i=1Ai×Bi 此题有坑,反正据 ...

  2. “玲珑杯”ACM比赛 Round #13 B -- 我也不是B,倍增+二分!

    B  我也不是B   这个题做了一下午,比赛两个小时还是没做出来,比完赛才知道要用一个倍增算法确定区间,然后再二分右端点.   题意:定义一个序列的混乱度为累加和:b[i]*v[i],b[i]为这个序 ...

  3. “玲珑杯”ACM比赛 Round #1

    Start Time:2016-08-20 13:00:00 End Time:2016-08-20 18:00:00 Refresh Time:2017-11-12 19:51:52 Public ...

  4. “玲珑杯”ACM比赛 Round #12题解&源码

    我能说我比较傻么!就只能做一道签到题,没办法,我就先写下A题的题解&源码吧,日后补上剩余题的题解&源码吧!                                     A ...

  5. “玲珑杯”ACM比赛 Round #19题解&源码【A,规律,B,二分,C,牛顿迭代法,D,平衡树,E,概率dp】

    A -- simple math problem Time Limit:2s Memory Limit:128MByte Submissions:1599Solved:270 SAMPLE INPUT ...

  6. “玲珑杯”ACM比赛 Round #19 B -- Buildings (RMQ + 二分)

    “玲珑杯”ACM比赛 Round #19 Start Time:2017-07-29 14:00:00 End Time:2017-07-29 16:30:00 Refresh Time:2017-0 ...

  7. “玲珑杯”ACM比赛 Round #18

    “玲珑杯”ACM比赛 Round #18 Start Time:2017-07-15 12:00:00 End Time:2017-07-15 15:46:00 A -- 计算几何你瞎暴力 Time ...

  8. “玲珑杯”ACM比赛 Round #1 题解

    A:DESCRIPTION Eric has an array of integers a1,a2,...,ana1,a2,...,an. Every time, he can choose a co ...

  9. 玲珑杯”ACM比赛 Round #4 1054 - String cut 暴力。学到了扫描的另一种思想

    http://www.ifrog.cc/acm/problem/1054 问删除一个字符后的最小循环节是多少. 比赛的时候想不出,不知道怎么暴力. 赛后看了别人代码才晓得.唉,还以为自己字符串还不错, ...

随机推荐

  1. java怎么计算散列码hashcode?

    转自:https://blog.csdn.net/qq_21430549/article/details/52225801 1.从HashMap说起 我们知道Map以键值对的形式来存储数据.有一点值得 ...

  2. PHP 防止恶意用户快速刷新页面

    <?php /** * 防止快速刷新 * */ session_start(); header("Content-type:text/html;charset=utf-8") ...

  3. Java Concurrency In Practice

    线程安全 定义 A class is thread-safe if it behaves correctly when accessed from multiple threads, regardle ...

  4. scrapy 原理,结构,基本命令,item,spider,selector简述

    原理,结构,基本命令,item,spider,selector简述 原理 (1)结构 (2)运行流程 实操 (1) scrapy命令: 注意先把python安装目录的scripts文件夹添加到环境变量 ...

  5. 手动拼写出来的sp_who结果集

    SELECT SPID = er.session_id  ,STATUS = ses.STATUS  ,[Login] = ses.login_name  ,Host = ses.host_name  ...

  6. export,import ,export default的区别

    首先要知道export,import ,export default是什么 ES6模块主要有两个功能:export和importexport用于对外输出本模块(一个文件可以理解为一个模块)变量的接口i ...

  7. 初识python (一)

    初识Python(一) python2和python3的一些区别 Python2 和 Python3 区别汇总:http://www.cnblogs.com/bigtreei/p/7806288.ht ...

  8. Learn How To Cross Over The Wall

    1.一个proxy的实现 http://blog.codingnow.com/2011/05/xtunnel.html 2.SOCK5 RFC http://www.faqs.org/rfcs/rfc ...

  9. Python(数据库之约束表的关系)

    一.约束 约束条件与数据类型的宽度一样,都是可选参数 作用:用于保证数据的完整性和一致性 主要分为: RIMARY KEY (PK) 标识该字段为该表的主键,可以唯一的标识记录 FOREIGN KEY ...

  10. 用仿ActionScript的语法来编写html5——第五篇,Graphics绘图

    用仿ActionScript的语法来编写html5——第五篇,Graphics绘图 canvas本身就是一个Graphics,可以直接进行绘图在actionscript里面,每个Sprite都有一个G ...