当前位置: 首页 > news >正文

每日一题之宝石组合

问题描述

在一个神秘的森林里,住着一个小精灵名叫小蓝。有一天,他偶然发现了一个隐藏在树洞里的宝藏,里面装满了闪烁着美丽光芒的宝石。这些宝石都有着不同的颜色和形状,但最引人注目的是它们各自独特的 “闪亮度” 属性。每颗宝石都有一个与生俱来的特殊能力,可以发出不同强度的闪光。小蓝共找到了 NN 枚宝石,第 ii 枚宝石的 “闪亮度” 属性值为 HiHi​,小蓝将会从这 NN 枚宝石中选出三枚进行组合,组合之后的精美程度 SS 可以用以下公式来衡量:

S=HaHbHc⋅LCM⁡(Ha,Hb,Hc)LCM⁡(Ha,Hb)⋅LCM⁡(Ha,Hc)⋅LCM⁡(Hb,Hc)S=Ha​Hb​Hc​⋅LCM(Ha​,Hb​)⋅LCM(Ha​,Hc​)⋅LCM(Hb​,Hc​)LCM(Ha​,Hb​,Hc​)​

其中 LCMLCM 表示的是最小公倍数函数。

小蓝想要使得三枚宝石组合后的精美程度 SS 尽可能的高,请你帮他找出精美程度最高的方案。如果存在多个方案 SS 值相同,优先选择按照 HH 值升序排列后字典序最小的方案。

输入格式

第一行包含一个整数 NN 表示宝石个数。

第二行包含 NN 个整数表示 NN 个宝石的 “闪亮度”。

输出格式

输出一行包含三个整数表示满足条件的三枚宝石的 “闪亮度”。

最直接的做法就是暴力,不过只能骗骗分。

#include <iostream>
#include <vector>
#include <algorithm> // for max function
#include <climits>   // for INT_MIN
using namespace std;

// 计算GCD
long long gcd(long long a, long long b) {
    if (b == 0) return a;
    return gcd(b, a % b);
}

// 计算LCM
long long lcm(long long a, long long b) {
    return (a / gcd(a, b)) * b; // 先除后乘,避免溢出
}

// 计算三个数的LCM
long long lcm_three(long long a, long long b, long long c) {
    return lcm(lcm(a, b), c);
}

signed main() {
    long long n;
    cin >> n;
    vector<long long> a(n);
    for (long long i = 0; i < n; i++) {
        cin >> a[i];
    }
    sort(a.begin(),a.end());
    long long max_value = INT_MIN; // 初始化为极小值
    long long result_a = 0, result_b = 0, result_c = 0;

    // 遍历所有三元组
    for (long long i = 0; i < n; i++) {
        for (long long j = i + 1; j < n; j++) {
            for (long long k = j + 1; k < n; k++) {
                // 计算三个数的LCM
                long long lcm1 = lcm_three(a[i], a[j], a[k]);
                long long lcm2=a[i]*a[j]*a[k];
                long long lcm3=lcm(a[i],a[j])*lcm(a[j],a[k])*lcm(a[i],a[k]);
                long long current_lcm=lcm1*lcm2/lcm3;
                // 更新最大值
                if (current_lcm > max_value) {
                    max_value = current_lcm;
                    result_a = a[i];
                    result_b = a[j];
                    result_c = a[k];
                }
            }
        }
    }

    // 输出结果
    cout << result_a << " " << result_b << " " << result_c << endl;
    return 0;
}

1. 公式分析

公式为:

S=HaHbHc⋅LCM(Ha,Hb,Hc)LCM(Ha,Hb)⋅LCM(Ha,Hc)⋅LCM(Hb,Hc)

最小公倍数与最大公约数(GCD)的关系:

LCM(a,b)=a⋅b/GCD(a,b)

步骤 1:将 LCM 转换为 GCD

根据 LCM 和 GCD 的关系:

LCM(Ha,Hb)=Ha⋅HbGCD(Ha,Hb)

LCM(Ha,Hc)=Ha⋅HcGCD(Ha,Hc)

LCM(Hb​,Hc​)=GCD(Hb​,Hc​)Hb​⋅Hc

​​LCM(Ha,Hb,Hc)=Ha⋅Hb⋅HcGCD(Ha,Hb,Hc)

步骤 2:代入公式

将上述表达式代入原公式:

S=HaHbHc⋅Ha⋅Hb⋅HcGCD(Ha,Hb,Hc)(Ha⋅HbGCD(Ha,Hb))⋅(Ha⋅HcGCD(Ha,Hc))⋅(Hb⋅HcGCD(Hb,Hc))

