我仍在尝试理解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)。
有人可以给我建议吗?
答案 0 :(得分:2)
算法的复杂性为O(n + m)
,其中n
是输入的大小,m
是“峰值”的数量,无论这些是什么。
原因是核心算法由两个循环组成,大约运行n
次。其中一个循环包含一个内循环。我们需要计算内循环体的执行次数。
当您遇到峰值时会运行内部循环,并且它看起来像循环的 body 执行的总次数大致是您拥有的峰值数。循环嵌套并不重要:对于复杂度计算,正文的总迭代次数是重要的。
(通常嵌套循环的迭代计数是相乘而不是相加,因为它在外循环的每次迭代中都是完整执行的,但这不是这里的情况。你逻辑上是从第一个峰值迭代(在内循环中)到最后;注意你跟踪(使用p
)你break
从外部迭代之间的内部循环中的位置,并在返回到p
时从{{1}}开始内循环。)
答案 1 :(得分:1)
我错过了什么吗?我宁愿你刚刚提供了伪代码/或必要的代码行,但看起来它只是两行for循环,而第二个for循环有一个嵌套循环。那将是O(N ^ 2),不是吗?随着数据集的增加,时间复杂度应与数据集的平方成正比...
但是,如果有人想要纠正我,我从未掌握过这个主题。编辑 - 请参阅评论,了解这是错误的原因..!