hdu 2988 Strange fuction【模擬退火】
阿新 • • 發佈:2018-11-30
計算:給出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; }