运动员最佳配对问题
西安交大 軟件53 蔡少斐
題號:6_5
題目敘述:
羽毛球隊有男女運動員各n人。
給定2個n×n矩陣P和Q。P[i][j]是男運動員i和女運動員j配對組成混合雙打的男運動員競賽優勢;Q[i][j]是女運動員i和男運動員j配合的女運動員競賽優勢。
?
由于技術配合和心理狀態等各種因素影響,P[i][j]不一定等于Q[j][i]。男運動員i和女運動員j配對組成混合雙打的男女雙方競賽優勢為P[i][j]*Q[j][i]。
?
設計一個算法,計算男女運動員最佳配對法,使各組男女雙方競賽優勢的總和達到最大。
?
編程任務:設計一個算法,對于給定的男女運動員競賽優勢,計算男女運動員最佳配對法,使各組男女雙方競賽優勢的總和達到最大。
輸入格式
輸入數據第一行有1 個正整數n (1≤n≤10)。接下來的2n行,每行n個數。前n行是P,后n行是Q。
?
輸出格式
將計算出的男女雙方競賽優勢的總和的最大值輸出。
輸入樣例
3?
10 2 3?
2 3 4?
3 4 5?
2 2 2?
3 5 3?
4 5 1?
輸出樣例
52?
題目解答:
?????? 這道題目共有n!種配對情況,也就是相當于固定男運動員,然后對女運動員進行一次全排列,并求出對應的優勢之和的最大值,本題可以用回溯法,也可以用分支限界法,在使用分支限界法的時候,關鍵是在于設計上界函數。
?????? 在這里,我們把上界函數定義為:剩下的未配對的女運動員(不考慮男運動員配對情況下)所能達到的優勢最大值之和(記為r)與當前配對已達到的優勢(記為sum)之和。在程序里體現如下:
在一開始時候,r被初始化為:
其中maxout的定義為:
?????? 在使用分支限界法的時候,一旦有一個葉節點出來,那么就立即結束算法,因為最先出來的葉節點必定是最優解。
代碼實現:
#include<algorithm> #include<iostream> #include<queue> #include<vector> #include <cstdio> usingnamespace std; const int MAX = 20; int P[MAX][MAX]; int Q[MAX][MAX]; int maxout[MAX]; int n; structnode {int id;int sum;int r;int up;int *x; }; structcmp {bool operator()( node* a, node* b ){return(a->up < b->up);} }; voidsolve() {int ans= 0;priority_queue<node*,vector<node*>, cmp> que;node *E = new node();E->id =1;E->sum =0;E->r =0;E->up =0;for ( int i = 1; i <= n; i++ ){E->r += maxout[i];}E->up =E->r;E->x =new int[n + 1];for ( int i = 1; i <= n; i++ ){E->x[i] = i;}while ( E->id != n + 1 ){for ( int i = E->id; i <= n;i++ ){node* nE = new node();nE->id = E->id + 1;nE->x = new int[n + 1];for ( int t = 1; t <= n;t++ ){nE->x[t] =E->x[t];}nE->x[E->id] = E->x[i];nE->x[i] = E->x[E->id];nE->sum = E->sum +P[E->id][nE->x[E->id]] * Q[nE->x[E->id]][E->id];nE->r = E->r - maxout[E->id];nE->up = nE->sum + nE->r;que.push( nE );}if ( !que.empty() ){E= que.top();que.pop();}else {ans = 0;break;}}ans = E->sum;cout << ans << endl; }intmain() {cin >> n;for ( int i = 1; i <= n; i++ ){for ( int j = 1; j <= n; j++ ){cin >> P[i][j];}}for ( int i = 1; i <= n; i++ ){for ( int j = 1; j <= n; j++ ){cin >> Q[i][j];}}for ( int i = 1; i <= n; i++ ){int ma = 0;for ( int j = 1; j <= n; j++ ){ma = max( ma, Q[i][j] *P[j][i] );}maxout[i] = ma;}solve();return(0); }/** 3* 10 2 3* 2 3 4* 3 4 5* 2 2 2* 3 5 3* 4 5 1*/運行結果:
總結
- 上一篇: 装习惯了ATX,想来试试ITX?老司机都
- 下一篇: 动态规划练习1 [导弹拦截]