[HDU 5107][2014 ACM/ICPC Asia Regional Xi'an Online]Ellipsoid(爬山)

55 篇文章 0 订阅
15 篇文章 0 订阅

题目链接

http://acm.hdu.edu.cn/showproblem.php?pid=5017

思路

裸的爬山啊!!网上一堆题解,一边代码里写的明明是爬山,一边标题上白纸黑字说是退火的是什么心态!代码还都差不多,我了个去
弄得我很naive地去写了个真退火,秒WA

代码

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

#define INF 1e12
#define EPS 1e-8

using namespace std;

double a,b,c,d,e,f;

double dist(double x,double y,double z) //求(x,y,z)到原点距离
{
    return sqrt(x*x+y*y+z*z);
}

double calc(double x,double y) //已知x,y坐标,求z坐标
{
    double delta=(d*y+e*x)*(d*y+e*x)-4*c*(a*x*x+b*y*y+f*x*y-1);
    if(delta<0) //delta<0则方程无解
        return INF;
    double z1=(-d*y-e*x+sqrt(delta))/(2*c),z2=(-d*y-e*x-sqrt(delta))/(2*c);
    if(dist(x,y,z1)<dist(x,y,z2)) return z1; //返回两个z中更优的那个
    return z2;
}

int xx[]={0,0,1,-1,1,-1,1,-1},yy[]={1,-1,1,-1,-1,1,0,0};

double SA()
{
    double T=1,x=0,y=0,z=sqrt(1/c); //当前最优解为(x,y,z)
    double minans=dist(x,y,z);
    while(T>EPS)
    {
        for(int dir=0;dir<8;dir++)
        {
            double newx=x+xx[dir]*T,newy=y+yy[dir]*T;
            double newz=calc(newx,newy);
            double nowans=dist(newx,newy,newz);
            if(nowans<minans) //新的解更优,移动到新解上
            {
                minans=nowans;
                x=newx,y=newy,z=newz;
            }
        }
        T*=0.99;
    }
    return minans;
}

int main()
{
    while(scanf("%lf%lf%lf%lf%lf%lf",&a,&b,&c,&d,&e,&f)!=EOF)
        printf("%.7lf\n",SA());
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值