Ideone.com

fork download

copy

from itertools import product, combinations_with_replacement
import time
import threading
from collections import defaultdict
import math
 
# 配置参数（可根据需要修改）
TARGET = 168000  # 目标值
BASE_VALUES = [36.5, 41.5, 59,68.5, 74, 91.5]  # 基础系数列表
FLUCTUATION = 1.0  # 系数波动范围
MAX_SOLUTIONS = 3  # 每个组合的最大解数量
SOLVER_TIMEOUT = 180  # 求解超时时间(秒)
THREE_VAR_THRESHOLD = 220000  # 使用三个变量的阈值
PRODUCT_RANGE_THRESHOLD = 148000  # 乘积范围限制阈值
HIGH_TARGET_THRESHOLD = 260000  # 更高目标值阈值
SHOW_PROGRESS = True  # 是否显示进度
MAX_SOLUTIONS_PER_COMB = 100  # 每个组合的最大解数量，用于提前终止
 
def is_valid_product(p):
    """检查单个乘积是否在有效范围内"""
    if TARGET > PRODUCT_RANGE_THRESHOLD:
        if TARGET > HIGH_TARGET_THRESHOLD:
            return p <= 115000  # 目标值超过260000时，仅限制最大值
        else:
            return 74000 <= p <= 115000  # 目标值在148000-260000之间时，使用原范围
    else:
        return True  # 目标值较小时不限制乘积范围
 
def find_single_variable_solutions(values):
    """查找单个数的解（a*x = TARGET）"""
    solutions = []
    for a in values:
        x = TARGET / a
        if x.is_integer() and 1 <= x <= 10000 and is_valid_product(a * x):
            solutions.append((a, x))
            if len(solutions) >= MAX_SOLUTIONS:
                return solutions  # 提前终止
    return solutions
 
def find_two_variable_solutions(values):
    """查找两个变量的解（a*x + b*y = TARGET）"""
    solutions = defaultdict(list)
 
    # 生成所有可能的(a, b)组合
    for a in values:
        for b in values:
            if a == b:  # 避免重复
                continue
 
            # 计算x的有效范围
            min_x = max(1, math.ceil((TARGET - b * 10000) / a))
            max_x = min(math.floor((TARGET - 1) / a), 10000)
 
            if max_x < min_x:
                continue
 
            for x in range(min_x, max_x + 1):
                remainder = TARGET - a * x
 
                # 检查remainder是否在可能的范围内
                if remainder < 1 or remainder > b * 10000:
                    continue
 
                # 检查remainder是否能被b整除
                if remainder % b == 0:
                    y = remainder // b
                    if 1 <= y <= 10000 and is_valid_product(b * y):
                        solutions[(a, b)].append((a, x, b, y))
                        if len(solutions[(a, b)]) >= MAX_SOLUTIONS_PER_COMB:
                            break  # 提前终止当前组合的搜索
 
    return solutions
 
