Project Euler_Problem 190_Maximising a Weighted Product_最优化问题,梯度方向

文章讨论了解决条件极值问题中,如何通过迭代调整实数值,使得特定连乘表达式达到最大值的算法,利用梯度和数值优化的方法,通过编程实现求解过程。
摘要由CSDN通过智能技术生成

原题目:

题目大意:条件极值, k个实数加起来=k,问当怎么取值的时候,x_1*x_2^2....x_k^k乘起来最大

解题思路: 条件极值问题,求那个连乘的梯度方向,然后从低的往高的挪值,重复个几十万次,就ok了

代码:

void solve() {
    ll i, j,k,x,y,z,p,q,u,v;
    ll N = 15,NN=1000;
    double a, b, c,d;

    for (i = 2; i <= N; i++) {

        ll T = 200000;
        for (j = 1; j <= i; j++) {
            AA[j] = 1.0;
        }
        d = 1.0;
        while (T) {


            for (k = 1; k <= i; k++) {
                c = k;
                for (p = 1; p <= i; p++) {
                    if (p == k) {
                        if (p == 1)continue;
                        c = c * pow(AA[p], p - 1);
                    }
                    else {
                        c = c * pow(AA[p], p);
                    }
                }
                BB[k] = c;
            }

            a = T / 800000.0;

            double M = BB[1], m = BB[1];
            ll M1 = 1, m2 = 1;
            for (k = 2; k <= i; k++) {
                if (BB[k] > M) {
                    M = BB[k];
                    M1 = k;
                }
                if (BB[k] < m) {
                    m = BB[k];
                    m2 = k;
                }

            }
            a = min(AA[m2], a);
            AA[M1] = AA[M1] + a;
            AA[m2] = AA[m2] - a;
            double sum1 = 1;
            for (k = 1; k <= i; k++) {
                sum1 = sum1 * pow(AA[k], k);
            }  
            if (sum1 > d) {
                d = sum1;
            }
            T--;
        }
        printf("%lf\n", d);
        ans1 = ans1 + d;
    }
    
    printf("%lld\n",ans1);
}

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值