这个算法(Big-O)的时间复杂度是多少?

时间:2013-12-31 20:17:58

标签: c algorithm time big-o

我仍在尝试理解BigO符号和时间复杂度,但是,我真的不确定这个算法的时间复杂度(我的代码)。

// 03_BeaverConstructions.c
// Created for FIKS on 28/12/2013 by Dominik Hadl
//
// Time complexity: O(N+M)
// Space complexity: O(N)
// ------------------------------------------
// LICENSE (MIT)
// Copyright (c) 2013 Dominik Hadl
//
// Permission is hereby granted, free of charge, to any person obtaining a copy
// of this software and associated documentation files (the "Software"), to deal
// in the Software without restriction, including without limitation the rights
// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
// copies of the Software, and to permit persons to whom the Software is
// furnished to do so, subject to the following conditions:
//
// The above copyright notice and this permission notice shall be included in
// all copies or substantial portions of the Software.
//
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
// THE SOFTWARE.
// ------------------------------------------

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

// ------------------------------------------
// Setup
// ------------------------------------------

#define MAX_PROFILE_LENGTH 10000
#define NO_VALUE -99999

#define BLOCK_VOLUME 1
#define SLOPE_VOLUME 0.5

const char kSlopeRise = '/';
const char kSlopeLower = '\\';
const char kSlopeStay = '_';

// ------------------------------------------
// Structs
// ------------------------------------------

typedef struct
{
    int start_elevation;
    float current_volume;
} Lake;

typedef struct
{
    int location;
    int elevation;
} Peak;

// ------------------------------------------
// Declarations
// ------------------------------------------

int main(int argc, char const *argv[]);
float get_water_volume_of_profile(char const *hill_profile);

// ------------------------------------------
// Main
// ------------------------------------------

int main(int argc, char const *argv[])
{
    // Get the profile
    char hill_profile[MAX_PROFILE_LENGTH + 1];
    fgets(hill_profile, MAX_PROFILE_LENGTH + 1, stdin);

    // Calculate the volume
    float volume = get_water_volume_of_profile(hill_profile);

    // Print it!
    printf("%0.1f\n", volume);

    return 0;
}

// ------------------------------------------
// Calculation
// ------------------------------------------

float get_water_volume_of_profile(char const *hill_profile)
{
    float total_volume = 0; 
    int current_elevation = 0, number_of_peaks = 0, last_peak_index = 0;

    // Get the actual length of the hill profile    
    int profile_length = strlen(hill_profile);

    // Prepare the peaks and lakes in the hill profile
    Peak peaks[profile_length / 2];
    Lake lake = {NO_VALUE, 0};

    // First, get all the peaks
    for (int i = 0; i < profile_length; i++)
    {
        char current_char = hill_profile[i];
        char next_char = hill_profile[i + 1];

        switch (current_char)
        {
            case kSlopeRise:
                current_elevation += 1;
                break;
            case kSlopeLower:
                current_elevation -= 1;
                break;
            case kSlopeStay:
                break;
        }

        if (next_char == '\n')
        {
            peaks[number_of_peaks].location = i + 1;
            peaks[number_of_peaks].elevation = current_elevation;
            number_of_peaks++;
            break;
        }

        if (current_char == kSlopeRise &&
            (next_char == kSlopeLower || next_char == kSlopeStay))
        {
            peaks[number_of_peaks].location = i + 1;
            peaks[number_of_peaks].elevation = current_elevation;
            number_of_peaks++;
        }
    }

    // Now, go through the profile and get the water volume
    current_elevation = 0;
    for (int i = 0; i < profile_length; i++)
    {
        // Get current char and decide what to do
        char current_char = hill_profile[i];
        switch (current_char)
        {
            case kSlopeRise:
            {
                if (lake.start_elevation != NO_VALUE &&
                    lake.start_elevation > current_elevation)
                {
                    lake.current_volume += SLOPE_VOLUME;
                }

                // Increase the elevation
                current_elevation++;

                if (lake.start_elevation == current_elevation)
                {
                    total_volume += lake.current_volume;

                    lake.start_elevation = NO_VALUE;
                    lake.current_volume = 0;
                    break;
                }

                if (lake.start_elevation != NO_VALUE)
                {
                    int elevation_diff = abs(lake.start_elevation - current_elevation);
                    if (elevation_diff > 1)
                    {
                        lake.current_volume += (elevation_diff - 1) * BLOCK_VOLUME;
                    }
                }

                break;
            }
            case kSlopeLower:
            {
                current_elevation--; // Lower the elevation

                // Set elevation where water starts if not already set
                if (lake.start_elevation == NO_VALUE)
                {
                    for (int p = last_peak_index; p < number_of_peaks; p++)
                    {
                        if (peaks[p].elevation >= current_elevation + 1 &&
                            peaks[p].location > i)
                        {
                            lake.start_elevation = current_elevation + 1;
                            last_peak_index = p;
                            break;
                        }
                    }
                    if (lake.start_elevation == NO_VALUE) 
                    {
                        break;
                    }
                }

                lake.current_volume += SLOPE_VOLUME;

                int elevation_diff = abs(lake.start_elevation - current_elevation);
                if (elevation_diff > 1)
                {
                    lake.current_volume += elevation_diff * BLOCK_VOLUME;
                }

                break;
            }
            case kSlopeStay:
            {
                if (lake.start_elevation != NO_VALUE)
                {
                    int elevation_diff = abs(lake.start_elevation - current_elevation);
                    lake.current_volume += elevation_diff * BLOCK_VOLUME;
                }   
                break;
            }
        }
    }

    // Return the total water volume
    return total_volume;
}

我不确定它是否是O(N),但我不这么认为,因为第二个for循环中有一个嵌套循环。但是,它可能不是O(N ^ 2)......更像O((N ^ 2)/ 2)。

有人可以给我建议吗?

2 个答案:

答案 0 :(得分:2)

算法的复杂性为O(n + m),其中n是输入的大小,m是“峰值”的数量,无论这些是什么。

原因是核心算法由两个循环组成,大约运行n次。其中一个循环包含一个内循环。我们需要计算内循环体的执行次数。

当您遇到峰值时会运行内部循环,并且它看起来像循环的 body 执行的总次数大致是您拥有的峰值数。循环嵌套并不重要:对于复杂度计算,正文的迭代次数是重要的。

(通常嵌套循环的迭代计数是相乘而不是相加,因为它在外循环的每次迭代中都是完整执行的,但这不是这里的情况。你逻辑上是从第一个峰值迭代(在内循环中)到最后;注意你跟踪(使用p)你break从外部迭代之间的内部循环中的位置,并在返回到p时从{{1}}开始内循环。)

答案 1 :(得分:1)

我错过了什么吗?我宁愿你刚刚提供了伪代码/或必要的代码行,但看起来它只是两行for循环,而第二个for循环有一个嵌套循环。那将是O(N ^ 2),不是吗?随着数据集的增加,时间复杂度应与数据集的平方成正比...

但是,如果有人想要纠正我,我从未掌握过这个主题。

编辑 - 请参阅评论,了解这是错误的原因..!