POJ 1018 Communication System

这道题跟那天做到的jury很像，第一感觉是开dp，dp[i][j]，j参数记录题目要求的B P的其中一个，然后dp[i][j]的值来表示另外一个。当时觉得最后算max(Bmin/Psum)的时候还需要再开一个path[][]记录路径。开始写之后发现其实用不到。
dp的核心还是在于当前状态与“上一”状态的关系，换句话说真正用得到的就只有两个状态。这里虽然看起来每次都要比较min{b1,b2,b3,b4…}，但是如果我们用dp[i][j]里的j来表示minB的话，其实每一次只要与dp[i-1]时的j比较就可以，也正因此不需要path[][]来回溯之前所有选过的b求最小值。如果用j来表示Psum的话是一样的，也不用回溯一次次地相加。
最开始我是用j来表示Psum的值，但是因为题目只说positive integer而没有更明确的限制，事实上如果我们来表示所有可能的p和的数值其实是不行的（后来也证明不行）样例的数据范围都很小，所以能过。
后来去搜题解，题解说可以一点点试试出一个大概的数据范围（学到了），试出来b的话是[1,400]。所以如果dp[i][j]表示选取到第i个物品且其Bmin为j时Psum的值，就可以了。
不过用j来表示Psum和Bmin其实核心的思想都是一样的，但是要注意表示Bmin时，dp[i][j]的值指的是当前的Psum，因为我们要求（Bmin/Psum）的最大值所以Psum要取min。这里最开始我想错了，一直取的max。
还有一个点就是最后循环求ansmax的时候，注意分数不要写反分子分母。

这是用j表示Psum的dp代码。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <iomanip>
using namespace std;

int dp[110][10000], path[110][10000], b[110][110], p[110][110], num[110];
int t;
//if b,p [1, 100]
int main()
{
    cin >> t;
    while(t)
    {
        t--;
        int n;
        cin >> n;
        for(int i = 1;i <= n;i++)
        {
            cin >> num[i];
            for(int j = 1;j <= num[i];j++)
                cin >> b[i][j] >> p[i][j];
        }

        //cout << "input fin" << endl;

        memset(dp, -1, sizeof(dp));
        memset(path, 0, sizeof(path));

        //dp[i][j]挑选至第i个物品，p和为j时，最小的b
        dp[0][0] = 0;

        for(int i = 1;i <= n;i++)
        {
            //cout << "into loop" << endl;
            for(int j = 0;j <= 10000;j++)
            {
                //cout << "j loop" << endl;
                if(dp[i-1][j] >= 0)//前一状态存在
                {
                    //cout << "meet need" << endl;

                    for(int k = 1;k <= num[i];k++)
                    {
                        //cout << "k loop" << endl;
                        if(dp[i-1][j] < b[i][k] && dp[i-1][j] > 0)
                            dp[i][j+p[i][k]] = dp[i-1][j];

                        else
                            dp[i][j+p[i][k]] = b[i][k];

                        //printf("dp[%d][%d] = %d\n", i, j+p[i][k], dp[i][j+p[i][k]]);

                    }

                }
            }
        }

        //cout << "ans compare" << endl;
        double ans = 0;
        for(int i = 1;i <= 10000;i++)
        {
            if(dp[n][i] > 0)
                ans = max(double(dp[n][i]) / double(i), ans);
        }

        cout << fixed << setprecision(3) << ans << endl;
    }

    return 0;
}

这是用j表示Bmin的AC代码。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <iomanip>
using namespace std;

int dp[110][10000], path[110][10000], b[110][110], p[110][110], num[110];
int t;
//if b,p [1, 100]
int main()
{
    cin >> t;
    while(t)
    {
        t--;
        int n;
        cin >> n;
        for(int i = 1;i <= n;i++)
        {
            cin >> num[i];
            for(int j = 1;j <= num[i];j++)
                cin >> b[i][j] >> p[i][j];
        }

        //cout << "input fin" << endl;

        memset(dp, -1, sizeof(dp));
        memset(path, 0, sizeof(path));

        //dp[i][j]挑选至第i个物品,bmin = j时,psum的值
        dp[0][0] = 0;

        for(int i = 1;i <= n;i++)
        {
            //cout << "into loop" << endl;
            for(int j = 0;j <= 10000;j++)
            {
                //cout << "j loop" << endl;
                if(dp[i-1][j] >= 0)//前一状态存在
                {
                    //cout << "meet need" << endl;

                    for(int k = 1;k <= num[i];k++)
                    {
                        //cout << "k loop" << endl;
                        if(b[i][k] <= j || j == 0)
                        {
                            if(dp[i][b[i][k]] >= 0)
                                dp[i][b[i][k]] = min(dp[i][b[i][k]], dp[i-1][j]+p[i][k]);
                            else
                                dp[i][b[i][k]] = dp[i-1][j]+p[i][k];
                            //printf("con1  dp[%d][%d] = %d\n", i, b[i][k], dp[i][b[i][k]]);
                        }

                        else
                        {
                            if(dp[i][j] >= 0)
                                dp[i][j] = min(dp[i][j], dp[i-1][j]+p[i][k]);
                            else
                                dp[i][j] = dp[i-1][j]+p[i][k];
                            //printf("con2  dp[%d][%d] = %d\n", i, j, dp[i][j]);
                        }
                    }
                }
            }
        }

        //cout << "ans compare" << endl;
        double ans = 0;
        for(int i = 1;i <= 10000;i++)
        {
            if(dp[n][i] > 0)//dp[n][i] bmin=i时 sump
                ans = max(double(i) / double(dp[n][i]), ans);
        }

        cout << fixed << setprecision(3) << ans << endl;
    }

    return 0;
}