SGU 485 Arrays

485. Arrays

Time limit per test: 1.75 second(s)
Memory limit: 262144 kilobytes

input: standard
output: standard

You are given a sequence of 3· N integers (X₁, X₂, ·s, X_3· N). Create three sequences (A₁, A₂, ·s, A_N), (B₁, B₂, ·s, B_N) and (C₁, C₂, ·s, C_N) such that:

• each of the integers from 1 to 3· N belongs to exactly one of the sequences A, B or C;
• the value of  is the largest possible.

Input

Constraints on N	Constraints on T
1 ≤ N ≤ 10	1 ≤ T ≤ 1000
11 ≤ N ≤ 15	1 ≤ T ≤ 100
16 ≤ N ≤ 20	1 ≤ T ≤ 10
21 ≤ N ≤ 25	T = 1

The input file contains T test cases, all having the same value of N. The first line of the input file contains the integers T and N, constrained as shown in the adjacent table. Each of the following T lines describes one test case and contains 3· N integers, the members of the sequence X. All these values are in the range from 0 to 1000.

Output

The output file should consist of T lines. Each line should contain the largest possible value of S for the corresponding test case from the input file.

Example(s)

sample input	sample output
1 2 4 1 8 2 0 5	46

Note. The maximal value is attained by taking A = (1, 3), B = (2, 5), C = (4, 6).

题意：给出一组数，将这些数分成三组，记为A,B,C，求出满足sigma[(Ai-Bi)*Ci]的最大值。

sl: 首先考虑下B ,很显然B中的元素应该是最小的N个，在考虑A,C很容易看出A,C应该是满足

A>=C的关系。

所以我们只需要枚举A,C就好了。

其中包含2个优化

优化1：在计算的过程中应满足全局平均值最小。

优化2：不等式 （ai-bi）*ci+(aj-bj)*cj - {(ai-bi)*aj+(ci-bj)*cj   } >0成立 即满足：

（ci-aj）*(ai-bi-cj)>0 成立。

所以ci>aj 时 ai-bi>cj

ci<aj时 ai-bi<cj

详见代码。

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <algorithm>
 4 using namespace std;
 5 typedef long long LL;
 6 const int maxn = 76;
 7 int T, N;
 8 int ans, a[maxn], va[maxn], vc[maxn];
 9 bool vis[maxn];
 bool test(int cur) {
     for (int i=1; i<cur; ++i) {
         if (abs(va[cur])>abs(vc[i]) && abs(vc[cur])<abs(va[i])-abs(a[3*N-i+1])) {
             return false;
         }
         if (abs(va[cur])<abs(vc[i]) && abs(vc[cur])>abs(va[i])-abs(a[3*N-i+1])) {
             return false;
         }
     }
     return true;
 }
 void dfs(int cur, int last, int val) {
     if (cur==N+1) {
         ans = max(ans, val);
         return;
     }
     for (int i=cur; i<=2*N; ++i) {
         if (!vis[i]) {
             vis[i] = 1;
             va[cur] = a[i];
             for (int j=max(last+1, i+1); j<=N*2; ++j) {
                 if (!vis[j]) {
                     vc[cur] = a[j];
                     vis[j] = 1;
                     int netVal = val+(a[i]-a[3*N-cur+1])*a[j];
                     if (netVal*N>ans*cur) {
                         if (test(cur)) {
                             dfs(cur+1, j, netVal);
                         }
                     }
                     vis[j] = 0;
                 }
             }
             vis[i] = 0;
             break;
         }
     }
 }
 int main()
 {
     scanf ("%d%d", &T, &N);
     while (T--) {
         for (int i=1; i<=3*N; ++i) {
             scanf ("%d", &a[i]);
             a[i] = -a[i];
         }
         sort(a+1, a+3*N+1);
         ans = 0;
         dfs(1, 0, 0);
         printf ("%d\n", ans);
     }
     return 0;

62 }