HDU5074

Hatsune Miku

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others)
Total Submission(s): 619    Accepted Submission(s): 445

Problem Description

Hatsune Miku is a popular virtual singer. It is very popular in both Japan and China. Basically it is a computer software that allows you to compose a song on your own using the vocal package.

Today you want to compose a song, which is just a sequence of notes. There are only m different notes provided in the package. And you want to make a song with n notes.

Also, you know that there is a system to evaluate the beautifulness of a song. For each two consecutive notes a and b, if b comes after a, then the beautifulness for these two notes is evaluated as score(a, b).

So the total beautifulness for a song consisting of notes a1, a2, . . . , an, is simply the sum of score(ai, ai+1) for 1 ≤ i ≤ n - 1.

Now, you find that at some positions, the notes have to be some specific ones, but at other positions you can decide what notes to use. You want to maximize your song’s beautifulness. What is the maximum beautifulness you can achieve?

 

 

Input

The first line contains an integer T (T ≤ 10), denoting the number of the test cases.

For each test case, the first line contains two integers n(1 ≤ n ≤ 100) and m(1 ≤ m ≤ 50) as mentioned above. Then m lines follow, each of them consisting of m space-separated integers, the j-th integer in the i-th line for score(i, j)( 0 ≤ score(i, j) ≤ 100). The next line contains n integers, a1, a2, . . . , an (-1 ≤ ai ≤ m, ai ≠ 0), where positive integers stand for the notes you cannot change, while negative integers are what you can replace with arbitrary notes. The notes are named from 1 to m.

 

 

Output

For each test case, output the answer in one line.

 

 

Sample Input

25 383 86 7715 93 3586 92 493 3 3 1 210 536 11 68 67 2982 30 62 23 6735 29 2 22 5869 67 93 56 1142 29 73 21 19-1 -1 5 -1 4 -1 -1 -1 4 -1

 

 

Sample Output

270625

解题思路：

其实听学长讲完就基本理解了，暴力DP。可能网上有更优的代码，你们开心就好。

题意为有m种音符，编号1到m，我们要用这m种音符来创造一首带有n个音符的曲子（当然，一种音符可以用多次）a[n],假设有两个连续的音符 i ,j ，那么定义score(i,j)为这两个音符的得分，题目中预先给出所有的score(i,j)   1<=i,j<=m, 那么我们创造出的n个音符的曲子的得分为  score( a[i] , a[i+1] ) + score (a[i+1] ,a[i+2) +......score(a[n-1],a[n]) ， i从1开始。 a [i] 代表n个音符的曲子中第i个音符是第几种音符， 1<=a[i]<=m， 比如 a[i]=3，就表示n个音符中第i个位置用的是第3种音符（一共有m种）,预先给出 这n个音符，a[1] 到a[n] ，其中a[ i] 或者等于-1 ，或者 大于等于1小于等于m，对于后者，该位置的音符不能改变，对于前者，该位置可以换成任意的音符j(1<=j<=m)， 问 这n个音符所获得的最大得分是多少。

第一个例子：

5 3 

83 86 77 

15 93 35 

86 92 49 

3 3 3 1 2

3*3的矩阵，5个数3 3 3 1 2的意思是依次执行4个操作

Score[3][3]+score[3][3]+score[3][1]+score[1][2];

49+49+86+86=270;

 

第二个例子：

10 536 11 68 67 2982 30 62 23 6735 29 2 22 5869 67 93 56 1142 29 73 21 19-1 -1 5 -1 4 -1 -1 -1 4 -1

依次执行 下面为对应的dp[i][j]的值

0    0    0    0   0

82   67   93   67  67

0    0    0    0   151

193  180  224  172  170

0    0    0   260   0

329  327  353  316  271

409  383  409  396  411

465  463  489  476  467

0    0    0   532  0

601  599  625  588  543

结果显然为625

代码思路：

定义dp [i ][j]，为前i个音符且第i个音符为第j种音符所获得的最大得分。

分情况讨论：

当a[i] >0 （i>=2)时：

当a[i-1]>0时： 

dp[i][a[i]]=dp[i-1][a[i-1]]+s[a[i-1]][a[i]];

当a[i-1]<0时： 

dp[i][a[i]]=max(dp[i][a[i]],dp[i-1][j]+s[j][a[i]]);

   1<=j<=m,枚举第i-1个音符是第几种音符

当a[i]<0 (i>=2)时：

当a[i-1]>0时：

dp[i][j]=max(dp[i][j],dp[i-1][a[i-1]]+s[a[i-1]][j]);

 1<=j<=m,枚举第i个音符是第几种音符

当note[i-1]<0时：

dp[i][j]=max(dp[i][j],dp[i-1][k]+s[k][j]); 

 1<=j,k<=m, 枚举第i个音符和第i-1个音符是第几种音符

解题代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<vector>
#include<queue>
#include<set>
#include<cmath>
#define mst(ss,b) memset((ss),(b),sizeof(ss))
#define maxn 0x3f3f3f3f
using namespace std;
typedef __int64 ll;
typedef long L;
int dp[105][55];
int s[55][55];
int a[105];
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        int n,m;
        mst(dp,0);
        cin>>n>>m;
        for(int i=1; i<=m; i++)
        {
            for(int j=1; j<=m; j++)
            {
                cin>>s[i][j];
            }
        }
        for(int i=1; i<=n; i++)
        {
            cin>>a[i];
        }
        for(int i=2; i<=n; i++)
        {
            if(a[i]!=-1)
            {
                if(a[i-1]!=-1)
                {
                    dp[i][a[i]]=dp[i-1][a[i-1]]+s[a[i-1]][a[i]];
                }
                else
                {
                    for(int j=1; j<=m; j++)
                    {
                        dp[i][a[i]]=max(dp[i][a[i]],dp[i-1][j]+s[j][a[i]]);
                    }
                }
            }
            else
            {
                if(a[i-1]!=-1)
                {
                    for(int j=1; j<=m; j++)
                    {
                        dp[i][j]=max(dp[i][j],dp[i-1][a[i-1]]+s[a[i-1]][j]);
                    }
                }
                else
                {
                    for(int j=1; j<=m; j++)
                    {
                        for(int k=1; k<=m; k++)
                        {
                            dp[i][j]=max(dp[i][j],dp[i-1][k]+s[k][j]);
                        }
                    }
                }
            }
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {
                cout<<dp[i][j]<<" ";
            }
            cout<<endl;
        }
        int ans=-1;
        for(int i=1; i<=m; i++)
        {
            if(ans<dp[n][i])
                ans=dp[n][i];
        }
        //cout<<ans<<endl;
    }
    return 0;
}