Fzu 2255 过河【Dp】好题！

 Problem C 过河

Accept: 3    Submit: 5
Time Limit: 1000 mSec    Memory Limit : 32768 KB

 Problem Description

遥远的YS星球上，生活着两种物种名为Yellow和Star，简称物种Y和物种S。

神奇的是物种Y的体重始终为50kg，物种S的体重始终为100kg。

这天，n只生物结伴出行，在路途中，它们被一条长长的河流挡住了去路。所幸岸边停靠着一条船，船上写着：负重m kg。

YS星球上的物种有很强的时间观念，它们需要找出一种最快的方案过河：

1、要开船每次至少要有一只生物在船上

2、载着的生物总重量不能超过船的负重m

3、无论载多少只生物，船每次从岸边到另一岸边需要的时间相同，并且不考虑中间换乘时间(换句话说，要求的是最少的开船次数)

请你帮助这些Yellow、Star们，找出最少的开船次数，并且求出最少次数下的有多少种不同的方案数。

当存在某轮开船时，两个方案乘客的集合不同，认为这两个方案是不同的。

 Input

包含多组测试数据。

第一行为n, m。

接下来n个数字，数字要么是50，要么是100。

1≤n≤50

200≤m≤1000

 Output

输出两行，第一行为最少开船次数，第二行为不同方案数，由于方案数可能很大，对1000000007(10^9+7)取模后输出。

 Sample Input

3 100

50 50 100

 Sample Output

5

2

 Hint

最优的两种情况为：

1和2过去，1回来，3过去，2回来，1和2过去。

1和2过去，2回来，3过去，1回来，1和2过去。

思路：

奇数时间我们是从河这边过到那边，人数是减少的，偶数时间我们是从河那边过到河这边，人数是增加的。。

所以我们直接Dp统计方案数的话，会存在后效性。

那么我们再增加一个时间状态，就能取消后效性。

设定Dp【k】【i】【j】表示时间为k的时候 ，50重量的人剩余i个，100重量的人剩余j个的方案数。

那么状态转移分成时间奇偶两种情况即可：

①k%2==1  Dp【k】【i】【j】+=Dp【k】【i+x】【j+y】*C（i+x，x）*C（j+y，y）；

②k%2==.0 Dp【k】【i】【j】+=Dp【k】【i-x】【j-y】*C（conta-i+x，x）*C（contb-j+y，y）；

注意，我们虽然贪心最优的时候肯定是选择一个人回来，但是这里问的是方案数，所以，我们回来的时候可能还一次回来多个人。（这里我傻逼了为了得到最优 ，在Dp过程中还贪心了一波。。。。。导致最终答案有问题Wa了好多发。。。）

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define ll __int64
const ll mod=1e9+7;
ll dp[250][55][55];
ll c[55][55];
ll a[55];
void init()
{
    memset(c,0,sizeof(c));
    c[0][0]=1%mod;
    c[1][0]=1%mod;
    c[1][1]=1%mod;
    for(ll i=2;i<=52;i++)
    {
        c[i][0]=1%mod;
        for(ll j=1;j<=i;j++)
        {
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
        }
    }
}
ll C(ll x,ll y)
{
    if(y==0)return 1;
    return c[x][y];
}
int main()
{
    init();
    ll n,m;
    while(~scanf("%I64d%I64d",&n,&m))
    {
        ll conta=0,contb=0;
        for(ll i=1;i<=n;i++)
        {
            scanf("%I64d",&a[i]);
            if(a[i]==50)conta++;
            else contb++;
        }
        memset(dp,0,sizeof(dp));
        dp[0][conta][contb]=1;
        for(ll k=1;k<=130;k++)
        {
            if(k%2==1)
            {
                for(ll i=0;i<=conta;i++)
                {
                    for(ll j=0;j<=contb;j++)
                    {
                        for(ll x=0;x<=conta;x++)
                        {
                            if(x*50>m||i+x>conta)break;
                            for(ll y=0;y<=contb;y++)
                            {
                                if(x*50+y*100>m||j+y>contb)break;
                                if(x+y==0)continue;
                                dp[k][i][j]+=dp[k-1][i+x][j+y]*C(i+x,x)*C(j+y,y);
                                dp[k][i][j]%=mod;
                            }
                        }
                    }
                }
            }
            else
            {
                for(ll i=0;i<=conta;i++)
                {
                    for(ll j=0;j<=contb;j++)
                    {
                        for(ll x=0;x<=conta;x++)
                        {
                            if(x*50>m||i-x<0)break;
                            for(ll y=0;y<=contb;y++)
                            {
                                if(x*50+y*100>m||j-y<0)break;
                                if(x+y==0)continue;
                                dp[k][i][j]+=dp[k-1][i-x][j-y]*C(conta-i+x,x)*C(contb-j+y,y);
                                dp[k][i][j]%=mod;
                            }
                        }
                    }
                }
            }
        }
        for(ll i=1;i<=130;i++)
        {
            if(dp[i][0][0]>0)
            {
                printf("%I64d\n%I64d\n",i,dp[i][0][0]);
                break;
            }
        }
    }
}