步骤 3:简化公式

将分母和分子展开:

S=HaHbHc⋅HaHbHcGCD(Ha,Hb,Hc)HaHbHaHcHbHcGCD(Ha,Hb)⋅GCD(Ha,Hc)⋅GCD(Hb,Hc)

进一步简化:

S=HaHbHc⋅GCD(Ha,Hb)⋅GCD(Ha,Hc)⋅GCD(Hb,Hc)GCD(Ha,Hb,Hc)⋅HaHbHc

步骤 4:最终简化

分子和分母中的 HaHbHcHa​Hb​Hc​ 可以约去:

S=GCD(Ha,Hb)⋅GCD(Ha,Hc)⋅GCD(Hb,Hc)GCD(Ha,Hb,Hc)


4. 公式的数学意义

最终推导出的公式为:

S=GCD(Ha,Hb)⋅GCD(Ha,Hc)⋅GCD(Hb,Hc)/GCD(Ha,Hb,Hc)

其实最后S就是一个和GCD(Ha,Hb,Hc)正相关的函数

既然S和最大公约数相关,那么我们就去找可能的最大公约数,从最大开始找(贪心的思路),直到找到满足的三元组

#include<stdio.h>
const int h=1e5;
int main(){
    int n;
    scanf("%d",&n);
    int mp[h+1]={0};//初始化宝石闪亮度统计表
    for(int i=0;i<n;i++){
        int t;
        scanf("%d",&t);
        mp[t]++;//统计亮度为t的宝石数量
    }
    //这里我们另辟蹊径,直接枚举精美程度
    for(int i=h;i>=1;i--){//枚举精美程度i
        int ans=0,now=0;//ans表示寻找到了几个宝石,now表示现在数组有几个宝石
        int num[3];//初始化枚举到的宝石
        for(int j=i;j<=h;j+=i){//对于每个精美度i,我们都需要寻找闪亮度为i,2i,3i...的宝石并统计数量
            ans+=mp[j];//把寻找到的宝石数量统计起来
            for(int k=0;k<mp[j]&&now<3;k++){//把统计到的宝石放到数组
                num[now]=j;
                now++;
            }
            if(ans>=3){//如果找到了三个以上的宝石,说明存在三个宝石使其精美度为i
                for(int k=0;k<3;k++){
                    printf("%d ",num[k]);
                }//输出找到的三个宝石
                printf("\n");
                return 0;
            }
        }
    }
}

相关文章:

  • Docker + Vue2 热重载:为什么需要 CHOKIDAR_USEPOLLING=true?
  • 健康饮食,健康早餐
  • 低功耗抄表方案-支持Modbus、DL/T645 及 DL/T698 协议‌电表
  • 简易的微信聊天网页版【项目测试报告】
  • Spring Boot 自动装配深度解析与实践指南
  • Libgdx游戏开发系列教程(2)——接水滴游戏实现
  • 知识周汇|SAP脚本自动化-淋过雨的人更懂得伞的价值
  • Elasticsearch:解锁深度匹配,运用Elasticsearch DSL构建闪电般的高效模糊搜索体验
  • CentOS 7 中安装 Docker和Docker Compose
  • 实战 Elasticsearch:快速上手与深度实践-2.2.3案例:电商订单日志每秒10万条写入优化
  • 基于OFDR的层压陆相页岩油储层中非对称裂缝群传播的分布式光纤监测
  • 可终身授权的外国工具,不限次数使用!PDF转CAD的软件
  • WeakAuras Lua Script TOC
  • .h264/.h265文件 前端直接播放
  • iBeacon数据包全解析:读懂BLE广播包中的定位密码
  • 网页复制小妙招
  • Electron、Tauri及其它跨平台方案终极对比
  • 用不同语言写力扣题的思考:如何选择最适合的编程语言
  • Spark核心之02:常用算子详解
  • 软考高级信息系统项目管理师笔记-第8章项目整合管理
  • 证据公布!菲律宾6人非法登上铁线礁活动
  • 69岁朱自强被查,曾任南京地铁总经理
  • 总有黑眼圈是因为“虚”吗?怎么睡才能改善?
  • 第一集丨《无尽的尽头》值得关注,《榜上佳婿》平平无奇
  • 传染病防治法修订草案提请三审,拟加强医疗机构疾控能力建设
  • “下山虎”张名扬一回合摘下“狮心”:你们再嘘一个给我听听