【UVA】10003-Cutting Sticks（动态规划、矩阵链乘）

一道动态规划题，不过似乎可以用回溯水过去，回溯的话效率很烂的。

13988658

10003

Cutting Sticks

Accepted

C++

1.882

2014-08-04 09:26:49

AC代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<list>
#include<string>
#include<sstream>
#include<ctime>
using namespace std;
#define _PI acos(-1.0)
#define INF 1 << 10
#define esp 1e-6
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> pill;
/*===========================================*/
#define MAXD 1000 + 10
int n,L;
int arr[MAXD];
int vis[MAXD];
int _dp[MAXD][MAXD];
int dp(int x,int y,int have_cut){ /*当前木棍区间，以及砍了几段木头了*/
    int ans = INF;
    int ok  = 0;
    if(_dp[x][y] > 0) return _dp[x][y];
    for(int i = x ; i <= y ; i++)
        for(int k = 0 ; k < n ; k++)if(!vis[k] && arr[k] >= x && arr[k] <= y){ /*该点未切割*/
            int pos = arr[k];
            ok = 1;
            vis[k] = 1;
            ans = min(ans,dp(x,pos,have_cut + 1) + dp(pos,y,have_cut + 1));
            vis[k] = 0;
    }
    if(!ok)
        return 0;
    else
        return _dp[x][y] = ans + (y - x);
}
int main(){
    while(scanf("%d",&L) && L){
        memset(vis,0,sizeof(vis));
        memset(_dp,-1,sizeof(_dp));
        scanf("%d",&n);
        for(int i = 0 ; i < n ; i++)
            scanf("%d",&arr[i]);
        int ans = dp(0,L,0);
        printf("The minimum cutting is %d.\n",ans);
    }
    return 0;
}

之后是动规递推

13989008

10003

Cutting Sticks

Accepted

C++

0.102

2014-08-04 10:56:30

可以看出来，这效率差的真不是一点半点的。

利用dp[i][j]表示切点i,j之间的线段分割的话需要产生的最小花费，那么可以有递推式：

dp[i][j] = dp[i][k] + dp[k][j] + array[j] - array[i]

所以如果 i = j - 1的话,dp[i][j] = 0;

递推的时候，区间长度大的可以由小的推导得到，所以按照区间长度的大小进行递推

AC代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<list>
#include<string>
#include<sstream>
#include<ctime>
using namespace std;
#define _PI acos(-1.0)
#define INF 1 << 20
#define esp 1e-6
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> pill;
/*===========================================*/
#define MAXD 50 + 5
int L,n,array[MAXD];
int dp[MAXD][MAXD];
bool Input(){
    scanf("%d",&L);
    if(!L)
        return false;
    scanf("%d",&n);
    for(int i = 1 ; i <= n ;i++)
        scanf("%d",&array[i]);
    array[0] = 0;
    array[n + 1] = L;
    return true;
}
void init(){
    for(int i = 0 ; i < n + 1; i++)
        for(int j = i + 1; j <= n + 1; j++)
            if(i + 1 == j)
                dp[i][j] = 0;
            else
                dp[i][j] = INF;
    return ;
}
void DP(){
    init();
    /*按照区间长度由小到大的长度进行递推*/
    for(int j = 2 ; j <= n + 1; j++)
        for(int i = 0; i + j <= n + 1 ; i++){
            for(int k = i + 1; k < i + j; k++)
                dp[i][i + j] = min(dp[i][i + j],dp[i][k] + dp[k][i + j] + array[i + j] - array[i]);
    }
    printf("The minimum cutting is %d.\n",dp[0][n + 1]);
}
int main(){
    while(Input()){
        DP();
    }
    return 0;
}

一下是一个矩阵链乘的链接，感觉还不错：

http://blog.csdn.net/liufeng_king/article/details/8497607