Advertisement
Morass

11158

Jan 31st, 2017
167
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.75 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define PB push_back
  4. #define ZERO (1e-10)
  5. #define INF (1<<29)
  6. #define CL(A,I) (memset(A,I,sizeof(A)))
  7. #define DEB printf("DEB!\n");
  8. #define D(X) cout<<"  "<<#X": "<<X<<endl;
  9. #define EQ(A,B) (A+ZERO>B&&A-ZERO<B)
  10. typedef long long ll;
  11. typedef long double ld;
  12. typedef pair<ll,ll> pll;
  13. typedef vector<int> vi;
  14. typedef pair<int,int> ii;
  15. typedef vector<ii> vii;
  16. #define IN(n) int n;scanf("%d",&n);
  17. #define FOR(i, m, n) for (int i(m); i < n; i++)
  18. #define REP(i, n) FOR(i, 0, n)
  19. #define F(n) REP(i, n)
  20. #define FF(n) REP(j, n)
  21. #define FT(m, n) FOR(k, m, n)
  22. #define aa first
  23. #define bb second
  24. void ga(int N,int *A){F(N)scanf("%d",A+i);}
  25. int brt(int*A,int N){
  26.     sort(A,A+N);
  27.     int S=0,B=0;
  28.     do{
  29.         S=0;
  30.         F(N-1)S+=abs(A[i]-A[i+1]);
  31.         B=max(B,S);
  32.     }while(next_permutation(A,A+N));
  33.     return B;
  34. }
  35. #define MX (55)
  36. #define DV (10000)
  37. #define TEMP 1e5
  38. #define CRT .99995
  39. #define PRC (4)
  40. int nxt(int*A,int N,int&I){
  41.     I=rand()%(N-1);
  42.     int S=0;
  43.     if(I)S-=abs(A[I]-A[I-1]);
  44.     if(I<N-2)S-=abs(A[I+1]-A[I+2]);
  45.     swap(A[I],A[I+1]);
  46.     if(I)S+=abs(A[I]-A[I-1]);
  47.     if(I<N-2)S+=abs(A[I+1]-A[I+2]);
  48.     return S;
  49. }
  50. int solANL(int*A,int N){
  51.     if(N<11)return brt(A,N);
  52.     double T=TEMP,d,c=CRT,a=1e-5;
  53.     int S=0,C=0,B,I;
  54.     random_shuffle(A,A+N);
  55.     F(N-1)C-=abs(A[i]-A[i+1]);
  56.     while(T>a){
  57.         B=C-nxt(A,N,I);    
  58.         d=B-C;
  59.         if(d<0||exp(-d/T)>(rand()%DV*1./DV))
  60.             S=min(S,C=B);
  61.         else swap(A[I],A[I+1]);
  62.         T*=c;
  63.     }
  64.     return -S;
  65. }
  66. int main(void){
  67.     int N,A[MX];
  68.     srand(time(0));
  69.     IN(_)F(_){
  70.         scanf("%d",&N),ga(N,A);
  71.         printf("Case %d: %d\n",i+1,solANL(A,N));
  72.     }
  73.     return 0;
  74. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement