UVa OJ 1451 - Average

原创 2016年08月29日 10:25:11

UVa OJ 1451 - Average

Problem

给定一个长度为n的01串,选一个长度至少为L的连续子串,使得子串中数字的平均值最 大。如果有多解,子串长度应尽量小;如果仍有多解,起点编号尽量小。序列中的字符编号 为1~n,因此[1,n]就是完整的字符串。1≤n≤100000,1≤L≤1000。

例如,对于如下长度为17的序列00101011011011010,如果L=7,最大平均值为6/8(子 序列为[7,14],其长度为8);如果L=5,子序列[7,11]的平均值最大,为4/5。

Input

Your program is to read from standard input. The input consists of T test cases. The number of test cases T is given in the first line of the input. Each test case starts with a line containing two integers n (1 ≤ n ≤ 100, 000) and L (1 ≤ L ≤ 1, 000) which are the length of a binary sequence and a length lower bound, respectively. In the next line, a string, binary sequence, of length n is given.

Output

Your program is to write to standard output. Print the starting and ending index of the subsequence.

Sample Input

2
17 5
00101011011011010
20 4
11100111100111110000

Sample Output

7 11
6 9

Solution

说来惭愧,这道题目我WA了快一面了。一开始看错题目,提交了好多次都错了,自己还没反应过来,一个劲傻傻地改。后来发现的时候,已经没有做下去的兴致了。

不过最终还是把题目完成了的。总的来说,这道题目还是让我有些收获的。题目的解题思路是通过将数列转化成坐标轴上的图像,平均值这时候也就变成了斜率了。然后一个个点去维护下凸函数的单调数列,找到最优解。

一开始我的cntAverage()这个函数只有三个参数,使用的是除法,提交之后花了0.1s,自己感觉慢了。于是把除法改成了现在的乘法,时间缩短一半,变成了0.05s。可见乘法和除法在计算的效率上还是相差很多的。

#include <iostream>
#include <string>
using namespace std;

const int maxn = 100005;
int n, L, start, ending,temp;
double maxd;
int DNA[maxn],cav[maxn];
string str;

inline int cntAverage(int L, int r, int LL, int rr, int DNA[]){
    return (DNA[r]-DNA[L])*(rr-LL) - (DNA[rr]-DNA[LL])*(r-L);
}

inline void changePoint(int pits,int bump){
    double temp;
    temp = cntAverage(pits, bump, start-1, ending, DNA);
    if (temp < 0) return;
    if (temp || (bump-pits) < (ending-start+1)){
        maxd = temp;
        start = pits + 1;
        ending = bump;
    }
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);

    int cas;
    cin >> cas;
    while (cas--){
        int find = 0, cor = -1;
        cin >> n >> L; cin.get();
        getline(cin, str);
        for (int i = 1; i <= n; ++i)
            DNA[i] = DNA[i-1] + (str[i-1] == '1');
        maxd = DNA[L] / L, start = 1, ending = L;
        for (int i = L; i <= n; ++i){
            temp = i - L;
            while (find < cor && cntAverage(cav[cor], temp, cav[cor - 1], cav[cor], DNA) <= 0)
                --cor;
            cav[++cor] = temp;
            while (find < cor && cntAverage(cav[find], i, cav[find + 1], i, DNA) <= 0)
                ++find;
            changePoint(cav[find],i);
        }
        cout << start << ' ' << ending << '\n';
    }
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

UVa 1451:Average(数形结合)

题目链接:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=844&page=show_pr...

UVa1451 - Average

UVa1451 - AverageA DNA sequence consists of four letters, A, C, G, and T. The GC-ratio of a DNA sequ...

uva 1451 Average 斜率优化

题意:给出一个长度为n的字符串,要求找到长度至少为k的连续子串的最大平均值。n #include #include #include #include #include #in...

1451 - Average(数形结合)

该题表面让我们求一个字符串的问题,但是却可以转化成求斜率的问题, 紫书上已经说的很清楚了,我这里就不再赘述  。 代码如下 : #include using namespace std; cons...

1451 - Average

A DNA sequenceconsists of four letters, A, C, G, and T. The GC-ratio of a DNA sequence is thenumber ...

1451 - Average 快速求平均值

如何快速求取一段区间的平均值 用前缀的思想来看 非常简单 但是 本题题意要求的是 大于等于一段长度的区间的平均值的最大值 而且给出的数据范围很大 O(n*L)的直白比较算法 用于解决此问题不合适 这种...

九度OJ 1451 不容易系列之一 -- 动态规划

题目地址:http://ac.jobdu.com/problem.php?pid=1451 题目描述: 大家常常感慨,要做好一件事情真的不容易,确实,失败比成功容易多了! 做好“一件”事情尚且...

九度OJ 1451-1460(10/10)

1451#include int main() { int i, j, n; long long f[21]; while(scanf("%d", &n) != EOF) ...

九度OJ-1451:不容易系列之一

错排问题。思考过程写在书上了。 题目描述: 大家常常感慨,要做好一件事情真的不容易,确实,失败比成功容易多了! 做好“一件”事情尚且不易,若想永远成功而总从不失败,那更是难上加难了,...

uva1451 平均值 数形结合

1. 数形结合是一种便于分析的手段,可用于表达多个变量间的函数关系. 2. 例如利用斜率来表示平均数(构造两个变量和第三个变量的关系)。 3. 比较斜率时将除法转换为乘法。 4. 通过i和i+1 j-...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)