题意

求长度为n的01串中1占总长(大于L)的比例最大的一个子串起点和终点。

分析

前缀和s[i]保存前i个数有几个1,[j+1,i] 这段区间1的比例就是(s[i]-s[j])/(i-j),于是问题转换为找斜率最大的两个点。

如图,加入j时,就要去掉b1、b2,才能维护斜率的单调递增。

以队列里的点做起点,i 结尾的线段斜率最大的是 i和队列里点组成的下凹线的切线。切点前的点就不会再用到了,因为i后面的点和他们的斜率也将不如和这个切点的斜率。

数形结合,斜率优化,单调队列。

代码

#include<deque>
#include<cstdio> using namespace std; deque<int> q; int s[];
int ansl,ansr; int great(int a,int b,int c,int d)//求ab斜率是否大于cd斜率
{
return (s[a]-s[b])*(c-d) - (s[c]-s[d])*(a-b);
} int main()
{
int t,L,n,a;
scanf("%d",&t);
while(t--)
{
scanf("%d%d ",&n,&L);
char cc;
for(int i=;i<=n;i++){
cc=getchar();
s[i]=s[i-]+cc-'';//这样读才不会超时
}
q.clear();
ansl=;
ansr=L;
for(int i=L; i<=n; i++)//以i做线段右端点
{
int j=i-L;
while(q.size()>)//把j作为线段左端点,加入单调队列
{//(单调指斜率单调递增)
int b1=q[q.size()-];//倒数第1个
int b2=q[q.size()-];//倒数第2个
if(great(b1,b2,j,b1)>)//如果b1b2斜率比jb1更大
q.pop_back();//弹出b1,维护单调性
else break;//
}
q.push_back(j);//j入队
while(q.size()>)//去掉队前头不优的起点(线段左端点)
{//因为和i斜率最大的是下凹线的切点,切点前的点不优
if(great(i,q[],i,q[])<=)//以i做终点 q[0]i的斜率小于q[1]i的斜率
q.pop_front();//就弹出队头
else break;
}
int tmp=great(i,q[],ansr,ansl);//i和切点的斜率,也就是最大斜率
if(tmp> || tmp== && i-q[]<ansr-ansl)
{
ansl=q[];//左端点更新
ansr=i;//右端点更新
}
}
printf("%d %d\n",ansl+,ansr);//存的是左端点右1的点,故输出时+1
}
return ;
}

  

【UVA 1451】Average的更多相关文章

  1. 【巧妙算法系列】【Uva 11464】 - Even Parity 偶数矩阵

    偶数矩阵(Even Parity, UVa 11464) 给你一个n×n的01矩阵(每个元素非0即1),你的任务是把尽量少的0变成1,使得每个元素的上.下.左.右的元素(如果存在的话)之和均为偶数.比 ...

  2. 【贪心+中位数】【UVa 11300】 分金币

    (解方程建模+中位数求最短累积位移) 分金币(Spreading the Wealth, UVa 11300) 圆桌旁坐着n个人,每人有一定数量的金币,金币总数能被n整除.每个人可以给他左右相邻的人一 ...

  3. 【UVa 10881】Piotr's Ants

    Piotr's Ants Porsition:Uva 10881 白书P9 中文改编题:[T^T][FJUT]第二届新生赛真S题地震了 "One thing is for certain: ...

  4. 【UVa 116】Unidirectional TSP

    [Link]:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...

  5. 【UVa 1347】Tour

    [Link]:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...

  6. 【UVA 437】The Tower of Babylon(记忆化搜索写法)

    [题目链接]:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...

  7. 【uva 1025】A Spy in the Metro

    [题目链接]:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...

  8. 【Uva 11584】Partitioning by Palindromes

    [Link]:https://cn.vjudge.net/contest/170078#problem/G [Description] 给你若干个只由小写字母组成的字符串; 问你,这个字符串,最少能由 ...

  9. 【Uva 11400】Lighting System Design

    [Link]: [Description] 你要构建一个供电系统; 给你n种灯泡来构建这么一个系统; 每种灯泡有4个参数 1.灯泡的工作电压 2.灯泡的所需的电源的花费(只要买一个电源就能供这种灯泡的 ...

随机推荐

  1. 用Navicat更新数据库表中的某一字段

    最近需要在A表中根据B表的某一值来进行排序输出,无奈SQL技术不够,不知道怎么连接才能达到目标,于是想到在A表中添加B表的目标值字段,然后通过更新A表从而使A表有目标字段,进而进行排名....够不够纠 ...

  2. WWW压缩解压缩

    unity的WWW参考文档:http://game.ceeger.com/Script/WWW/WWW.html 在unity中把资源打包成Assetbundle其实把资源通过 LZMA 压缩成二进制 ...

  3. java10-2 toString()方法

    public String toString():返回该对象的字符串表示. Integer类下的一个静态方法: public static String toHexString(int i):把一个整 ...

  4. cotangent Laplacian

    几何网格处理经常用到 cotangent laplacian矩阵.前几天把这个功能整合到我的Maya 转 Matlab插件了. 这里发一个利用cotangent laplacian计算特征向量并显示的 ...

  5. 或得的一个div是变量时

  6. mousewheel 模拟滚动

    div{ box-sizing:border-box; } .father{ width:500px; height:400px; margin:auto; margin-top: 50px; bor ...

  7. salt-ssh

    1.生产---测试--开发 salt-runners salt-run manage.status#显示当前up或down的minion salt-run manage.down salt-run m ...

  8. C语言 详解多级指针与指针类型的关系

    //V推论①:指针变量的步长只与‘指针变量的值’的类型有关(指针的值的类型 == 指针指向数据的类型) //指针类型跟指针的值有关,指针是占据4个字节大小的内存空间,但是指针的类型却是各不相同的 // ...

  9. C语言 预处理一(文件包含--#include)

    //预处理命令不需要在结束末尾加":" //#inlude可以包含任意类型的文件 //#inlude 将一个源文件的全部内容包含到另一个源文件中,成为它的一个部分,文件包含的一般格 ...

  10. WorldWind源码剖析系列:WorldWind如何确定与视点相关的地形数据的LOD层级与范围

    1.WorldWind如何确定与视点相关的地形数据的LOD层级与范围? 问题描述:WW中是如何判断LOD层次的呢,即在什么情况下获得哪一层级的数据?是否只通过相机视点的高度进行判断? 问题切入:要解决 ...