1. 程式人生 > >hdu 2988 Strange fuction【模擬退火】

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;
}