luogu P2210 Haywire

版权声明:thanks for your copy! https://blog.csdn.net/weixin_43907802/article/details/88926984

在这里插入图片描述

analysis

模拟退火大法好!!!
让我们一起念:

P ( d E ) = d e t l a t e m p e r a t u r e P(dE)=\frac{detla}{temperature}

好吧,思路很简单,每次交换两头奶牛的位置,然后评价,决定是否保留,若不保留就换回去,然后又换,又评价…
这就是SA的步骤了:

  • 得到一个状态s1
  • 随便你怎么转移到另一个状态s2
  • 评价状态s1是否可以向更优转移,若不能,用上面的概率公式来判断是否暂且接受较差解s2,若能,更新s1更优状态s2
  • 循环直到温度降低到精度范围以内结束

code

#include<bits/stdc++.h>
using namespace std;
#define loop(i,start,end) for(register int i=start;i<=end;++i)
#define anti_loop(i,start,end) for(register int i=start;i>=end;--i)
#define clean(arry,num); memset(arry,num,sizeof(arry));
#define max(a,b) ((a>b)?a:b)
#define min(a,b) ((a<b)?a:b)
#define asc(a) (a-'A')
int n;
const int maxn=15;
int pos[maxn];
int f[maxn][3];
bool con[maxn][maxn];
template<typename T>void read(T &x)
{
    x=0;char r=getchar();T neg=1;
    while(r>'9'||r<'0'){if(r=='-')neg=-1;r=getchar();}
    while(r>='0'&&r<='9'){x=(x<<3)+(x<<1)+r-'0';r=getchar();}
    x*=neg;
}
inline int abs(int x){return ((x>=0)?x:-x);}
inline int judge()
{
    clean(con,false);
    int res=0;
    loop(i,1,n)
    {
        loop(j,0,2)
        {
            if(con[f[i][j]][i])continue;
            res+=abs(pos[i]-pos[f[i][j]]);
            con[f[i][j]][i]=true;
        }
    }
    return res;
}
inline void SA()
{
    int res=judge();
    srand(19260817);//东方八位神秘质数
    double Temp=1e10,eps=1e-17,detla=0.99;
    for(;Temp>=eps;Temp*=detla)
    {
        int swap1=rand()%n+1;
        int swap2=rand()%n+1;
        if(swap1==swap2)continue;
        swap(pos[swap1],pos[swap2]);
        int output=judge();
        if(output<res)
        {
            res=output;
        }
        else if(res<=output&&exp((res-output)/Temp)>=((double)rand()/(double)RAND_MAX))
        {
            res=output;
        }
        else swap(pos[swap1],pos[swap2]);
    }
    printf("%d",res>>1);
}
int main()
{
    #ifndef ONLINE_JUDGE
    freopen("datain.txt","r",stdin);
    #endif
    read(n);
    loop(i,1,n)pos[i]=i,read(f[i][0]),read(f[i][1]),read(f[i][2]);
    SA();
    return 0;
}

猜你喜欢

转载自blog.csdn.net/weixin_43907802/article/details/88926984