def find_three_variable_solutions(values):
    """优化的三变量求解算法"""
    solutions = defaultdict(list)
 
    # 对系数进行排序，便于剪枝
    sorted_values = sorted(values)
 
    # 预计算每个系数的有效范围
    value_ranges = {}
    for a in sorted_values:
        if TARGET > PRODUCT_RANGE_THRESHOLD:
            if TARGET > HIGH_TARGET_THRESHOLD:
                min_x = max(1, math.ceil(1 / a))  # 取消下限，最小为1
                max_x = min(10000, math.floor(115000 / a))
            else:
                min_x = max(1, math.ceil(74000 / a))
                max_x = min(10000, math.floor(115000 / a))
        else:
            min_x = 1
            max_x = 10000
        value_ranges[a] = (min_x, max_x)
 
    total_combinations = len(sorted_values) * (len(sorted_values) - 1) * (len(sorted_values) - 2) // 6
    processed_combinations = 0
 
    # 三重循环，但添加了更多剪枝条件
    for i, a in enumerate(sorted_values):
        min_x, max_x = value_ranges[a]
 
        # 计算可能的x值数量，决定步长
        x_count = max_x - min_x + 1
        x_step = max(1, x_count // 1000)  # 自适应步长
 
        for x in range(min_x, max_x + 1, x_step):
            ax = a * x
            if not is_valid_product(ax):
                continue
 
            remainder1 = TARGET - ax
            if TARGET > PRODUCT_RANGE_THRESHOLD:
                if TARGET > HIGH_TARGET_THRESHOLD:
                    if remainder1 < 0:
                        continue
                else:
                    if remainder1 < 2 * 74000:
                        continue
 
            # 限制b的选择范围，避免重复组合
            for j in range(i + 1, len(sorted_values)):
                b = sorted_values[j]
 
                if TARGET > PRODUCT_RANGE_THRESHOLD:
                    if TARGET > HIGH_TARGET_THRESHOLD:
                        min_y = max(1, math.ceil(1 / b))  # 取消下限，最小为1
                        max_y = min(10000, math.floor(remainder1 / b))
                    else:
                        min_y = max(1, math.ceil(74000 / b))
                        max_y = min(10000, math.floor((remainder1 - 74000) / b))
                else:
                    min_y = 1
                    max_y = math.floor(remainder1 / b)
 
                if max_y < min_y:
                    continue
 
                # 计算可能的y值数量，决定步长
                y_count = max_y - min_y + 1
                y_step = max(1, y_count // 100)  # 自适应步长
 
                for y in range(min_y, max_y + 1, y_step):
                    by = b * y
                    if not is_valid_product(by):
                        continue
 
                    remainder2 = remainder1 - by
                    if TARGET > PRODUCT_RANGE_THRESHOLD:
                        if TARGET > HIGH_TARGET_THRESHOLD:
                            if remainder2 < 0 or remainder2 > 115000:
                                continue
                        else:
                            if remainder2 < 74000 or remainder2 > 115000:
                                continue
 
                    # 限制c的选择范围
                    found_solution = False
                    for k in range(j + 1, len(sorted_values)):
                        c = sorted_values[k]
 
                        # 检查remainder2是否能被c整除
                        if remainder2 % c != 0:
                            continue
 
                        z = remainder2 // c
                        if 1 <= z <= 10000 and is_valid_product(c * z):
                            key = (a, b, c)
                            solutions[key].append((a, x, b, y, c, z))
                            found_solution = True
                            if len(solutions[key]) >= MAX_SOLUTIONS_PER_COMB:
                                break  # 提前终止当前组合的搜索
 
                    # 如果找到解且达到最大数量，跳出y循环
                    if found_solution and len(solutions[key]) >= MAX_SOLUTIONS_PER_COMB:
                        break  # 跳出y循环
 
                # 进度显示
                if SHOW_PROGRESS and j % 10 == 0:
                    print(f"\r三变量组合进度: {processed_combinations}/{total_combinations} 组", end='')
                processed_combinations += 1
 
    if SHOW_PROGRESS:
        print(f"\r三变量组合进度: {processed_combinations}/{total_combinations} 组 - 完成")
 
    return solutions
 
def find_balanced_solutions(solutions, var_count, num=2):
    """从所有解中筛选出最平衡的解"""
    if var_count == 1 or not solutions:
        return solutions  # 单变量或无解时直接返回
 
    # 当目标值超过220000时，不计算平衡解
    if TARGET > THREE_VAR_THRESHOLD:
        return []
 
    # 计算解的平衡性（变量间差异最小）
    balanced = []
    for sol in solutions:
        vars = sol[1::2]  # 提取x, y, z
        diff = max(vars) - min(vars)
        balanced.append((diff, sol))
 
    # 按差异排序，取前num个
    return [s for _, s in sorted(balanced, key=lambda x: x[0])[:num]]
 
def find_original_solutions(solutions, balanced_solutions, num=3):
    """从剩余解中获取原始顺序的解"""
    if not solutions:
        return []
 
    # 排除已在平衡解中的项
    remaining = [s for s in solutions if s not in balanced_solutions]
    return remaining[:num]
 
def display_solutions(solutions_dict, var_count):
    """优化的解显示函数"""
    if not solutions_dict:
        return
 
    print(f"\n找到 {len(solutions_dict)} 组{var_count}变量解：")
 
    for i, (coeffs, pair_solutions) in enumerate(sorted(solutions_dict.items()), 1):
        # 当目标值超过220000时，不显示平衡解
        if TARGET > THREE_VAR_THRESHOLD:
            all_display = pair_solutions[:MAX_SOLUTIONS]
            print_tag = "[原始解]"
        else:
            # 计算平衡解和原始解
            balanced = find_balanced_solutions(pair_solutions, var_count)
            original = find_original_solutions(pair_solutions, balanced)
            all_display = balanced + original
 
        if var_count == 1:
            a = coeffs
            print(f"\n{i}. 组合: a={a} ({len(pair_solutions)} 个有效解)")
        elif var_count == 2:
            a, b = coeffs
            print(f"\n{i}. 组合: a={a}, b={b} ({len(pair_solutions)} 个有效解)")
        else:  # var_count == 3
            a, b, c = coeffs
            print(f"\n{i}. 组合: a={a}, b={b}, c={c} ({len(pair_solutions)} 个有效解)")
 
        for j, sol in enumerate(all_display, 1):
            if TARGET > THREE_VAR_THRESHOLD:
                tag = print_tag
            else:
                tag = "[平衡解]" if j <= len(balanced) else "[原始解]"
 
            if var_count == 1:
                a, x = sol
                print(f"  {j}. x={x}, a*x={a*x:.1f}, 总和={a*x:.1f} {tag}")
            elif var_count == 2:
                a, x, b, y = sol
                print(f"  {j}. x={x}, y={y}, a*x={a*x:.1f}, b*y={b*y:.1f}, 总和={a*x + b*y:.1f} {tag}")
            else:  # var_count == 3
                a, x, b, y, c, z = sol
                print(f"  {j}. x={x}, y={y}, z={z}, "
                      f"a*x={a*x:.1f}, b*y={b*y:.1f}, c*z={c*z:.1f}, "
                      f"总和={a*x + b*y + c*z:.1f} {tag}")
 
def run_with_timeout(func, args=(), kwargs=None, timeout=SOLVER_TIMEOUT):
    """运行函数并设置超时限制"""
    if kwargs is None:
        kwargs = {}
 
    result = []
    error = []
 
    def wrapper():
        try:
            result.append(func(*args, **kwargs))
        except Exception as e:
            error.append(e)
 
    thread = threading.Thread(target=wrapper)
    thread.daemon = True
    thread.start()
    thread.join(timeout)
 
    if thread.is_alive():
        print(f"警告: {func.__name__} 超时（{timeout}秒），跳过此方法")
        return None
 
    if error:
        raise error[0]
 
    return result[0]
 
def main():
    print(f"目标值: {TARGET}")
 
    # 生成波动后的系数
    FLUCTUATED_VALUES = [round(v - FLUCTUATION, 1) for v in BASE_VALUES]
 
    # 尝试基础系数
    print(f"\n==== 尝试基础系数 ====")
 
    # 检查目标值是否超过阈值
    if TARGET > THREE_VAR_THRESHOLD:
        print(f"目标值 {TARGET} 超过阈值 {THREE_VAR_THRESHOLD}，只尝试三变量解")
        base_solutions = {
            'single': [],
            'two': [],
            'three': run_with_timeout(find_three_variable_solutions, args=(BASE_VALUES,))
        }
 
        # 显示三变量解
        if base_solutions['three'] and len(base_solutions['three']) > 0:
            display_solutions(base_solutions['three'], 3)
            print(f"\n使用基础系数列表，共找到有效解")
            return
        else:
            print("\n基础系数三变量无解，尝试波动系数")
    else:
        # 目标值未超过阈值，按顺序尝试单、双、三变量解
        base_solutions = {
            'single': run_with_timeout(find_single_variable_solutions, args=(BASE_VALUES,)),
            'two': run_with_timeout(find_two_variable_solutions, args=(BASE_VALUES,)),
            'three': []
        }
 
        # 检查是否有解
        has_solution = False
 
        # 显示单变量解
        if base_solutions['single']:
            has_solution = True
            display_solutions({a: [sol] for a, sol in zip(BASE_VALUES, base_solutions['single']) if sol}, 1)
 
        # 显示双变量解
        if base_solutions['two'] and len(base_solutions['two']) > 0:
            has_solution = True
            display_solutions(base_solutions['two'], 2)
 
        # 单变量和双变量都无解时，尝试三变量解
        if not has_solution:
            print(f"\n==== 单变量和双变量无解，尝试三变量解 ====")
            base_solutions['three'] = run_with_timeout(find_three_variable_solutions, args=(BASE_VALUES,))
 
            if base_solutions['three'] and len(base_solutions['three']) > 0:
                has_solution = True
                display_solutions(base_solutions['three'], 3)
 
        # 如果找到解，退出程序
        if has_solution:
            print(f"\n使用基础系数列表，共找到有效解")
            return
 
    # 如果基础系数没有找到解，尝试波动系数
    print(f"\n==== 尝试波动系数 ====")
 
    # 检查目标值是否超过阈值
    if TARGET > THREE_VAR_THRESHOLD:
        print(f"目标值 {TARGET} 超过阈值 {THREE_VAR_THRESHOLD}，只尝试三变量解")
        fluctuated_solutions = {
            'single': [],
            'two': [],
            'three': run_with_timeout(find_three_variable_solutions, args=(FLUCTUATED_VALUES,))
        }
 
        # 显示三变量解
        if fluctuated_solutions['three'] and len(fluctuated_solutions['three']) > 0:
            display_solutions(fluctuated_solutions['three'], 3)
            print(f"\n使用波动系数列表，共找到有效解")
            return
    else:
        # 目标值未超过阈值，按顺序尝试单、双、三变量解
        fluctuated_solutions = {
            'single': run_with_timeout(find_single_variable_solutions, args=(FLUCTUATED_VALUES,)),
            'two': run_with_timeout(find_two_variable_solutions, args=(FLUCTUATED_VALUES,)),
            'three': []
        }
 
        # 重置标志
        has_solution = False
 
        # 显示单变量解
        if fluctuated_solutions['single']:
            has_solution = True
            display_solutions({a: [sol] for a, sol in zip(FLUCTUATED_VALUES, fluctuated_solutions['single']) if sol}, 1)
 
        # 显示双变量解
        if fluctuated_solutions['two'] and len(fluctuated_solutions['two']) > 0:
            has_solution = True
            display_solutions(fluctuated_solutions['two'], 2)
 
        # 单变量和双变量都无解时，尝试三变量解
        if not has_solution:
            print(f"\n==== 单变量和双变量无解，尝试三变量解 ====")
            fluctuated_solutions['three'] = run_with_timeout(find_three_variable_solutions, args=(FLUCTUATED_VALUES,))
 
            if fluctuated_solutions['three'] and len(fluctuated_solutions['three']) > 0:
                has_solution = True
                display_solutions(fluctuated_solutions['three'], 3)
 
        # 如果找到解，退出程序
        if has_solution:
            print(f"\n使用波动系数列表，共找到有效解")
            return
 
    # 如果所有系数集都没有找到解
    print("\n没有找到符合条件的解，即使使用波动后的系数列表。")
 
if __name__ == "__main__":
    start_time = time.time()
    main()
    print(f"\n总耗时: {time.time() - start_time:.2f}秒")

ZnJvbSBpdGVydG9vbHMgaW1wb3J0IHByb2R1Y3QsIGNvbWJpbmF0aW9uc193aXRoX3JlcGxhY2VtZW50CmltcG9ydCB0aW1lCmltcG9ydCB0aHJlYWRpbmcKZnJvbSBjb2xsZWN0aW9ucyBpbXBvcnQgZGVmYXVsdGRpY3QKaW1wb3J0IG1hdGgKCiMg6YWN572u5Y+C5pWw77yI5Y+v5qC55o2u6ZyA6KaB5L+u5pS577yJClRBUkdFVCA9IDE2ODAwMCAgIyDnm67moIflgLwKQkFTRV9WQUxVRVMgPSBbMzYuNSwgNDEuNSwgNTksNjguNSwgNzQsIDkxLjVdICAjIOWfuuehgOezu+aVsOWIl+ihqApGTFVDVFVBVElPTiA9IDEuMCAgIyDns7vmlbDms6LliqjojIPlm7QKTUFYX1NPTFVUSU9OUyA9IDMgICMg5q+P5Liq57uE5ZCI55qE5pyA5aSn6Kej5pWw6YePClNPTFZFUl9USU1FT1VUID0gMTgwICAjIOaxguino+i2heaXtuaXtumXtCjnp5IpClRIUkVFX1ZBUl9USFJFU0hPTEQgPSAyMjAwMDAgICMg5L2/55So5LiJ5Liq5Y+Y6YeP55qE6ZiI5YC8ClBST0RVQ1RfUkFOR0VfVEhSRVNIT0xEID0gMTQ4MDAwICAjIOS5mOenr+iMg+WbtOmZkOWItumYiOWAvApISUdIX1RBUkdFVF9USFJFU0hPTEQgPSAyNjAwMDAgICMg5pu06auY55uu5qCH5YC86ZiI5YC8ClNIT1dfUFJPR1JFU1MgPSBUcnVlICAjIOaYr+WQpuaYvuekuui/m+W6pgpNQVhfU09MVVRJT05TX1BFUl9DT01CID0gMTAwICAjIOavj+S4que7hOWQiOeahOacgOWkp+ino+aVsOmHj++8jOeUqOS6juaPkOWJjee7iOatogoKZGVmIGlzX3ZhbGlkX3Byb2R1Y3QocCk6CiAgICAiIiLmo4Dmn6XljZXkuKrkuZjnp6/mmK/lkKblnKjmnInmlYjojIPlm7TlhoUiIiIKICAgIGlmIFRBUkdFVCA+IFBST0RVQ1RfUkFOR0VfVEhSRVNIT0xEOgogICAgICAgIGlmIFRBUkdFVCA+IEhJR0hfVEFSR0VUX1RIUkVTSE9MRDoKICAgICAgICAgICAgcmV0dXJuIHAgPD0gMTE1MDAwICAjIOebruagh+WAvOi2hei/hzI2MDAwMOaXtu+8jOS7hemZkOWItuacgOWkp+WAvAogICAgICAgIGVsc2U6CiAgICAgICAgICAgIHJldHVybiA3NDAwMCA8PSBwIDw9IDExNTAwMCAgIyDnm67moIflgLzlnKgxNDgwMDAtMjYwMDAw5LmL6Ze05pe277yM5L2/55So5Y6f6IyD5Zu0CiAgICBlbHNlOgogICAgICAgIHJldHVybiBUcnVlICAjIOebruagh+WAvOi+g+Wwj+aXtuS4jemZkOWItuS5mOenr+iMg+WbtAoKZGVmIGZpbmRfc2luZ2xlX3ZhcmlhYmxlX3NvbHV0aW9ucyh2YWx1ZXMpOgogICAgIiIi5p+l5om+5Y2V5Liq5pWw55qE6Kej77yIYSp4ID0gVEFSR0VU77yJIiIiCiAgICBzb2x1dGlvbnMgPSBbXQogICAgZm9yIGEgaW4gdmFsdWVzOgogICAgICAgIHggPSBUQVJHRVQgLyBhCiAgICAgICAgaWYgeC5pc19pbnRlZ2VyKCkgYW5kIDEgPD0geCA8PSAxMDAwMCBhbmQgaXNfdmFsaWRfcHJvZHVjdChhICogeCk6CiAgICAgICAgICAgIHNvbHV0aW9ucy5hcHBlbmQoKGEsIHgpKQogICAgICAgICAgICBpZiBsZW4oc29sdXRpb25zKSA+PSBNQVhfU09MVVRJT05TOgogICAgICAgICAgICAgICAgcmV0dXJuIHNvbHV0aW9ucyAgIyDmj5DliY3nu4jmraIKICAgIHJldHVybiBzb2x1dGlvbnMKCmRlZiBmaW5kX3R3b192YXJpYWJsZV9zb2x1dGlvbnModmFsdWVzKToKICAgICIiIuafpeaJvuS4pOS4quWPmOmHj+eahOino++8iGEqeCArIGIqeSA9IFRBUkdFVO+8iSIiIgogICAgc29sdXRpb25zID0gZGVmYXVsdGRpY3QobGlzdCkKICAgIAogICAgIyDnlJ/miJDmiYDmnInlj6/og73nmoQoYSwgYinnu4TlkIgKICAgIGZvciBhIGluIHZhbHVlczoKICAgICAgICBmb3IgYiBpbiB2YWx1ZXM6CiAgICAgICAgICAgIGlmIGEgPT0gYjogICMg6YG/5YWN6YeN5aSNCiAgICAgICAgICAgICAgICBjb250aW51ZQogICAgICAgICAgICAgICAgCiAgICAgICAgICAgICMg6K6h566XeOeahOacieaViOiMg+WbtAogICAgICAgICAgICBtaW5feCA9IG1heCgxLCBtYXRoLmNlaWwoKFRBUkdFVCAtIGIgKiAxMDAwMCkgLyBhKSkKICAgICAgICAgICAgbWF4X3ggPSBtaW4obWF0aC5mbG9vcigoVEFSR0VUIC0gMSkgLyBhKSwgMTAwMDApCiAgICAgICAgICAgIAogICAgICAgICAgICBpZiBtYXhfeCA8IG1pbl94OgogICAgICAgICAgICAgICAgY29udGludWUKICAgICAgICAgICAgICAgIAogICAgICAgICAgICBmb3IgeCBpbiByYW5nZShtaW5feCwgbWF4X3ggKyAxKToKICAgICAgICAgICAgICAgIHJlbWFpbmRlciA9IFRBUkdFVCAtIGEgKiB4CiAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICMg5qOA5p+lcmVtYWluZGVy5piv5ZCm5Zyo5Y+v6IO955qE6IyD5Zu05YaFCiAgICAgICAgICAgICAgICBpZiByZW1haW5kZXIgPCAxIG9yIHJlbWFpbmRlciA+IGIgKiAxMDAwMDoKICAgICAgICAgICAgICAgICAgICBjb250aW51ZQogICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgIyDmo4Dmn6VyZW1haW5kZXLmmK/lkKbog73ooqti5pW06ZmkCiAgICAgICAgICAgICAgICBpZiByZW1haW5kZXIgJSBiID09IDA6CiAgICAgICAgICAgICAgICAgICAgeSA9IHJlbWFpbmRlciAvLyBiCiAgICAgICAgICAgICAgICAgICAgaWYgMSA8PSB5IDw9IDEwMDAwIGFuZCBpc192YWxpZF9wcm9kdWN0KGIgKiB5KToKICAgICAgICAgICAgICAgICAgICAgICAgc29sdXRpb25zWyhhLCBiKV0uYXBwZW5kKChhLCB4LCBiLCB5KSkKICAgICAgICAgICAgICAgICAgICAgICAgaWYgbGVuKHNvbHV0aW9uc1soYSwgYildKSA+PSBNQVhfU09MVVRJT05TX1BFUl9DT01COgogICAgICAgICAgICAgICAgICAgICAgICAgICAgYnJlYWsgICMg5o+Q5YmN57uI5q2i5b2T5YmN57uE5ZCI55qE5pCc57SiCiAgICAKICAgIHJldHVybiBzb2x1dGlvbnMKCmRlZiBmaW5kX3RocmVlX3ZhcmlhYmxlX3NvbHV0aW9ucyh2YWx1ZXMpOgogICAgIiIi5LyY5YyW55qE5LiJ5Y+Y6YeP5rGC6Kej566X5rOVIiIiCiAgICBzb2x1dGlvbnMgPSBkZWZhdWx0ZGljdChsaXN0KQogICAgCiAgICAjIOWvueezu+aVsOi/m+ihjOaOkuW6j++8jOS+v+S6juWJquaenQogICAgc29ydGVkX3ZhbHVlcyA9IHNvcnRlZCh2YWx1ZXMpCiAgICAKICAgICMg6aKE6K6h566X5q+P5Liq57O75pWw55qE5pyJ5pWI6IyD5Zu0CiAgICB2YWx1ZV9yYW5nZXMgPSB7fQogICAgZm9yIGEgaW4gc29ydGVkX3ZhbHVlczoKICAgICAgICBpZiBUQVJHRVQgPiBQUk9EVUNUX1JBTkdFX1RIUkVTSE9MRDoKICAgICAgICAgICAgaWYgVEFSR0VUID4gSElHSF9UQVJHRVRfVEhSRVNIT0xEOgogICAgICAgICAgICAgICAgbWluX3ggPSBtYXgoMSwgbWF0aC5jZWlsKDEgLyBhKSkgICMg5Y+W5raI5LiL6ZmQ77yM5pyA5bCP5Li6MQogICAgICAgICAgICAgICAgbWF4X3ggPSBtaW4oMTAwMDAsIG1hdGguZmxvb3IoMTE1MDAwIC8gYSkpCiAgICAgICAgICAgIGVsc2U6CiAgICAgICAgICAgICAgICBtaW5feCA9IG1heCgxLCBtYXRoLmNlaWwoNzQwMDAgLyBhKSkKICAgICAgICAgICAgICAgIG1heF94ID0gbWluKDEwMDAwLCBtYXRoLmZsb29yKDExNTAwMCAvIGEpKQogICAgICAgIGVsc2U6CiAgICAgICAgICAgIG1pbl94ID0gMQogICAgICAgICAgICBtYXhfeCA9IDEwMDAwCiAgICAgICAgdmFsdWVfcmFuZ2VzW2FdID0gKG1pbl94LCBtYXhfeCkKICAgIAogICAgdG90YWxfY29tYmluYXRpb25zID0gbGVuKHNvcnRlZF92YWx1ZXMpICogKGxlbihzb3J0ZWRfdmFsdWVzKSAtIDEpICogKGxlbihzb3J0ZWRfdmFsdWVzKSAtIDIpIC8vIDYKICAgIHByb2Nlc3NlZF9jb21iaW5hdGlvbnMgPSAwCiAgICAKICAgICMg5LiJ6YeN5b6q546v77yM5L2G5re75Yqg5LqG5pu05aSa5Ymq5p6d5p2h5Lu2CiAgICBmb3IgaSwgYSBpbiBlbnVtZXJhdGUoc29ydGVkX3ZhbHVlcyk6CiAgICAgICAgbWluX3gsIG1heF94ID0gdmFsdWVfcmFuZ2VzW2FdCiAgICAgICAgCiAgICAgICAgIyDorqHnrpflj6/og73nmoR45YC85pWw6YeP77yM5Yaz5a6a5q2l6ZW/CiAgICAgICAgeF9jb3VudCA9IG1heF94IC0gbWluX3ggKyAxCiAgICAgICAgeF9zdGVwID0gbWF4KDEsIHhfY291bnQgLy8gMTAwMCkgICMg6Ieq6YCC5bqU5q2l6ZW/CiAgICAgICAgCiAgICAgICAgZm9yIHggaW4gcmFuZ2UobWluX3gsIG1heF94ICsgMSwgeF9zdGVwKToKICAgICAgICAgICAgYXggPSBhICogeAogICAgICAgICAgICBpZiBub3QgaXNfdmFsaWRfcHJvZHVjdChheCk6CiAgICAgICAgICAgICAgICBjb250aW51ZQogICAgICAgICAgICAgICAgCiAgICAgICAgICAgIHJlbWFpbmRlcjEgPSBUQVJHRVQgLSBheAogICAgICAgICAgICBpZiBUQVJHRVQgPiBQUk9EVUNUX1JBTkdFX1RIUkVTSE9MRDoKICAgICAgICAgICAgICAgIGlmIFRBUkdFVCA+IEhJR0hfVEFSR0VUX1RIUkVTSE9MRDoKICAgICAgICAgICAgICAgICAgICBpZiByZW1haW5kZXIxIDwgMDoKICAgICAgICAgICAgICAgICAgICAgICAgY29udGludWUKICAgICAgICAgICAgICAgIGVsc2U6CiAgICAgICAgICAgICAgICAgICAgaWYgcmVtYWluZGVyMSA8IDIgKiA3NDAwMDoKICAgICAgICAgICAgICAgICAgICAgICAgY29udGludWUKICAgICAgICAgICAgICAgIAogICAgICAgICAgICAjIOmZkOWItmLnmoTpgInmi6nojIPlm7TvvIzpgb/lhY3ph43lpI3nu4TlkIgKICAgICAgICAgICAgZm9yIGogaW4gcmFuZ2UoaSArIDEsIGxlbihzb3J0ZWRfdmFsdWVzKSk6CiAgICAgICAgICAgICAgICBiID0gc29ydGVkX3ZhbHVlc1tqXQogICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICBpZiBUQVJHRVQgPiBQUk9EVUNUX1JBTkdFX1RIUkVTSE9MRDoKICAgICAgICAgICAgICAgICAgICBpZiBUQVJHRVQgPiBISUdIX1RBUkdFVF9USFJFU0hPTEQ6CiAgICAgICAgICAgICAgICAgICAgICAgIG1pbl95ID0gbWF4KDEsIG1hdGguY2VpbCgxIC8gYikpICAjIOWPlua2iOS4i+mZkO+8jOacgOWwj+S4ujEKICAgICAgICAgICAgICAgICAgICAgICAgbWF4X3kgPSBtaW4oMTAwMDAsIG1hdGguZmxvb3IocmVtYWluZGVyMSAvIGIpKQogICAgICAgICAgICAgICAgICAgIGVsc2U6CiAgICAgICAgICAgICAgICAgICAgICAgIG1pbl95ID0gbWF4KDEsIG1hdGguY2VpbCg3NDAwMCAvIGIpKQogICAgICAgICAgICAgICAgICAgICAgICBtYXhfeSA9IG1pbigxMDAwMCwgbWF0aC5mbG9vcigocmVtYWluZGVyMSAtIDc0MDAwKSAvIGIpKQogICAgICAgICAgICAgICAgZWxzZToKICAgICAgICAgICAgICAgICAgICBtaW5feSA9IDEKICAgICAgICAgICAgICAgICAgICBtYXhfeSA9IG1hdGguZmxvb3IocmVtYWluZGVyMSAvIGIpCiAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICBpZiBtYXhfeSA8IG1pbl95OgogICAgICAgICAgICAgICAgICAgIGNvbnRpbnVlCiAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAjIOiuoeeul+WPr+iDveeahHnlgLzmlbDph4/vvIzlhrPlrprmraXplb8KICAgICAgICAgICAgICAgIHlfY291bnQgPSBtYXhfeSAtIG1pbl95ICsgMQogICAgICAgICAgICAgICAgeV9zdGVwID0gbWF4KDEsIHlfY291bnQgLy8gMTAwKSAgIyDoh6rpgILlupTmraXplb8KICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgZm9yIHkgaW4gcmFuZ2UobWluX3ksIG1heF95ICsgMSwgeV9zdGVwKToKICAgICAgICAgICAgICAgICAgICBieSA9IGIgKiB5CiAgICAgICAgICAgICAgICAgICAgaWYgbm90IGlzX3ZhbGlkX3Byb2R1Y3QoYnkpOgogICAgICAgICAgICAgICAgICAgICAgICBjb250aW51ZQogICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICByZW1haW5kZXIyID0gcmVtYWluZGVyMSAtIGJ5CiAgICAgICAgICAgICAgICAgICAgaWYgVEFSR0VUID4gUFJPRFVDVF9SQU5HRV9USFJFU0hPTEQ6CiAgICAgICAgICAgICAgICAgICAgICAgIGlmIFRBUkdFVCA+IEhJR0hfVEFSR0VUX1RIUkVTSE9MRDoKICAgICAgICAgICAgICAgICAgICAgICAgICAgIGlmIHJlbWFpbmRlcjIgPCAwIG9yIHJlbWFpbmRlcjIgPiAxMTUwMDA6CiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29udGludWUKICAgICAgICAgICAgICAgICAgICAgICAgZWxzZToKICAgICAgICAgICAgICAgICAgICAgICAgICAgIGlmIHJlbWFpbmRlcjIgPCA3NDAwMCBvciByZW1haW5kZXIyID4gMTE1MDAwOgogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbnRpbnVlCiAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICMg6ZmQ5Yi2Y+eahOmAieaLqeiMg+WbtAogICAgICAgICAgICAgICAgICAgIGZvdW5kX3NvbHV0aW9uID0gRmFsc2UKICAgICAgICAgICAgICAgICAgICBmb3IgayBpbiByYW5nZShqICsgMSwgbGVuKHNvcnRlZF92YWx1ZXMpKToKICAgICAgICAgICAgICAgICAgICAgICAgYyA9IHNvcnRlZF92YWx1ZXNba10KICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICMg5qOA5p+lcmVtYWluZGVyMuaYr+WQpuiDveiiq2PmlbTpmaQKICAgICAgICAgICAgICAgICAgICAgICAgaWYgcmVtYWluZGVyMiAlIGMgIT0gMDoKICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbnRpbnVlCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgeiA9IHJlbWFpbmRlcjIgLy8gYwogICAgICAgICAgICAgICAgICAgICAgICBpZiAxIDw9IHogPD0gMTAwMDAgYW5kIGlzX3ZhbGlkX3Byb2R1Y3QoYyAqIHopOgogICAgICAgICAgICAgICAgICAgICAgICAgICAga2V5ID0gKGEsIGIsIGMpCiAgICAgICAgICAgICAgICAgICAgICAgICAgICBzb2x1dGlvbnNba2V5XS5hcHBlbmQoKGEsIHgsIGIsIHksIGMsIHopKQogICAgICAgICAgICAgICAgICAgICAgICAgICAgZm91bmRfc29sdXRpb24gPSBUcnVlCiAgICAgICAgICAgICAgICAgICAgICAgICAgICBpZiBsZW4oc29sdXRpb25zW2tleV0pID49IE1BWF9TT0xVVElPTlNfUEVSX0NPTUI6CiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYnJlYWsgICMg5o+Q5YmN57uI5q2i5b2T5YmN57uE5ZCI55qE5pCc57SiCiAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgIyDlpoLmnpzmib7liLDop6PkuJTovr7liLDmnIDlpKfmlbDph4/vvIzot7Plh7p55b6q546vCiAgICAgICAgICAgICAgICAgICAgaWYgZm91bmRfc29sdXRpb24gYW5kIGxlbihzb2x1dGlvbnNba2V5XSkgPj0gTUFYX1NPTFVUSU9OU19QRVJfQ09NQjoKICAgICAgICAgICAgICAgICAgICAgICAgYnJlYWsgICMg6Lez5Ye6eeW+queOrwogICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAjIOi/m+W6puaYvuekugogICAgICAgICAgICAgICAgaWYgU0hPV19QUk9HUkVTUyBhbmQgaiAlIDEwID09IDA6CiAgICAgICAgICAgICAgICAgICAgcHJpbnQoZiJccuS4ieWPmOmHj+e7hOWQiOi/m+W6pjoge3Byb2Nlc3NlZF9jb21iaW5hdGlvbnN9L3t0b3RhbF9jb21iaW5hdGlvbnN9IOe7hCIsIGVuZD0nJykKICAgICAgICAgICAgICAgIHByb2Nlc3NlZF9jb21iaW5hdGlvbnMgKz0gMQogICAgCiAgICBpZiBTSE9XX1BST0dSRVNTOgogICAgICAgIHByaW50KGYiXHLkuInlj5jph4/nu4TlkIjov5vluqY6IHtwcm9jZXNzZWRfY29tYmluYXRpb25zfS97dG90YWxfY29tYmluYXRpb25zfSDnu4QgLSDlrozmiJAiKQogICAgCiAgICByZXR1cm4gc29sdXRpb25zCgpkZWYgZmluZF9iYWxhbmNlZF9zb2x1dGlvbnMoc29sdXRpb25zLCB2YXJfY291bnQsIG51bT0yKToKICAgICIiIuS7juaJgOacieino+S4reetm+mAieWHuuacgOW5s+ihoeeahOinoyIiIgogICAgaWYgdmFyX2NvdW50ID09IDEgb3Igbm90IHNvbHV0aW9uczoKICAgICAgICByZXR1cm4gc29sdXRpb25zICAjIOWNleWPmOmHj+aIluaXoOino+aXtuebtOaOpei/lOWbngogICAgCiAgICAjIOW9k+ebruagh+WAvOi2hei/hzIyMDAwMOaXtu+8jOS4jeiuoeeul+W5s+ihoeinowogICAgaWYgVEFSR0VUID4gVEhSRUVfVkFSX1RIUkVTSE9MRDoKICAgICAgICByZXR1cm4gW10KICAgIAogICAgIyDorqHnrpfop6PnmoTlubPooaHmgKfvvIjlj5jph4/pl7Tlt67lvILmnIDlsI/vvIkKICAgIGJhbGFuY2VkID0gW10KICAgIGZvciBzb2wgaW4gc29sdXRpb25zOgogICAgICAgIHZhcnMgPSBzb2xbMTo6Ml0gICMg5o+Q5Y+WeCwgeSwgegogICAgICAgIGRpZmYgPSBtYXgodmFycykgLSBtaW4odmFycykKICAgICAgICBiYWxhbmNlZC5hcHBlbmQoKGRpZmYsIHNvbCkpCiAgICAKICAgICMg5oyJ5beu5byC5o6S5bqP77yM5Y+W5YmNbnVt5LiqCiAgICByZXR1cm4gW3MgZm9yIF8sIHMgaW4gc29ydGVkKGJhbGFuY2VkLCBrZXk9bGFtYmRhIHg6IHhbMF0pWzpudW1dXQoKZGVmIGZpbmRfb3JpZ2luYWxfc29sdXRpb25zKHNvbHV0aW9ucywgYmFsYW5jZWRfc29sdXRpb25zLCBudW09Myk6CiAgICAiIiLku47liankvZnop6PkuK3ojrflj5bljp/lp4vpobrluo/nmoTop6MiIiIKICAgIGlmIG5vdCBzb2x1dGlvbnM6CiAgICAgICAgcmV0dXJuIFtdCiAgICAKICAgICMg5o6S6Zmk5bey5Zyo5bmz6KGh6Kej5Lit55qE6aG5CiAgICByZW1haW5pbmcgPSBbcyBmb3IgcyBpbiBzb2x1dGlvbnMgaWYgcyBub3QgaW4gYmFsYW5jZWRfc29sdXRpb25zXQogICAgcmV0dXJuIHJlbWFpbmluZ1s6bnVtXQoKZGVmIGRpc3BsYXlfc29sdXRpb25zKHNvbHV0aW9uc19kaWN0LCB2YXJfY291bnQpOgogICAgIiIi5LyY5YyW55qE6Kej5pi+56S65Ye95pWwIiIiCiAgICBpZiBub3Qgc29sdXRpb25zX2RpY3Q6CiAgICAgICAgcmV0dXJuCiAgICAKICAgIHByaW50KGYiXG7mib7liLAge2xlbihzb2x1dGlvbnNfZGljdCl9IOe7hHt2YXJfY291bnR95Y+Y6YeP6Kej77yaIikKICAgIAogICAgZm9yIGksIChjb2VmZnMsIHBhaXJfc29sdXRpb25zKSBpbiBlbnVtZXJhdGUoc29ydGVkKHNvbHV0aW9uc19kaWN0Lml0ZW1zKCkpLCAxKToKICAgICAgICAjIOW9k+ebruagh+WAvOi2hei/hzIyMDAwMOaXtu+8jOS4jeaYvuekuuW5s+ihoeinowogICAgICAgIGlmIFRBUkdFVCA+IFRIUkVFX1ZBUl9USFJFU0hPTEQ6CiAgICAgICAgICAgIGFsbF9kaXNwbGF5ID0gcGFpcl9zb2x1dGlvbnNbOk1BWF9TT0xVVElPTlNdCiAgICAgICAgICAgIHByaW50X3RhZyA9ICJb5Y6f5aeL6KejXSIKICAgICAgICBlbHNlOgogICAgICAgICAgICAjIOiuoeeul+W5s+ihoeino+WSjOWOn+Wni+inowogICAgICAgICAgICBiYWxhbmNlZCA9IGZpbmRfYmFsYW5jZWRfc29sdXRpb25zKHBhaXJfc29sdXRpb25zLCB2YXJfY291bnQpCiAgICAgICAgICAgIG9yaWdpbmFsID0gZmluZF9vcmlnaW5hbF9zb2x1dGlvbnMocGFpcl9zb2x1dGlvbnMsIGJhbGFuY2VkKQogICAgICAgICAgICBhbGxfZGlzcGxheSA9IGJhbGFuY2VkICsgb3JpZ2luYWwKICAgICAgICAKICAgICAgICBpZiB2YXJfY291bnQgPT0gMToKICAgICAgICAgICAgYSA9IGNvZWZmcwogICAgICAgICAgICBwcmludChmIlxue2l9LiDnu4TlkIg6IGE9e2F9ICh7bGVuKHBhaXJfc29sdXRpb25zKX0g5Liq5pyJ5pWI6KejKSIpCiAgICAgICAgZWxpZiB2YXJfY291bnQgPT0gMjoKICAgICAgICAgICAgYSwgYiA9IGNvZWZmcwogICAgICAgICAgICBwcmludChmIlxue2l9LiDnu4TlkIg6IGE9e2F9LCBiPXtifSAoe2xlbihwYWlyX3NvbHV0aW9ucyl9IOS4quacieaViOinoykiKQogICAgICAgIGVsc2U6ICAjIHZhcl9jb3VudCA9PSAzCiAgICAgICAgICAgIGEsIGIsIGMgPSBjb2VmZnMKICAgICAgICAgICAgcHJpbnQoZiJcbntpfS4g57uE5ZCIOiBhPXthfSwgYj17Yn0sIGM9e2N9ICh7bGVuKHBhaXJfc29sdXRpb25zKX0g5Liq5pyJ5pWI6KejKSIpCiAgICAgICAgCiAgICAgICAgZm9yIGosIHNvbCBpbiBlbnVtZXJhdGUoYWxsX2Rpc3BsYXksIDEpOgogICAgICAgICAgICBpZiBUQVJHRVQgPiBUSFJFRV9WQVJfVEhSRVNIT0xEOgogICAgICAgICAgICAgICAgdGFnID0gcHJpbnRfdGFnCiAgICAgICAgICAgIGVsc2U6CiAgICAgICAgICAgICAgICB0YWcgPSAiW+W5s+ihoeino10iIGlmIGogPD0gbGVuKGJhbGFuY2VkKSBlbHNlICJb5Y6f5aeL6KejXSIKICAgICAgICAgICAgCiAgICAgICAgICAgIGlmIHZhcl9jb3VudCA9PSAxOgogICAgICAgICAgICAgICAgYSwgeCA9IHNvbAogICAgICAgICAgICAgICAgcHJpbnQoZiIgIHtqfS4geD17eH0sIGEqeD17YSp4Oi4xZn0sIOaAu+WSjD17YSp4Oi4xZn0ge3RhZ30iKQogICAgICAgICAgICBlbGlmIHZhcl9jb3VudCA9PSAyOgogICAgICAgICAgICAgICAgYSwgeCwgYiwgeSA9IHNvbAogICAgICAgICAgICAgICAgcHJpbnQoZiIgIHtqfS4geD17eH0sIHk9e3l9LCBhKng9e2EqeDouMWZ9LCBiKnk9e2IqeTouMWZ9LCDmgLvlkow9e2EqeCArIGIqeTouMWZ9IHt0YWd9IikKICAgICAgICAgICAgZWxzZTogICMgdmFyX2NvdW50ID09IDMKICAgICAgICAgICAgICAgIGEsIHgsIGIsIHksIGMsIHogPSBzb2wKICAgICAgICAgICAgICAgIHByaW50KGYiICB7an0uIHg9e3h9LCB5PXt5fSwgej17en0sICIKICAgICAgICAgICAgICAgICAgICAgIGYiYSp4PXthKng6LjFmfSwgYip5PXtiKnk6LjFmfSwgYyp6PXtjKno6LjFmfSwgIgogICAgICAgICAgICAgICAgICAgICAgZiLmgLvlkow9e2EqeCArIGIqeSArIGMqejouMWZ9IHt0YWd9IikKCmRlZiBydW5fd2l0aF90aW1lb3V0KGZ1bmMsIGFyZ3M9KCksIGt3YXJncz1Ob25lLCB0aW1lb3V0PVNPTFZFUl9USU1FT1VUKToKICAgICIiIui/kOihjOWHveaVsOW5tuiuvue9rui2heaXtumZkOWItiIiIgogICAgaWYga3dhcmdzIGlzIE5vbmU6CiAgICAgICAga3dhcmdzID0ge30KICAgIAogICAgcmVzdWx0ID0gW10KICAgIGVycm9yID0gW10KICAgIAogICAgZGVmIHdyYXBwZXIoKToKICAgICAgICB0cnk6CiAgICAgICAgICAgIHJlc3VsdC5hcHBlbmQoZnVuYygqYXJncywgKiprd2FyZ3MpKQogICAgICAgIGV4Y2VwdCBFeGNlcHRpb24gYXMgZToKICAgICAgICAgICAgZXJyb3IuYXBwZW5kKGUpCiAgICAKICAgIHRocmVhZCA9IHRocmVhZGluZy5UaHJlYWQodGFyZ2V0PXdyYXBwZXIpCiAgICB0aHJlYWQuZGFlbW9uID0gVHJ1ZQogICAgdGhyZWFkLnN0YXJ0KCkKICAgIHRocmVhZC5qb2luKHRpbWVvdXQpCiAgICAKICAgIGlmIHRocmVhZC5pc19hbGl2ZSgpOgogICAgICAgIHByaW50KGYi6K2m5ZGKOiB7ZnVuYy5fX25hbWVfX30g6LaF5pe277yIe3RpbWVvdXR956eS77yJ77yM6Lez6L+H5q2k5pa55rOVIikKICAgICAgICByZXR1cm4gTm9uZQogICAgCiAgICBpZiBlcnJvcjoKICAgICAgICByYWlzZSBlcnJvclswXQogICAgCiAgICByZXR1cm4gcmVzdWx0WzBdCgpkZWYgbWFpbigpOgogICAgcHJpbnQoZiLnm67moIflgLw6IHtUQVJHRVR9IikKICAgIAogICAgIyDnlJ/miJDms6LliqjlkI7nmoTns7vmlbAKICAgIEZMVUNUVUFURURfVkFMVUVTID0gW3JvdW5kKHYgLSBGTFVDVFVBVElPTiwgMSkgZm9yIHYgaW4gQkFTRV9WQUxVRVNdCiAgICAKICAgICMg5bCd6K+V5Z+656GA57O75pWwCiAgICBwcmludChmIlxuPT09PSDlsJ3or5Xln7rnoYDns7vmlbAgPT09PSIpCiAgICAKICAgICMg5qOA5p+l55uu5qCH5YC85piv5ZCm6LaF6L+H6ZiI5YC8CiAgICBpZiBUQVJHRVQgPiBUSFJFRV9WQVJfVEhSRVNIT0xEOgogICAgICAgIHByaW50KGYi55uu5qCH5YC8IHtUQVJHRVR9IOi2hei/h+mYiOWAvCB7VEhSRUVfVkFSX1RIUkVTSE9MRH3vvIzlj6rlsJ3or5XkuInlj5jph4/op6MiKQogICAgICAgIGJhc2Vfc29sdXRpb25zID0gewogICAgICAgICAgICAnc2luZ2xlJzogW10sCiAgICAgICAgICAgICd0d28nOiBbXSwKICAgICAgICAgICAgJ3RocmVlJzogcnVuX3dpdGhfdGltZW91dChmaW5kX3RocmVlX3ZhcmlhYmxlX3NvbHV0aW9ucywgYXJncz0oQkFTRV9WQUxVRVMsKSkKICAgICAgICB9CiAgICAgICAgCiAgICAgICAgIyDmmL7npLrkuInlj5jph4/op6MKICAgICAgICBpZiBiYXNlX3NvbHV0aW9uc1sndGhyZWUnXSBhbmQgbGVuKGJhc2Vfc29sdXRpb25zWyd0aHJlZSddKSA+IDA6CiAgICAgICAgICAgIGRpc3BsYXlfc29sdXRpb25zKGJhc2Vfc29sdXRpb25zWyd0aHJlZSddLCAzKQogICAgICAgICAgICBwcmludChmIlxu5L2/55So5Z+656GA57O75pWw5YiX6KGo77yM5YWx5om+5Yiw5pyJ5pWI6KejIikKICAgICAgICAgICAgcmV0dXJuCiAgICAgICAgZWxzZToKICAgICAgICAgICAgcHJpbnQoIlxu5Z+656GA57O75pWw5LiJ5Y+Y6YeP5peg6Kej77yM5bCd6K+V5rOi5Yqo57O75pWwIikKICAgIGVsc2U6CiAgICAgICAgIyDnm67moIflgLzmnKrotoXov4fpmIjlgLzvvIzmjInpobrluo/lsJ3or5XljZXjgIHlj4zjgIHkuInlj5jph4/op6MKICAgICAgICBiYXNlX3NvbHV0aW9ucyA9IHsKICAgICAgICAgICAgJ3NpbmdsZSc6IHJ1bl93aXRoX3RpbWVvdXQoZmluZF9zaW5nbGVfdmFyaWFibGVfc29sdXRpb25zLCBhcmdzPShCQVNFX1ZBTFVFUywpKSwKICAgICAgICAgICAgJ3R3byc6IHJ1bl93aXRoX3RpbWVvdXQoZmluZF90d29fdmFyaWFibGVfc29sdXRpb25zLCBhcmdzPShCQVNFX1ZBTFVFUywpKSwKICAgICAgICAgICAgJ3RocmVlJzogW10KICAgICAgICB9CiAgICAgICAgCiAgICAgICAgIyDmo4Dmn6XmmK/lkKbmnInop6MKICAgICAgICBoYXNfc29sdXRpb24gPSBGYWxzZQogICAgICAgIAogICAgICAgICMg5pi+56S65Y2V5Y+Y6YeP6KejCiAgICAgICAgaWYgYmFzZV9zb2x1dGlvbnNbJ3NpbmdsZSddOgogICAgICAgICAgICBoYXNfc29sdXRpb24gPSBUcnVlCiAgICAgICAgICAgIGRpc3BsYXlfc29sdXRpb25zKHthOiBbc29sXSBmb3IgYSwgc29sIGluIHppcChCQVNFX1ZBTFVFUywgYmFzZV9zb2x1dGlvbnNbJ3NpbmdsZSddKSBpZiBzb2x9LCAxKQogICAgICAgIAogICAgICAgICMg5pi+56S65Y+M5Y+Y6YeP6KejCiAgICAgICAgaWYgYmFzZV9zb2x1dGlvbnNbJ3R3byddIGFuZCBsZW4oYmFzZV9zb2x1dGlvbnNbJ3R3byddKSA+IDA6CiAgICAgICAgICAgIGhhc19zb2x1dGlvbiA9IFRydWUKICAgICAgICAgICAgZGlzcGxheV9zb2x1dGlvbnMoYmFzZV9zb2x1dGlvbnNbJ3R3byddLCAyKQogICAgICAgIAogICAgICAgICMg5Y2V5Y+Y6YeP5ZKM5Y+M5Y+Y6YeP6YO95peg6Kej5pe277yM5bCd6K+V5LiJ5Y+Y6YeP6KejCiAgICAgICAgaWYgbm90IGhhc19zb2x1dGlvbjoKICAgICAgICAgICAgcHJpbnQoZiJcbj09PT0g5Y2V5Y+Y6YeP5ZKM5Y+M5Y+Y6YeP5peg6Kej77yM5bCd6K+V5LiJ5Y+Y6YeP6KejID09PT0iKQogICAgICAgICAgICBiYXNlX3NvbHV0aW9uc1sndGhyZWUnXSA9IHJ1bl93aXRoX3RpbWVvdXQoZmluZF90aHJlZV92YXJpYWJsZV9zb2x1dGlvbnMsIGFyZ3M9KEJBU0VfVkFMVUVTLCkpCiAgICAgICAgICAgIAogICAgICAgICAgICBpZiBiYXNlX3NvbHV0aW9uc1sndGhyZWUnXSBhbmQgbGVuKGJhc2Vfc29sdXRpb25zWyd0aHJlZSddKSA+IDA6CiAgICAgICAgICAgICAgICBoYXNfc29sdXRpb24gPSBUcnVlCiAgICAgICAgICAgICAgICBkaXNwbGF5X3NvbHV0aW9ucyhiYXNlX3NvbHV0aW9uc1sndGhyZWUnXSwgMykKICAgICAgICAKICAgICAgICAjIOWmguaenOaJvuWIsOino++8jOmAgOWHuueoi+W6jwogICAgICAgIGlmIGhhc19zb2x1dGlvbjoKICAgICAgICAgICAgcHJpbnQoZiJcbuS9v+eUqOWfuuehgOezu+aVsOWIl+ihqO+8jOWFseaJvuWIsOacieaViOinoyIpCiAgICAgICAgICAgIHJldHVybgogICAgCiAgICAjIOWmguaenOWfuuehgOezu+aVsOayoeacieaJvuWIsOino++8jOWwneivleazouWKqOezu+aVsAogICAgcHJpbnQoZiJcbj09PT0g5bCd6K+V5rOi5Yqo57O75pWwID09PT0iKQogICAgCiAgICAjIOajgOafpeebruagh+WAvOaYr+WQpui2hei/h+mYiOWAvAogICAgaWYgVEFSR0VUID4gVEhSRUVfVkFSX1RIUkVTSE9MRDoKICAgICAgICBwcmludChmIuebruagh+WAvCB7VEFSR0VUfSDotoXov4fpmIjlgLwge1RIUkVFX1ZBUl9USFJFU0hPTER977yM5Y+q5bCd6K+V5LiJ5Y+Y6YeP6KejIikKICAgICAgICBmbHVjdHVhdGVkX3NvbHV0aW9ucyA9IHsKICAgICAgICAgICAgJ3NpbmdsZSc6IFtdLAogICAgICAgICAgICAndHdvJzogW10sCiAgICAgICAgICAgICd0aHJlZSc6IHJ1bl93aXRoX3RpbWVvdXQoZmluZF90aHJlZV92YXJpYWJsZV9zb2x1dGlvbnMsIGFyZ3M9KEZMVUNUVUFURURfVkFMVUVTLCkpCiAgICAgICAgfQogICAgICAgIAogICAgICAgICMg5pi+56S65LiJ5Y+Y6YeP6KejCiAgICAgICAgaWYgZmx1Y3R1YXRlZF9zb2x1dGlvbnNbJ3RocmVlJ10gYW5kIGxlbihmbHVjdHVhdGVkX3NvbHV0aW9uc1sndGhyZWUnXSkgPiAwOgogICAgICAgICAgICBkaXNwbGF5X3NvbHV0aW9ucyhmbHVjdHVhdGVkX3NvbHV0aW9uc1sndGhyZWUnXSwgMykKICAgICAgICAgICAgcHJpbnQoZiJcbuS9v+eUqOazouWKqOezu+aVsOWIl+ihqO+8jOWFseaJvuWIsOacieaViOinoyIpCiAgICAgICAgICAgIHJldHVybgogICAgZWxzZToKICAgICAgICAjIOebruagh+WAvOacqui2hei/h+mYiOWAvO+8jOaMiemhuuW6j+WwneivleWNleOAgeWPjOOAgeS4ieWPmOmHj+inowogICAgICAgIGZsdWN0dWF0ZWRfc29sdXRpb25zID0gewogICAgICAgICAgICAnc2luZ2xlJzogcnVuX3dpdGhfdGltZW91dChmaW5kX3NpbmdsZV92YXJpYWJsZV9zb2x1dGlvbnMsIGFyZ3M9KEZMVUNUVUFURURfVkFMVUVTLCkpLAogICAgICAgICAgICAndHdvJzogcnVuX3dpdGhfdGltZW91dChmaW5kX3R3b192YXJpYWJsZV9zb2x1dGlvbnMsIGFyZ3M9KEZMVUNUVUFURURfVkFMVUVTLCkpLAogICAgICAgICAgICAndGhyZWUnOiBbXQogICAgICAgIH0KICAgICAgICAKICAgICAgICAjIOmHjee9ruagh+W/lwogICAgICAgIGhhc19zb2x1dGlvbiA9IEZhbHNlCiAgICAgICAgCiAgICAgICAgIyDmmL7npLrljZXlj5jph4/op6MKICAgICAgICBpZiBmbHVjdHVhdGVkX3NvbHV0aW9uc1snc2luZ2xlJ106CiAgICAgICAgICAgIGhhc19zb2x1dGlvbiA9IFRydWUKICAgICAgICAgICAgZGlzcGxheV9zb2x1dGlvbnMoe2E6IFtzb2xdIGZvciBhLCBzb2wgaW4gemlwKEZMVUNUVUFURURfVkFMVUVTLCBmbHVjdHVhdGVkX3NvbHV0aW9uc1snc2luZ2xlJ10pIGlmIHNvbH0sIDEpCiAgICAgICAgCiAgICAgICAgIyDmmL7npLrlj4zlj5jph4/op6MKICAgICAgICBpZiBmbHVjdHVhdGVkX3NvbHV0aW9uc1sndHdvJ10gYW5kIGxlbihmbHVjdHVhdGVkX3NvbHV0aW9uc1sndHdvJ10pID4gMDoKICAgICAgICAgICAgaGFzX3NvbHV0aW9uID0gVHJ1ZQogICAgICAgICAgICBkaXNwbGF5X3NvbHV0aW9ucyhmbHVjdHVhdGVkX3NvbHV0aW9uc1sndHdvJ10sIDIpCiAgICAgICAgCiAgICAgICAgIyDljZXlj5jph4/lkozlj4zlj5jph4/pg73ml6Dop6Pml7bvvIzlsJ3or5XkuInlj5jph4/op6MKICAgICAgICBpZiBub3QgaGFzX3NvbHV0aW9uOgogICAgICAgICAgICBwcmludChmIlxuPT09PSDljZXlj5jph4/lkozlj4zlj5jph4/ml6Dop6PvvIzlsJ3or5XkuInlj5jph4/op6MgPT09PSIpCiAgICAgICAgICAgIGZsdWN0dWF0ZWRfc29sdXRpb25zWyd0aHJlZSddID0gcnVuX3dpdGhfdGltZW91dChmaW5kX3RocmVlX3ZhcmlhYmxlX3NvbHV0aW9ucywgYXJncz0oRkxVQ1RVQVRFRF9WQUxVRVMsKSkKICAgICAgICAgICAgCiAgICAgICAgICAgIGlmIGZsdWN0dWF0ZWRfc29sdXRpb25zWyd0aHJlZSddIGFuZCBsZW4oZmx1Y3R1YXRlZF9zb2x1dGlvbnNbJ3RocmVlJ10pID4gMDoKICAgICAgICAgICAgICAgIGhhc19zb2x1dGlvbiA9IFRydWUKICAgICAgICAgICAgICAgIGRpc3BsYXlfc29sdXRpb25zKGZsdWN0dWF0ZWRfc29sdXRpb25zWyd0aHJlZSddLCAzKQogICAgICAgIAogICAgICAgICMg5aaC5p6c5om+5Yiw6Kej77yM6YCA5Ye656iL5bqPCiAgICAgICAgaWYgaGFzX3NvbHV0aW9uOgogICAgICAgICAgICBwcmludChmIlxu5L2/55So5rOi5Yqo57O75pWw5YiX6KGo77yM5YWx5om+5Yiw5pyJ5pWI6KejIikKICAgICAgICAgICAgcmV0dXJuCiAgICAKICAgICMg5aaC5p6c5omA5pyJ57O75pWw6ZuG6YO95rKh5pyJ5om+5Yiw6KejCiAgICBwcmludCgiXG7msqHmnInmib7liLDnrKblkIjmnaHku7bnmoTop6PvvIzljbPkvb/kvb/nlKjms6LliqjlkI7nmoTns7vmlbDliJfooajjgIIiKQoKaWYgX19uYW1lX18gPT0gIl9fbWFpbl9fIjoKICAgIHN0YXJ0X3RpbWUgPSB0aW1lLnRpbWUoKQogICAgbWFpbigpCiAgICBwcmludChmIlxu5oC76ICX5pe2OiB7dGltZS50aW1lKCkgLSBzdGFydF90aW1lOi4yZn3np5IiKSAgICA=

Success #stdin #stdout 0.06s 10700KB

stdin

copy

Standard input is empty

stdout

copy

目标值: 168000

==== 尝试基础系数 ====

找到 30 组2变量解：

1. 组合: a=36.5, b=41.5 (13 个有效解)
  1. x=2173, y=2137.0, a*x=79314.5, b*y=88685.5, 总和=168000.0 [平衡解]
  2. x=2090, y=2210.0, a*x=76285.0, b*y=91715.0, 总和=168000.0 [平衡解]
  3. x=1509, y=2721.0, a*x=55078.5, b*y=112921.5, 总和=168000.0 [原始解]
  4. x=1592, y=2648.0, a*x=58108.0, b*y=109892.0, 总和=168000.0 [原始解]
  5. x=1675, y=2575.0, a*x=61137.5, b*y=106862.5, 总和=168000.0 [原始解]

2. 组合: a=36.5, b=59 (10 个有效解)
  1. x=1816, y=1724.0, a*x=66284.0, b*y=101716.0, 总和=168000.0 [平衡解]
  2. x=1698, y=1797.0, a*x=61977.0, b*y=106023.0, 总和=168000.0 [平衡解]
  3. x=1462, y=1943.0, a*x=53363.0, b*y=114637.0, 总和=168000.0 [原始解]
  4. x=1580, y=1870.0, a*x=57670.0, b*y=110330.0, 总和=168000.0 [原始解]
  5. x=1934, y=1651.0, a*x=70591.0, b*y=97409.0, 总和=168000.0 [原始解]

3. 组合: a=36.5, b=68.5 (9 个有效解)
  1. x=1600, y=1600.0, a*x=58400.0, b*y=109600.0, 总和=168000.0 [平衡解]
  2. x=1463, y=1673.0, a*x=53399.5, b*y=114600.5, 总和=168000.0 [平衡解]
  3. x=1737, y=1527.0, a*x=63400.5, b*y=104599.5, 总和=168000.0 [原始解]
  4. x=1874, y=1454.0, a*x=68401.0, b*y=99599.0, 总和=168000.0 [原始解]
  5. x=2011, y=1381.0, a*x=73401.5, b*y=94598.5, 总和=168000.0 [原始解]

4. 组合: a=36.5, b=74 (7 个有效解)
  1. x=1588, y=1487.0, a*x=57962.0, b*y=110038.0, 总和=168000.0 [平衡解]
  2. x=1736, y=1414.0, a*x=63364.0, b*y=104636.0, 总和=168000.0 [平衡解]
  3. x=1884, y=1341.0, a*x=68766.0, b*y=99234.0, 总和=168000.0 [原始解]
  4. x=2032, y=1268.0, a*x=74168.0, b*y=93832.0, 总和=168000.0 [原始解]
  5. x=2180, y=1195.0, a*x=79570.0, b*y=88430.0, 总和=168000.0 [原始解]

5. 组合: a=36.5, b=91.5 (6 个有效解)
  1. x=1587, y=1203.0, a*x=57925.5, b*y=110074.5, 总和=168000.0 [平衡解]
  2. x=1770, y=1130.0, a*x=64605.0, b*y=103395.0, 总和=168000.0 [平衡解]
  3. x=1953, y=1057.0, a*x=71284.5, b*y=96715.5, 总和=168000.0 [原始解]
  4. x=2136, y=984.0, a*x=77964.0, b*y=90036.0, 总和=168000.0 [原始解]
  5. x=2319, y=911.0, a*x=84643.5, b*y=83356.5, 总和=168000.0 [原始解]

6. 组合: a=41.5, b=36.5 (13 个有效解)
  1. x=2137, y=2173.0, a*x=88685.5, b*y=79314.5, 总和=168000.0 [平衡解]
  2. x=2210, y=2090.0, a*x=91715.0, b*y=76285.0, 总和=168000.0 [平衡解]
  3. x=1334, y=3086.0, a*x=55361.0, b*y=112639.0, 总和=168000.0 [原始解]
  4. x=1407, y=3003.0, a*x=58390.5, b*y=109609.5, 总和=168000.0 [原始解]
  5. x=1480, y=2920.0, a*x=61420.0, b*y=106580.0, 总和=168000.0 [原始解]

7. 组合: a=41.5, b=59 (8 个有效解)
  1. x=1728, y=1632.0, a*x=71712.0, b*y=96288.0, 总和=168000.0 [平衡解]
  2. x=1610, y=1715.0, a*x=66815.0, b*y=101185.0, 总和=168000.0 [平衡解]
  3. x=1374, y=1881.0, a*x=57021.0, b*y=110979.0, 总和=168000.0 [原始解]
  4. x=1492, y=1798.0, a*x=61918.0, b*y=106082.0, 总和=168000.0 [原始解]
  5. x=1846, y=1549.0, a*x=76609.0, b*y=91391.0, 总和=168000.0 [原始解]

8. 组合: a=41.5, b=68.5 (7 个有效解)
  1. x=1465, y=1565.0, a*x=60797.5, b*y=107202.5, 总和=168000.0 [平衡解]
  2. x=1602, y=1482.0, a*x=66483.0, b*y=101517.0, 总和=168000.0 [平衡解]
  3. x=1328, y=1648.0, a*x=55112.0, b*y=112888.0, 总和=168000.0 [原始解]
  4. x=1739, y=1399.0, a*x=72168.5, b*y=95831.5, 总和=168000.0 [原始解]
  5. x=1876, y=1316.0, a*x=77854.0, b*y=90146.0, 总和=168000.0 [原始解]

9. 组合: a=41.5, b=74 (7 个有效解)
  1. x=1468, y=1447.0, a*x=60922.0, b*y=107078.0, 总和=168000.0 [平衡解]
  2. x=1320, y=1530.0, a*x=54780.0, b*y=113220.0, 总和=168000.0 [平衡解]
  3. x=1616, y=1364.0, a*x=67064.0, b*y=100936.0, 总和=168000.0 [原始解]
  4. x=1764, y=1281.0, a*x=73206.0, b*y=94794.0, 总和=168000.0 [原始解]
  5. x=1912, y=1198.0, a*x=79348.0, b*y=88652.0, 总和=168000.0 [原始解]

10. 组合: a=41.5, b=91.5 (5 个有效解)
  1. x=1398, y=1202.0, a*x=58017.0, b*y=109983.0, 总和=168000.0 [平衡解]
  2. x=1581, y=1119.0, a*x=65611.5, b*y=102388.5, 总和=168000.0 [平衡解]
  3. x=1764, y=1036.0, a*x=73206.0, b*y=94794.0, 总和=168000.0 [原始解]
  4. x=1947, y=953.0, a*x=80800.5, b*y=87199.5, 总和=168000.0 [原始解]
  5. x=2130, y=870.0, a*x=88395.0, b*y=79605.0, 总和=168000.0 [原始解]

11. 组合: a=59, b=36.5 (10 个有效解)
  1. x=1578, y=2052.0, a*x=93102.0, b*y=74898.0, 总和=168000.0 [平衡解]
  2. x=1505, y=2170.0, a*x=88795.0, b*y=79205.0, 总和=168000.0 [平衡解]
  3. x=921, y=3114.0, a*x=54339.0, b*y=113661.0, 总和=168000.0 [原始解]
  4. x=994, y=2996.0, a*x=58646.0, b*y=109354.0, 总和=168000.0 [原始解]
  5. x=1067, y=2878.0, a*x=62953.0, b*y=105047.0, 总和=168000.0 [原始解]

12. 组合: a=59, b=41.5 (8 个有效解)
  1. x=1549, y=1846.0, a*x=91391.0, b*y=76609.0, 总和=168000.0 [平衡解]
  2. x=1466, y=1964.0, a*x=86494.0, b*y=81506.0, 总和=168000.0 [平衡解]
  3. x=968, y=2672.0, a*x=57112.0, b*y=110888.0, 总和=168000.0 [原始解]
  4. x=1051, y=2554.0, a*x=62009.0, b*y=105991.0, 总和=168000.0 [原始解]
  5. x=1134, y=2436.0, a*x=66906.0, b*y=101094.0, 总和=168000.0 [原始解]

13. 组合: a=59, b=68.5 (5 个有效解)
  1. x=1366, y=1276.0, a*x=80594.0, b*y=87406.0, 总和=168000.0 [平衡解]
  2. x=1229, y=1394.0, a*x=72511.0, b*y=95489.0, 总和=168000.0 [平衡解]
  3. x=955, y=1630.0, a*x=56345.0, b*y=111655.0, 总和=168000.0 [原始解]
  4. x=1092, y=1512.0, a*x=64428.0, b*y=103572.0, 总和=168000.0 [原始解]
  5. x=1503, y=1158.0, a*x=88677.0, b*y=79323.0, 总和=168000.0 [原始解]

14. 组合: a=59, b=74 (9 个有效解)
  1. x=1232, y=1288, a*x=72688.0, b*y=95312.0, 总和=168000.0 [平衡解]
  2. x=1306, y=1229, a*x=77054.0, b*y=90946.0, 总和=168000.0 [平衡解]
  3. x=936, y=1524, a*x=55224.0, b*y=112776.0, 总和=168000.0 [原始解]
  4. x=1010, y=1465, a*x=59590.0, b*y=108410.0, 总和=168000.0 [原始解]
  5. x=1084, y=1406, a*x=63956.0, b*y=104044.0, 总和=168000.0 [原始解]

15. 组合: a=59, b=91.5 (4 个有效解)
  1. x=1095, y=1130.0, a*x=64605.0, b*y=103395.0, 总和=168000.0 [平衡解]
  2. x=1278, y=1012.0, a*x=75402.0, b*y=92598.0, 总和=168000.0 [平衡解]
  3. x=912, y=1248.0, a*x=53808.0, b*y=114192.0, 总和=168000.0 [原始解]
  4. x=1461, y=894.0, a*x=86199.0, b*y=81801.0, 总和=168000.0 [原始解]

16. 组合: a=68.5, b=36.5 (8 个有效解)
  1. x=1308, y=2148.0, a*x=89598.0, b*y=78402.0, 总和=168000.0 [平衡解]
  2. x=1235, y=2285.0, a*x=84597.5, b*y=83402.5, 总和=168000.0 [平衡解]
  3. x=797, y=3107.0, a*x=54594.5, b*y=113405.5, 总和=168000.0 [原始解]
  4. x=870, y=2970.0, a*x=59595.0, b*y=108405.0, 总和=168000.0 [原始解]
  5. x=943, y=2833.0, a*x=64595.5, b*y=103404.5, 总和=168000.0 [原始解]

17. 组合: a=68.5, b=41.5 (7 个有效解)
  1. x=1316, y=1876.0, a*x=90146.0, b*y=77854.0, 总和=168000.0 [平衡解]
  2. x=1233, y=2013.0, a*x=84460.5, b*y=83539.5, 总和=168000.0 [平衡解]
  3. x=818, y=2698.0, a*x=56033.0, b*y=111967.0, 总和=168000.0 [原始解]
  4. x=901, y=2561.0, a*x=61718.5, b*y=106281.5, 总和=168000.0 [原始解]
  5. x=984, y=2424.0, a*x=67404.0, b*y=100596.0, 总和=168000.0 [原始解]

18. 组合: a=68.5, b=59 (5 个有效解)
  1. x=1276, y=1366.0, a*x=87406.0, b*y=80594.0, 总和=168000.0 [平衡解]
  2. x=1158, y=1503.0, a*x=79323.0, b*y=88677.0, 总和=168000.0 [平衡解]
  3. x=804, y=1914.0, a*x=55074.0, b*y=112926.0, 总和=168000.0 [原始解]
  4. x=922, y=1777.0, a*x=63157.0, b*y=104843.0, 总和=168000.0 [原始解]
  5. x=1040, y=1640.0, a*x=71240.0, b*y=96760.0, 总和=168000.0 [原始解]

19. 组合: a=68.5, b=74 (4 个有效解)
  1. x=1140, y=1215.0, a*x=78090.0, b*y=89910.0, 总和=168000.0 [平衡解]
  2. x=1288, y=1078.0, a*x=88228.0, b*y=79772.0, 总和=168000.0 [平衡解]
  3. x=844, y=1489.0, a*x=57814.0, b*y=110186.0, 总和=168000.0 [原始解]
  4. x=992, y=1352.0, a*x=67952.0, b*y=100048.0, 总和=168000.0 [原始解]

20. 组合: a=68.5, b=91.5 (3 个有效解)
  1. x=1050, y=1050.0, a*x=71925.0, b*y=96075.0, 总和=168000.0 [平衡解]
  2. x=867, y=1187.0, a*x=59389.5, b*y=108610.5, 总和=168000.0 [平衡解]
  3. x=1233, y=913.0, a*x=84460.5, b*y=83539.5, 总和=168000.0 [原始解]

21. 组合: a=74, b=36.5 (8 个有效解)
  1. x=1268, y=2032.0, a*x=93832.0, b*y=74168.0, 总和=168000.0 [平衡解]
  2. x=1195, y=2180.0, a*x=88430.0, b*y=79570.0, 总和=168000.0 [平衡解]
  3. x=757, y=3068.0, a*x=56018.0, b*y=111982.0, 总和=168000.0 [原始解]
  4. x=830, y=2920.0, a*x=61420.0, b*y=106580.0, 总和=168000.0 [原始解]
  5. x=903, y=2772.0, a*x=66822.0, b*y=101178.0, 总和=168000.0 [原始解]

22. 组合: a=74, b=41.5 (6 个有效解)
  1. x=1198, y=1912.0, a*x=88652.0, b*y=79348.0, 总和=168000.0 [平衡解]
  2. x=1115, y=2060.0, a*x=82510.0, b*y=85490.0, 总和=168000.0 [平衡解]
  3. x=783, y=2652.0, a*x=57942.0, b*y=110058.0, 总和=168000.0 [原始解]
  4. x=866, y=2504.0, a*x=64084.0, b*y=103916.0, 总和=168000.0 [原始解]
  5. x=949, y=2356.0, a*x=70226.0, b*y=97774.0, 总和=168000.0 [原始解]

23. 组合: a=74, b=59 (9 个有效解)
  1. x=1229, y=1306, a*x=90946.0, b*y=77054.0, 总和=168000.0 [平衡解]
  2. x=1170, y=1380, a*x=86580.0, b*y=81420.0, 总和=168000.0 [平衡解]
  3. x=757, y=1898, a*x=56018.0, b*y=111982.0, 总和=168000.0 [原始解]
  4. x=816, y=1824, a*x=60384.0, b*y=107616.0, 总和=168000.0 [原始解]
  5. x=875, y=1750, a*x=64750.0, b*y=103250.0, 总和=168000.0 [原始解]

24. 组合: a=74, b=68.5 (4 个有效解)
  1. x=1215, y=1140.0, a*x=89910.0, b*y=78090.0, 总和=168000.0 [平衡解]
  2. x=1078, y=1288.0, a*x=79772.0, b*y=88228.0, 总和=168000.0 [平衡解]
  3. x=804, y=1584.0, a*x=59496.0, b*y=108504.0, 总和=168000.0 [原始解]
  4. x=941, y=1436.0, a*x=69634.0, b*y=98366.0, 总和=168000.0 [原始解]

25. 组合: a=74, b=91.5 (3 个有效解)
  1. x=1014, y=1016.0, a*x=75036.0, b*y=92964.0, 总和=168000.0 [平衡解]
  2. x=1197, y=868.0, a*x=88578.0, b*y=79422.0, 总和=168000.0 [平衡解]
  3. x=831, y=1164.0, a*x=61494.0, b*y=106506.0, 总和=168000.0 [原始解]

26. 组合: a=91.5, b=36.5 (6 个有效解)
  1. x=984, y=2136.0, a*x=90036.0, b*y=77964.0, 总和=168000.0 [平衡解]
  2. x=911, y=2319.0, a*x=83356.5, b*y=84643.5, 总和=168000.0 [平衡解]
  3. x=619, y=3051.0, a*x=56638.5, b*y=111361.5, 总和=168000.0 [原始解]
  4. x=692, y=2868.0, a*x=63318.0, b*y=104682.0, 总和=168000.0 [原始解]
  5. x=765, y=2685.0, a*x=69997.5, b*y=98002.5, 总和=168000.0 [原始解]

27. 组合: a=91.5, b=41.5 (5 个有效解)
  1. x=953, y=1947.0, a*x=87199.5, b*y=80800.5, 总和=168000.0 [平衡解]
  2. x=870, y=2130.0, a*x=79605.0, b*y=88395.0, 总和=168000.0 [平衡解]
  3. x=621, y=2679.0, a*x=56821.5, b*y=111178.5, 总和=168000.0 [原始解]
  4. x=704, y=2496.0, a*x=64416.0, b*y=103584.0, 总和=168000.0 [原始解]
  5. x=787, y=2313.0, a*x=72010.5, b*y=95989.5, 总和=168000.0 [原始解]

28. 组合: a=91.5, b=59 (4 个有效解)
  1. x=1012, y=1278.0, a*x=92598.0, b*y=75402.0, 总和=168000.0 [平衡解]
  2. x=894, y=1461.0, a*x=81801.0, b*y=86199.0, 总和=168000.0 [平衡解]
  3. x=658, y=1827.0, a*x=60207.0, b*y=107793.0, 总和=168000.0 [原始解]
  4. x=776, y=1644.0, a*x=71004.0, b*y=96996.0, 总和=168000.0 [原始解]

29. 组合: a=91.5, b=68.5 (3 个有效解)
  1. x=913, y=1233.0, a*x=83539.5, b*y=84460.5, 总和=168000.0 [平衡解]
  2. x=776, y=1416.0, a*x=71004.0, b*y=96996.0, 总和=168000.0 [平衡解]
  3. x=639, y=1599.0, a*x=58468.5, b*y=109531.5, 总和=168000.0 [原始解]

30. 组合: a=91.5, b=74 (3 个有效解)
  1. x=1016, y=1014.0, a*x=92964.0, b*y=75036.0, 总和=168000.0 [平衡解]
  2. x=868, y=1197.0, a*x=79422.0, b*y=88578.0, 总和=168000.0 [平衡解]
  3. x=720, y=1380.0, a*x=65880.0, b*y=102120.0, 总和=168000.0 [原始解]

使用基础系数列表，共找到有效解

总耗时: 0.02秒

https://ideone.com/GWD5aa

language:

Python 3 nbc (python 3.7.3)

created:

visibility:

public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!

Discover > Sphere Engine API

Discover > IDE Widget

Choose your language