本文主要是介绍hdu 2988 Strange fuction【模拟退火】,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
计算:给出y , 的最小值
wa到哭啊,简直上火了,是板子没有套对, nex=now+T*i;//新解 T*i是变化的范围加上原先的才是新解,
还有精确度,因为评估函数要平方,所以要高一些 最终15ms
#include<bits/stdc++.h>
#define mp make_pair
#define sz(x) int((x).size())
#define fin freopen("in.txt","r",stdin)
#define fout freopen("out.txt","w",stdout)
#define io ios::sync_with_stdio(0),cin.tie(0)
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int inf = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const int maxn = 1e4 + 5;
double n,x,y;
double now;
double f(double x)
{ //评估函数return 6.0*pow(x,7)+8.0*pow(x,6)+7.0*pow(x,3)+5.0*pow(x,2)-y*x;
}
double ans;
double sa()
{double ansp,nex;ans=1e111;//最开始的能量值,初始很大就可以,不用修改double T=100.0; //初始温度, (可以适当修改,最好和给的数据最大范围相同,或者缩小其原来0.1)double d=0.98; //降温系数 (可以适当修改,影响结果的精度和循环的次数,)double eps=1e-8; //最终温度 (要是因为精度问题,可以适当减小最终温度)double TT=0.99; //采纳新解的初始概率double dd=0.99; //(可以适当修改,采纳新解变更的概率)(这个概率下面新解更新的时候,最好和未采纳的新解更新的次数是一半一半)double res=f(now); //传入的初始默认解(now)下得到的评估能量值if (res<ans) ans=res;int num=1;while (T>eps){for (int i=-1;i<=1;++i)if (now+i*T<=100&&now+T*i>=0){nex=now+T*i;//新解double tmp=f(nex);//新解下的评估能量值if (tmp<ans) ans=tmp;//降温成功,更新当前最优解if (tmp<res) res=tmp,now=nex;// 降温成功,采纳新解else if (TT>rand()%10000/10000.0) res=tmp,now=nex;//,cout<<"======"<<endl;//没有 降温成功,但是以一定的概率采纳新解//else cout<<"="<<endl;//用于测试,设定的采纳新解的概率,是否为一半一半,可以适当修改降温参数dd}T*=d; TT*=dd;//num++;}//cout<<num<<endl;return ans;
}
int main()
{srand(time(0));scanf("%lf",&n);while(n--){cin>>y;now=50.0;sa();printf("%.4f\n",ans);}return 0;
}
这篇关于hdu 2988 Strange fuction【模拟退火】的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!