POJ 2609 Ferry Loading

双塔DP+输出路径。

由于内存限制，DP只能开滚动数组来记录。

我的写法比较渣，但是POJ能AC，但是ZOJ依旧MLE，更加奇怪的是Uva上无论怎么改都是WA，其他人POJ过的交到Uva也是WA。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<stack>
#include<algorithm>
using namespace std;

int dp[][+];
struct P
{
    int x,y;
    int f;
}pre[][+];
int sum;
int a[],tot;

void read()
{
    while()
    {
        scanf("%d",&a[tot+]);
        if(a[tot+]==) break;
        tot++;
    }

    for(int i=;i<=tot;i++)
        for(int j=;j<=*sum;j++)
            pre[i][j].f=-;

    for(int j=;j<=*sum;j++) dp[][j]=dp[][j]=0x7FFFFFFF;
}

void init()
{
    sum=sum*;
    tot=;
}

void work()
{
    int now=;
    dp[][a[]+sum]=;
    pre[][a[]+sum].f=;
    dp[][-a[]+sum]=;
    pre[][-a[]+sum].f=;

    for(int i=;i<=tot;i++)
    {
        now=now^;
        for(int j=;j<=*sum;j++) dp[now][j]=0x7FFFFFFF;
        for(int j=;j<=*sum;j++)
        {
            if(dp[now^][j]==0x7FFFFFFF) continue;

            int tmp=j-sum;

            if(tmp>=)
            {
                if(dp[now^][j]+tmp+a[i]<=sum&&dp[now^][j]<dp[now][tmp+a[i]+sum])
                {
                    dp[now][tmp+a[i]+sum]=dp[now^][j];
                    pre[i][tmp+a[i]+sum].x=i-;
                    pre[i][tmp+a[i]+sum].y=j;
                    pre[i][tmp+a[i]+sum].f=;
                }

                if(max(tmp,a[i])+dp[now^][j]<=sum&&dp[now^][j]+min(tmp,a[i])<dp[now][tmp-a[i]+sum])
                {
                    dp[now][tmp-a[i]+sum]=dp[now^][j]+min(tmp,a[i]);
                    pre[i][tmp-a[i]+sum].x=i-;
                    pre[i][tmp-a[i]+sum].y=j;
                    pre[i][tmp-a[i]+sum].f=;
                }
            }

            else
            {
                tmp=-tmp;
                if(dp[now^][j]+tmp+a[i]<=sum&&dp[now^][j]<dp[now][-(tmp+a[i])+sum])
                {
                    dp[now][-(tmp+a[i])+sum]=dp[now^][j];
                    pre[i][-(tmp+a[i])+sum].x=i-;
                    pre[i][-(tmp+a[i])+sum].y=j;
                    pre[i][-(tmp+a[i])+sum].f=;
                }

                if(max(tmp,a[i])+dp[now^][j]<=sum&&dp[now^][j]+min(tmp,a[i])<dp[now][a[i]-tmp+sum])
                {
                    dp[now][a[i]-tmp+sum]=dp[now^][j]+min(tmp,a[i]);
                    pre[i][a[i]-tmp+sum].x=i-;
                    pre[i][a[i]-tmp+sum].y=j;
                    pre[i][a[i]-tmp+sum].f=;
                }

            }
        }
    }
}

void output()
{
    bool flag=;
    int now_i,now_j;
    for(int i=tot;i>=;i--)
    {
        for(int j=*sum;j>=;j--)
        {
            if(pre[i][j].f!=-)
            {
                now_i=i;
                now_j=j;
                flag=;
                break;
            }
        }
        if(flag) break;
    }

    if(flag==) {printf("0\n");return;}
    printf("%d\n",now_i);

    stack<int>S;

    while()
    {
        if(pre[now_i][now_j].f==-) break;
        S.push(pre[now_i][now_j].f);

        int newx,newy;
        newx=pre[now_i][now_j].x;
        newy=pre[now_i][now_j].y;
        now_i=newx;
        now_j=newy;
    }

    int k=;
    int left=;
    int right=;

    while(!S.empty())
    {
        if(S.top()==)
        {
            printf("port\n");
            left=left+a[k];k++;
        }
        else
        {
            printf("starboard\n");
            right=right+a[k];k++;
        }
        S.pop();
    }
}

int main()
{
    scanf("%d",&sum);
    init();
    read();
    work();
    output();
    return ;
}