1. 程式人生 > >UVA - 348Optimal Array Multiplication Sequence(遞推)

UVA - 348Optimal Array Multiplication Sequence(遞推)

type tip track 而是 popu ret 滿足 -m lan

題目:Optimal Array Multiplication Sequence


題目大意:給出N個矩陣相乘。求這些矩陣相乘乘法次數最少的順序。


解題思路:矩陣相乘不滿足交換率但滿足結合率。dp【i】【j】 代表第1個矩陣到第j個矩陣之間的最少的乘法次數,轉移狀態方程:dp【i】【j】 = Min(dp【i】【k】 + dp【k + 1】【j】 + A[i - 1] * A[k] *A[j]) k>= i && k <= j - 1。A0A1A2A3..Ak | Ak+1AK+2Aj, 第i個矩陣的行Ai - 1 ,列Ai。上面那個式子的意思是在K這個地方斷開。兩邊的矩陣都先做相乘,這樣左右兩邊的相乘是互不影響的。這樣左右兩邊的相乘也是採取找最優的形式(最優的子結構)。全部就有上面的狀態轉移方程。最後是路徑輸出,我沒有記錄路徑而是直接用遞歸的寫法,可是這裏要註意僅僅要找一種最優序列就好了(有的矩陣相乘存在多種最優序列),在遞歸中少了一個break害我wa了好久,找不到原因。


代碼:

#include <cstdio>
#include <cstring>

typedef long long ll;

const int INF = 0x7fffffff;
const int N = 15;

ll A[N];
int n;
ll dp[N][N];

ll Min (const ll a, const ll b) { return a < b?

a: b;} void printf_ans (int s, int e) { if (s == e) { printf ("A%d", s); return; } for (int i = s; i < e; i++) { if ((dp[s][i] + dp[i + 1][e] + A[s - 1] * A[i] * A[e]) == dp[s][e]) { if (s != i) printf ("("); printf_ans(s, i); if (s != i) printf (")"); printf (" x "); if (i + 1 != e) printf ("("); printf_ans(i + 1, e); if (i + 1 != e) printf (")"); break; } } } int main () { int cas = 0; while (scanf ("%d", &n), n) { for (int i = 0; i < n; i++) scanf ("%lld%lld", &A[i], &A[i + 1]); for (int i = 1; i <= n; i++) dp[i][i] = 0; ll temp; for (int len = 1; len < n; len++) for (int i = 1; i + len <= n; i++) { temp = INF; for (int k = 0; k < len; k++) temp = Min (temp, dp[i][i + k] + dp[i + k + 1][i + len] + A[i - 1] * A[i + k] * A[i + len]); dp[i][i + len] = temp; } // printf ("%lld\n", dp[1][n]); printf ("Case %d: (", ++cas); printf_ans(1, n); printf (")\n"); } return 0; }



UVA - 348Optimal Array Multiplication Sequence(遞推)