C++博客-Dreams-随笔分类-DP

hdu 2372 El Dorado

DreamSky — Thu, 14 May 2009 11:49:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=2372

//1367544 2009-05-14 19:45:52 Accepted 2372 0MS 352K 678 B C++ no way
#include<iostream>
using namespace std;
const int N = 105;
__int64 dp[N][N];//dp[i][j]表示一第i个元素结尾的长度为j的单调序列个数
int main()
{
    int n,k;
    while(cin>>n>>k)
    {
        if(n == 0 && k == 0)
            break;
        int i,j,p;
        __int64 num[N],temp;
        for(i=1;i<=n;i++)
            scanf("%I64d",&num[i]);

        for(j=2;j<=k;j++)
            for(i=j;i<=n;i++)
                dp[i][j] = 0;

        for(i=1;i<=n;i++)
            dp[i][1] = 1;

        for(j=2;j<=k;j++)
            for(i=j;i<=n;i++)
            {
                for(p=j-1;p<i;p++)
                    if(num[i] > num[p])
                        dp[i][j] += dp[p][j-1];
            }

        for(temp=0,i=k;i<=n;i++)
            temp += dp[i][k];

        printf("%I64d\n",temp);
    }
    return 0;
}

DreamSky 2009-05-14 19:49 发表评论

01-package

DreamSky — Fri, 08 May 2009 13:48:00 GMT

http://info.zjfc.edu.cn/acm/contest/contest_problemDetail.aspx?pid=1005&cid=32
题目描述：
    给定一个背包的容量k，给定n个物品的体积和价值，物品不可分割，将n个物品中选若干个物品放入背包，求背包内物品的最大价值总和，在价值总和最大的前提下求背包内的最小物品个数c。
输出描述：
    第一行是一个整数t，表示测试数据的组数t。
   对于每组测试数据，第一行是两个整数n和k,表示物品的个数和背包的容量；
    接下来n行，每行两个整数，分别是物品的价值和体积。
输出描述：
输出背包内物品的最大价值v，在价值最大的前提下求背包内的最小物品个数c，中间用一个空格隔开。
样例输入：
1
3 10
4 5
6 5
10 10
样例输出：
10 1
作者：
xiewenxiu

//16829 2009-05-08 19:58:40 1005  Accepted 765MS 464K Visual C++ liyunsong
#include<iostream>
using namespace std;
struct Node
{
    int ns;//最小物品数
    int vs;//最大价值
}dp[2001];
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int v[2001],w[2001];
        int n,c;
        int i,j;
        cin>>n>>c;
        for(i=1;i<=n;i++)
            scanf("%d%d",&v[i],&w[i]);
        for(j=0 ; j < w[1];j++)
        {
            dp[j].vs = 0;
            dp[j].ns = 0;
        }
        for(; j <= c;j++)
        {
            dp[j].vs = v[1];
            dp[j].ns = 1;
        }
        for(i=2;i<=n;i++)
        {
            for(j=c;j>=w[i];j--)
            {
                if(dp[j].vs < dp[j-w[i]].vs + v[i])
                {
                    dp[j].vs = dp[j-w[i]].vs + v[i];
                    dp[j].ns = dp[j-w[i]].ns + 1;
                }
                else if(dp[j].vs == dp[j-w[i]].vs + v[i] && dp[j].ns > dp[j-w[i]].ns + 1)
                    dp[j].ns = dp[j-w[i]].ns + 1;
            }
        }
        cout<<dp[c].vs<<" "<<dp[c].ns<<endl;
    }
    return 0;
}

DreamSky 2009-05-08 21:48 发表评论

zju 1883 Tight Words

DreamSky — Fri, 08 May 2009 07:24:00 GMT

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=883

#include<stdio.h>
#include<math.h>
int main()
{
    int k,n;
    while(scanf("%d%d",&k,&n)!=EOF)
    {
        int i,t;
        double f[11][101];
        for(i=0;i<=k;i++)
            f[i][1] = 1;
        for(t = 2; t<=n;t++)
        {
            for(i=0;i<=k;i++)
            {
                f[i][t] = f[i][t-1] ;
                if(i+1<=k)
                   f[i][t] += f[i+1][t-1];
                if(i-1>=0)
                   f[i][t] += f[i-1][t-1];
            }
        }
        double sum=0.0,p;
        for(i=0;i<=k;i++)
            sum += f[i][n];
        p = sum / pow((double)(k+1),(double)n);
        printf("%.5lf\n",p*100);
    }
    return 0;
}

DreamSky 2009-05-08 15:24 发表评论

zju 3201 Tree of Tree

DreamSky — Thu, 07 May 2009 06:36:00 GMT

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=3326

//1858391 2009-05-07 14:27:05 Accepted  3201 C++ 0 228 aslys
#include<iostream>
#include<vector>
using namespace std;
const int N = 105;
int n,k;
int f[N][N],w[N];
vector<int>sum[N];
//学习使用vector
void dp(int x,int y)
{
    int i,j,temp[N];
    memset(temp,0,sizeof(temp));
    for(i=1;i<=k;i++)
    {
        for(j=1;j<=i;j++)
            if(f[x][j] + f[y][i-j] > temp[i])//合并出最大的
                temp[i] = f[x][j] + f[y][i-j];
    }
    for(i=1;i<=k;i++)
        if(f[x][i] < temp[i])
            f[x][i] = temp[i];
}

void dfs(int now,int pre)
{
    int len = sum[now].size();
    int i,cnt = 0;
    for(i=0;i<len;i++)
    {
        if(sum[now][i] != pre)
        {
            dfs(sum[now][i],now);
            f[now][1] = w[now];
            dp(now,sum[now][i]);
            cnt ++;
        }
    }
    if(cnt == 0)
        f[now][1] = w[now];
}

int main()
{
    while(scanf("%d%d",&n,&k)!=EOF)
    {
        int i,a,b,maxs;
        for(i=0;i<n;i++)
        {
            scanf("%d",&w[i]);
            sum[i].clear();
        }
        for(i=1;i<n;i++)
        {
            scanf("%d%d",&a,&b);
            sum[a].push_back(b);
            sum[b].push_back(a);
        }
        maxs = 0;
        memset(f,0,sizeof(f));
        dfs(0,-1);//从任一点出发都行
        for(i=0;i<n;i++)
            if(maxs < f[i][k])
                maxs = f[i][k];
            printf("%d\n",maxs);
    }
    return 0;
}

DreamSky 2009-05-07 14:36 发表评论

zju 2852 Deck of Cards

DreamSky — Wed, 06 May 2009 00:48:00 GMT

摘要: http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1852 #includeusing namespace std;int dp[201][22][22][22];//dp[k][r][s][t]表示前k张card里面三堆分别为r点、s点、t点时的最优值int&nbs... 阅读全文

DreamSky 2009-05-06 08:48 发表评论

hdu 2191 悼念512汶川大地震遇难同胞——珍惜现在，感恩生活

DreamSky — Mon, 04 May 2009 12:49:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=2191

//1337971 2009-05-04 20:46:09 Accepted 2191 15MS 268K 857 B C++ no way
#include<iostream>
using namespace std;
struct Node
{
    int p;//价格
    int w;//重量
    int c;//数量
};
int main()
{
    int t,N,m;
    cin>>t;
    while(t--)
    {
        int i,j,k,temp;
        int dp[101];//dp[j]表示钱为j时能买的最大重量
        Node infor[101];
        cin>>N>>m;
        for(i=1;i<=m;i++)
            scanf("%d%d%d",&infor[i].p,&infor[i].w,&infor[i].c);
        for(j=0;j<infor[1].p;j++)
            dp[j] = 0;
        for(;j<=N;j++)
        {
            temp = j/infor[1].p;
            if(temp <= infor[1].c)
                dp[j] = temp * infor[1].w;
            else
                dp[j] = infor[1].c * infor[1].w;
        }

        for(i=2;i<=m;i++)//米的种类
        {
            for(j=N;j>=infor[i].p;j--)//钱的数量
            {
                for(k = 0;k<=infor[i].c;k++)//米的代数
                {
                    temp = k * infor[i].p;
                    if(temp <= j )
                    {
                        if(dp[j-temp]+k*infor[i].w > dp[j])
                            dp[j] = dp[j-temp]+k*infor[i].w;
                    }
                    else
                        break;
                }
            }
        }
        cout<<dp[N]<<endl;
    }
    return 0;
}

DreamSky 2009-05-04 20:49 发表评论

hdu 2765 Recursively Palindromic Partitions

DreamSky — Sat, 02 May 2009 10:58:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=2765

//1331246 2009-05-02 18:53:43 Accepted 2765 15MS 320K 882 B C++ no way
#include<iostream>
using namespace std;
const int MAX = 2140000000;
int f[1000001];

void dfs(int p)
{
    int i,sum = 1,temp;
    if(p%2 == 1)
    {
        for(i=1;i<p;i+=2)
        {
            temp = (p-i)/2;
            if(f[temp]==MAX)
                dfs(temp);
            sum += f[temp];
        }
    }
    else //even
    {
        for(i=0;i<p;i+=2)
        {
            temp = (p-i)/2;
            if(f[temp]==MAX)
                dfs(temp);
            sum += f[temp];
        }
    }
    f[p] = sum;
}

int main()
{
    int t,n,i,temp;
    for(i=0;i<=10000;i++)
        f[i] = MAX;
    f[0] = 0;
    f[1] = 1;
    f[2] = 2;
    f[3] = 2;
    int cas = 1;
    cin>>t;
    while(t--)
    {
        cin>>n;
        if(f[n] == MAX)
            dfs(n);
        cout<<cas++<<" "<<f[n]<<endl;
    }
    return 0;
}

DreamSky 2009-05-02 18:58 发表评论

vijos 1313 金明的预算方案

DreamSky — Mon, 27 Apr 2009 09:41:00 GMT

摘要: http://www.vijos.cn/Problem_Show.asp?id=1313 #includeusing namespace std;struct Node{ int v;//该物品的价格 int p;//该物品... 阅读全文

DreamSky 2009-04-27 17:41 发表评论

vijos 1133 装箱问题

DreamSky — Mon, 27 Apr 2009 06:23:00 GMT

http://www.vijos.cn/Problem_Show.asp?id=1133

#include<iostream>
using namespace std;
int MIN(int x,int y)
{
    if(x<y)
        return x;
    else
        return y;
}
int main()
{
    int i,j,c,n;//c箱子容量，n表示物品数量
    while(cin>>c)
    {
        int v[31],dp[20001];//v[i]表示第i个物品的体积
        //dp[j]表示容量为j时的最优值
        scanf("%d",&n);
        for(i=1;i<=n;i++)
            scanf("%d",&v[i]);
        for(j=0;j<v[1];j++)
            dp[j] = c;
        for(;j<=c;j++)
            dp[j] = c-v[1];
        for(i=2;i<=n;i++)
            for(j=c;j>=v[i];j--)
                dp[j] = MIN(dp[j],dp[j-v[i]] - v[i]);
        cout<<dp[c]<<endl;
    }
}

DreamSky 2009-04-27 14:23 发表评论

vijos 1317 开心的金明

DreamSky — Mon, 27 Apr 2009 05:56:00 GMT

http://www.vijos.cn/Problem_Show.asp?id=1317

#include<iostream>
using namespace std;
int dp[30001];//dp[i]背包容量为i时从1……x个物品中选取的最优值
int MAX(int x,int y)
{
    if(x>y)
        return x;
    else
        return y;
}
int main()
{
    int c,m;//c总共钱数，m物品数
    while(cin>>c>>m)
    {
        int i,j;
        int v[25],p[25];//v[i]、p[i]分别表示第 i 件物品的价格和重要度
        for(i=1;i<=m;i++)
            scanf("%d%d",&v[i],&p[i]);

        for(j=0;j<v[1];j++)
            dp[j] = 0;
        for(;j<=c;j++)
            dp[j] = v[1] * p[1];

        for(i=2;i<=m;i++)
        {
            for(j=c;j>=v[i];j--)
                dp[j] = MAX(dp[j],dp[j-v[i]] + v[i] * p[i]);
        }
        cout<<dp[c]<<endl;
    }
    return 0;
}

DreamSky 2009-04-27 13:56 发表评论

hdu 2670 Girl Love Value

DreamSky — Sun, 26 Apr 2009 12:09:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=2670

//1311986 2009-04-26 19:58:30 Accepted 2670 78MS 4180K 934 B C++ no way
#include<iostream>
#include<algorithm>
#define MAX(x,y) x>y?x:y
using namespace std;
struct node
{
    int love;
    int dece;
};
bool comp(node a,node b)
{
    return a.dece > b.dece;
}
int dp[1001][1001];//dp[i][j] 表示前i个人中选j个的最优值
int lvs[1001][1001];
int main()
{
    int n,k;
    while(scanf("%d%d",&n,&k)!=EOF)
    {
        int i,j,maxs;
        node infor[1001];
        for(i=1;i<=n;i++)
            scanf("%d",&infor[i].love);
        for(i=1;i<=n;i++)
            scanf("%d",&infor[i].dece);
        sort(infor+1,infor+n+1,comp);
        ////////////////////////////////////////////////////////////////
        for(i=1;i<=n;i++)
            for(j=1;j<=k;j++)
                lvs[i][j] = infor[i].love - (j-1)*infor[i].dece;
        //这样处理时间并未节约，相反还浪费了空间~？？？？
        ////////////////////////////////////////////////////////////////
        maxs = 0;
        for(i=1;i<=n;i++)
        {
            if(infor[i].love > maxs)
                maxs = infor[i].love;
            dp[i][1] = maxs;
        }
        for(i=2;i<=n;i++)
        {
            for(j=2;j<=k && j<=i;j++)
            {
                if(i-1>=j)
                    dp[i][j] = MAX(dp[i-1][j],dp[i-1][j-1] + lvs[i][j]);
                else
                    dp[i][j] = dp[i-1][j-1] + lvs[i][j];
            }
        }
        cout<<dp[n][k]<<endl;
    }
    return 0;
}

DreamSky 2009-04-26 20:09 发表评论

hdu 1074 Doing Homework

DreamSky — Sun, 26 Apr 2009 11:30:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=1074

//1311853 2009-04-26 19:16:12 Accepted 1074 281MS 524K 1323 B C++ no way
#include<iostream>
#include<string>
using namespace std;
struct Node
{
    string course;
    int deadline;
    int need;
}infor[16];
int n;
int binary[16] = {1,2,4,8,16,32,64,128,256,512,1024,2048,4096,8192,16384,32768};
int flag[65536];
bool mark[16];
string str[16],outs[16];
void dfs(int days,int cost,int sum,int num)//num表示选择的作业数
{ //days表示写作业的开始日期，cost表示前面的花费，sum记录作业是否完成情况
    int i,temp;
    if(num == n)
    {
        if(flag[sum] == cost)
        {
            flag[sum] = cost;
            for(i=0;i<num;i++)
                outs[i] = str[i];
        }
        return ;
    }
    for(i=1;i<=n;i++)
    {
        if(mark[i] == false)
        {
            mark[i] = true;
            sum += binary[i];//要写第i门作业
           ///////////////////////////////////////////////////
            temp = days + infor[i].need - infor[i].deadline;
            if(temp < 0)
                temp = cost;
            else
                temp += cost;
           //  第i门作业完成后的代价 temp
           //////////////////////////////////////
            if(flag[sum] > temp)
            {
                flag[sum] = temp; //记录状态

                str[num++] = infor[i].course;

                dfs(days + infor[i].need,temp,sum,num);

                num --;
            }
            sum -= binary[i];
            mark[i] = false;
        }
    }
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int i,st=0;
        cin>>n;
        for(i=1;i<=n;i++)
            cin>>infor[i].course>>infor[i].deadline>>infor[i].need;
        for(i=0;i<65536;i++)
            flag[i] = 1000000;
        memset(mark,false,sizeof(mark));//标记状态
        dfs(0,0,0,0);
        for(i=1;i<=n;i++)
            st += binary[i];  //结果存放在下标为st的flag[st]值中
        cout<<flag[st]<<endl;
        for(i=0;i<n;i++)
            cout<<outs[i]<<endl;
    }
    return 0;
}

DreamSky 2009-04-26 19:30 发表评论

hdu 2059 龟兔赛跑

DreamSky — Sun, 26 Apr 2009 10:51:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=2059

#include<iostream>
using namespace std;
struct Node
{
    double tim;
    double ls;
};
int main()
{
    double S;
    while(cin>>S)
    {
        int n,i,j;
        double L,time,rv,tortoiseC,tortoiseQ,pos[103];
        Node dp[103],temp;
        double d,d1,d2,t1,t2;
        cin>>n>>L>>time;
        cin>>rv>>tortoiseQ>>tortoiseC;
        for(i=1;i<=n;i++)
            cin>>pos[i];
        if(tortoiseQ<=tortoiseC)//骑车不如蹬车
        {
            if(S/tortoiseC>S/rv)
                cout<<"Good job,rabbit!"<<endl;
            else
                cout<<"What a pity rabbit!"<<endl;
            continue;
        }
        dp[0].tim=0;
        dp[0].ls=L;
        pos[0]=0;//一定要赋初值，否则其默认初值为负数，因为是double类型
        pos[n+1]=S;
        for(i=1;i<=n+1;i++)
        {
            dp[i].tim=100000000;
            dp[i].ls=0;
            for(j=0;j<i;j++)//从出发点开始到i站之前选择一个最优的
            {
                d=pos[i]-pos[j];
                if(dp[j].ls>=d)//无需充电
                {
                    temp.tim=dp[j].tim+d/tortoiseQ;
                    temp.ls=dp[j].ls-d;
                }
                else
                {
                    t1=dp[j].tim+dp[j].ls/tortoiseQ+(d-dp[j].ls)/tortoiseC;
                    d1=0;
                    if(L>=d)//充电够用情况
                    {
                        t2=dp[j].tim+time+d/tortoiseQ;
                        d2=L-d;
                    }
                    else//充电不够用情况
                    {
                        t2=dp[j].tim+time+L/tortoiseQ+(d-L)/tortoiseC;
                        d2=0;
                    }

                    if(t1<t2)//选择时间最少的
                    {
                        temp.tim=t1;
                        temp.ls=d1;
                    }
                    else
                    {
                        temp.tim=t2;
                        temp.ls=d2;
                    }
                }
                if(temp.tim < dp[i].tim)
                {
                    dp[i].tim=temp.tim;
                    dp[i].ls=temp.ls;
                }
            }//for(j=0;j
        }//for(i=1;i<=n+1;i++)
        if(dp[n+1].tim>S/rv)
            cout<<"Good job,rabbit!"<<endl;
        else
            cout<<"What a pity rabbit!"<<endl;
    }
    return 0;
}

DreamSky 2009-04-26 18:51 发表评论

hdu 1114 Piggy-Bank

DreamSky — Sat, 25 Apr 2009 12:22:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=1114
原始代码：

//1308978 2009-04-25 19:39:15 Accepted 1114 546MS 312K 730 B C++ no way
#include<iostream>
using namespace std;
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int i,j,k;
        int c,e,f,n;
        int w[501],v[501],dp[10001];
        scanf("%d%d",&e,&f);
        c = f - e;
        scanf("%d",&n);
        for(i=1;i<=n;i++)
            scanf("%d%d",&v[i],&w[i]);

        for(i=1;i<=c;i++)
            dp[i] = 1000000;
        for(i=0; i*w[1] <= c; i++)
            dp[i*w[1]] = i*v[1];

        for(i=2;i<=n;i++)
        {
            for(k=0;k*w[i]<=c;k++)//所取物品数量
            {
                for(j=k*w[i];j<=c;j++)
                {
                    if(dp[j] > dp[j-k*w[i]] + k*v[i])
                        dp[j] = dp[j-k*w[i]] + k*v[i];
                }
            }
        }
        if(dp[c] == 1000000)
            cout<<"This is impossible."<<endl;
        else
            cout<<"The minimum amount of money in the piggy-bank is "<<dp[c]<<"."<<endl;
    }
    return 0;
}

优化代码：

//1309096 2009-04-25 20:12:49 Accepted 1114 93MS 384K 793 B C++ no way
#include<iostream>
struct node
{
    int w;
    int v;
}tt[10001];
using namespace std;
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int i,j,k=1;
        int c,e,f,n;
        int dp[10001];//dp[i]表容量为i的时候所装东西的最小价值//w[10001],v[10001],
        scanf("%d%d",&e,&f);
        c = f - e;
        scanf("%d",&n);
        for(i=1;i<=n;i++)
        {
            scanf("%d%d",&e,&f);
            tt[k].w = f;
            tt[k].v = e;
            while(1)
            {
                k ++;
                if(tt[k-1].w*2 > c)
                    break;
                tt[k].w = tt[k-1].w * 2;
                tt[k].v = tt[k-1].v * 2;
            }
        }
        for(i=1;i<=c;i++)
            dp[i] = 1000000;
        dp[0] = 0;//背包容量为 0 啥也不能装
        for(i=1;i < k;i++)
        {
            for(j=tt[i].w;j<=c;j++)
            {
                if(dp[j] > dp[j-tt[i].w] + tt[i].v)
                    dp[j] = dp[j-tt[i].w] + tt[i].v;
            }
        }
        if(dp[c] == 1000000)
            cout<<"This is impossible."<<endl;
        else
            cout<<"The minimum amount of money in the piggy-bank is "<<dp[c]<<"."<<endl;
    }
    return 0;
}

DreamSky 2009-04-25 20:22 发表评论

va家族的等级制

DreamSky — Fri, 24 Apr 2009 13:52:00 GMT

http://info.zjfc.edu.cn/acm/contest/contest_problemDetail.aspx?pid=1002&cid=29

//15344 2009-04-24 21:46:41 1002  Accepted 125MS 7984K Visual C++ xredman
#include<iostream>
#include<stdio.h>
#include<string>
using namespace std;
bool dp[2001][2001];//dp[i][j]标志从i到j这段串是否为回文
char str[2001];
int f[2001];//f[i]表示从1到i最少的回文数
int main()
{
    while(scanf("%s",str+ 1)!=EOF)
    {
        int len,i,j,mins;
        len = strlen(str+1);

        for(i=1;i<=len;i++)
            for(j=1;j<=len;j++)
                dp[i][j] = 0;

        for(i=1;i<=len;i++)
            dp[i][i] = 1;

        for(i=2;i<=len;i++)
            if(str[i] == str[i-1])
                dp[i-1][i] = 1;

        for(i=3;i<=len;i++)//长度
        {
            for(j=1;j<=len-i+1;j++)
            {
                if(str[j] == str[j+i-1] && dp[j+1][j+i-2] == 1)
                    dp[j][j+i-1] = 1;
            }
        }
        if(dp[1][len] == 1)
            printf("1\n");
        else
        {
            f[1] = 1;
            for(i=2;i<=len;i++)
            {
                mins = i;
                if(dp[1][i] == 1)
                {
                    f[i] = 1;
                    continue;
                }
                for(j=1;j<i;j++)
                {
                    if(dp[j+1][i] == 1 && mins > f[j] + 1)
                        mins = f[j] + 1;
                }
                f[i] = mins;
            }
            printf("%d\n",f[len]);
        }
    }
    return 0;
}

DreamSky 2009-04-24 21:52 发表评论

hdu 2512 一卡通大冒险

DreamSky — Mon, 20 Apr 2009 08:38:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=2512

//1290997 2009-04-20 16:28:41 Accepted 2512 62MS 14012K 545 B C++ no way
#include<iostream>
using namespace std;
const int N=2001;
int dp[N][N];//dp[i][j]表示i张卡片分成j堆的情况数
int main()
{
    int t,i,j;
    int num[N];

    for(i=1;i<=N-1;i++)
        dp[i][i] = dp[i][1] = 1;

    for(j=2;j<N-1;j++)//堆数
        for(i=j+1;i<=N-1;i++)//卡片数
            dp[i][j] =( (dp[i-1][j] * j)%10000 + dp[i-1][j-1] ) % 10000;

    for(i=1;i<=N-1;i++)
    {
        num[i] = 0;
        for(j=1;j<=i;j++)
        {
            num[i] += dp[i][j];
            num[i] %= 10000;
        }
        num[i] %= 1000;
    }
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&i);
        printf("%d\n",num[i]);
    }
    return 0;
}

DreamSky 2009-04-20 16:38 发表评论

hdu 2182 Frog

DreamSky — Fri, 17 Apr 2009 14:10:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=2182

//1281694 2009-04-17 22:05:25 Accepted 2182 0MS 304K 718 B C++ no way
#include<iostream>
#include<map>
using namespace std;
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n,a,b,k;
        int dp[101][101],num[101];
        int i,j,p,M,maxs;
        scanf("%d%d%d%d",&n,&a,&b,&k);
        for(i=1;i<=n;i++)
            scanf("%d",&num[i]);

        for(i=1;i<=n;i++) //跳0步肯定是0
            dp[i][0] = 0;

        M = dp[1][0] = num[1];//第1个位置的0步状态是有虫的

        for(j=1;j<=k;j++)//跳跃的步数
        {
            for(i=2;i<=n;i++)//青蛙所在的位置
            {
                maxs = -1;
                for(p=1;p<i;p++) //跟其前面的位置有关
                {
                    if(p+a<= i && p+b >=i && maxs < dp[p][j-1])
                        maxs = dp[p][j-1];
                }
                if(maxs == -1)//第j次跳不能跳到j这个位置
                    dp[i][j] = 0;
                else
                    dp[i][j] = maxs + num[i];
                if(M < dp[i][j])
                    M = dp[i][j];
            }
        }
        cout<<M<<endl;
    }
    return 0;
}

DreamSky 2009-04-17 22:10 发表评论

NK 1137 （hutc 1036）Pebble Merging（石子合并问题）

DreamSky — Fri, 17 Apr 2009 05:05:00 GMT

http://acm.nankai.edu.cn/p1137.html
http://172.20.16.106/JudgeOnline/showproblem?problem_id=1036

//7562 liyunsong B Accepted     G++   2009-04-17 12:57:36
#include<iostream>
using namespace std;
const int N = 101;
int dpmax[2*N][2*N];//dpmax[i][j]表示从第i对石子开始连续顺时针合并j堆的最大耗费
int dpmin[2*N][2*N];//dpmin[i][j]表示从第i对石子开始连续顺时针合并j堆的最小耗费
int mins[2*N][2*N];
int maxs[2*N][2*N];
int num[2*N];

int main()
{
    int n,i,j,k,r;
    while(cin>>n)
    {
        int t;
        for(i=1;i<=n;i++)
            scanf("%d",&num[i]);

        if(n==1)
        {
            printf("0\n0\n");
            continue;
        }

        for(i=1;i<=n;i++)
            num[i+n] = num[i];

        for(i=1;i<=2*n;i++)
            for(j=1;j<=2*n;j++)
            {
                dpmin[i][j] = 10000000;
                dpmax[i][j] = 0;
            }

        for(i=1;i<=2*n;i++)
        {
            dpmin[i][1] = dpmax[i][1] = 0;
            dpmin[i][i] = dpmax[i][i] = 0;
            mins[i][i] = maxs[i][i] = num[i];
        }

        for(r = 2;r<=n;r++)//合并堆数
        {
            for(i=1;i<2*n-r+1;i++)//起点
            {
                j = i+r-1;
                dpmin[i][j] = dpmin[i+1][j] + num[i] + mins[i+1][j];
                mins[i][j] = num[i] + mins[i+1][j];
                dpmax[i][j] = dpmax[i+1][j] + num[i] + maxs[i+1][j];
                maxs[i][j] = num[i] + maxs[i+1][j];
                for(k=i+1;k<j;k++)
                {
                    t = mins[i][k] + mins[k+1][j] + dpmin[i][k] + dpmin[k+1][j];
                    if(t < dpmin[i][j])
                    {
                        dpmin[i][j] = t;
                        mins[i][j] = mins[i][k] + mins[k+1][j];
                    }

                    t = maxs[i][k] + maxs[k+1][j] + dpmax[i][k] + dpmax[k+1][j];
                    if(t > dpmax[i][j])
                    {
                        dpmax[i][j] = t;
                        maxs[i][j] = maxs[i][k] + maxs[k+1][j];
                    }
                }
            }
        }
        int MAX = dpmax[1][n];
        int MIN = dpmin[1][n];
        for(i=2;i<=n;i++)
        {
            if(MAX < dpmax[i][i+n-1])
                MAX = dpmax[i][i+n-1];
            if(MIN > dpmin[i][i+n-1])
                MIN = dpmin[i][i+n-1];
        }
        printf("%d\n%d\n",MIN,MAX);
    }
    return 0;
}

DreamSky 2009-04-17 13:05 发表评论

hutc 1035 编辑距离问题

DreamSky — Thu, 16 Apr 2009 06:37:00 GMT

http://172.20.16.106/JudgeOnline/showproblem?problem_id=1035
http://acm.nankai.edu.cn/p1131.html
DP代码

#include<iostream>
#include<string>
using namespace std;
const int N = 2001;//用1001 Runtime Error，不要开得太小
int dp[N][N];//dp[i][j]表示串s的前i个字符匹配到串t的前j个字符的最小步数
int lens,lent;
char s[N],t[N];

int main()
{
    int i,j;
    while(scanf("%s%s",s,t)!=EOF)
    {
        lens = strlen(s);
        lent = strlen(t);

        for(i=0;i<=lens;i++)
            dp[i][0] = i;
        for(j=0;j<=lent;j++)
            dp[0][j] = j;//变换到第j个字符最多用j步

        for(i=1;i<=lens;i++)
        {
            for(j=1;j<=lent;j++)//目标串指针
            {
                if(s[i-1] == t[j-1])
                    dp[i][j] = dp[i-1][j-1];
                else //更改
                    dp[i][j] = dp[i-1][j-1] + 1;

                if(dp[i][j] > dp[i-1][j] + 1)//删除第i个
                    dp[i][j] = dp[i-1][j] + 1;

                if(dp[i][j] > dp[i][j-1] + 1)//插入第i个
                    dp[i][j] = dp[i][j-1] + 1;
            }
        }
        printf("%d\n",dp[lens][lent]);
    }
    return 0;
}

DFS代码（在南开JudgeOnline超时）

#include<iostream>
#include<string>
using namespace std;
const int M = 10000000;
const int N = 2001;//用1001 Runtime Error，不要开得太小
int dp[N][N];
int lens,lent;
char s[N],t[N];
int dfs(int i,int j)
{
    int p,q,r,mins=M;

    if(i == lens && j == lent )
        return 0;
    else if(i == lens && j != lent )
        return lent - j;
    else if(i != lens && j == lent )
        return lens - i;

    if(dp[i][j] != M) //做过的没有必要再做
        return dp[i][j];

    if(s[i] == t[j])
        mins = dfs(i+1,j+1);//不用改变
    else
    {
        p = dfs(i+1,j+1) + 1;//修改
        q = dfs(i+1,j) + 1;//删除
        r = dfs(i,j+1) + 1;//插入
        if(p > q)
            mins = q;
        else
            mins = p;
        if(mins > r)
            mins = r;
    }
    dp[i][j] = mins;
    return mins;
}

int main()
{
    int i,j;
    while(scanf("%s%s",s,t)!=EOF)
    {
        lens = strlen(s);
        lent = strlen(t);
        for(i=0;i<lens;i++)
            for(j=0;j<lent;j++)
                dp[i][j] = M;
        dfs(0,0);
        printf("%d\n",dp[0][0]);
    }
    return 0;
}

DreamSky 2009-04-16 14:37 发表评论

zju 1234 Chopsticks

DreamSky — Thu, 16 Apr 2009 02:04:00 GMT

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=234
http://acm.hdu.edu.cn/showproblem.php?pid=1500

#include<iostream>
using namespace std;
long dp[5001][1009];//dp[i][j]表示从第i个筷子到第n个筷子中取j对的最优值
int main()
{
    long t,k,n,i,j;
    scanf("%ld",&t);
    while(t--)
    {
        long chopsticks[5001],num[5001];
        //cin>>k>>n;
        scanf("%ld%ld",&k,&n);
        k += 8;
        for(i=1;i<=n;i++)
            scanf("%ld",&chopsticks[i]);

        for(i=1;i<n;i++)
            num[i] = (chopsticks[i+1] - chopsticks[i]) * (chopsticks[i+1] - chopsticks[i]);

        for(i=0;i<=n;i++)
            for(j=0;j<=k;j++)
                dp[i][j] = 1024000000;//初值一定要足够大

        for(i=0;i<=n;i++)//取0对筷子
            dp[i][0] = 0;

        for(i=n-2;i>=1;i--)//从后往前
        {
            for(j=1;3*j<=(n-i+1);j++)
            {
                dp[i][j] = dp[i+1][j];
                if(dp[i][j] > dp[i+2][j-1] + num[i])
                    dp[i][j] = dp[i+2][j-1] + num[i];
            }
        }
        //cout<
        printf("%ld\n",dp[1][k]);
    }
    return 0;
}

DreamSky 2009-04-16 10:04 发表评论

hdu 1421 搬寝室

DreamSky — Wed, 15 Apr 2009 11:01:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=1421

//1274154 2009-04-15 18:46:58 Accepted 1421 187MS 8100K 793 B C++ no way
/*
8 3
1 1 5 6 7 9 11 13
*/
#include<iostream>
#include<algorithm>
#define MIN(x,y) x
using namespace std;
bool comp(int a,int b)
{
    return a<b;
}
int dp[2001][1001];//dp[i][j]表示从1……i个物品中选取j对
int num[2001];
int main()
{
    int n,k,i,j;
    while(cin>>n>>k)
    {
        int weights[2001];
        for(i=1;i<=n;i++)
            cin>>weights[i];

        sort(weights+1,weights+n+1,comp);

        for(i=1;i<n;i++)
            num[i] = (weights[i+1] - weights[i]) * (weights[i+1] - weights[i]);

        for(i=0;i<=n;i++)//物品
            for(j=1;j<=k;j++)//所取物品对数
                dp[i][j] = 100000000;//要足够大

        for(i=0;i<=n;i++)
            dp[i][0] = 0;

        for(i=2;i<=n;i++)//物品
        {
            for(j=1;2*j<=i;j++)//所取物品对数
            {
                dp[i][j] = MIN(dp[i-2][j-1] + num[i-1] , dp[i-1][j]);
            }
        }
        cout<<dp[n][k]<<endl;
    }
    return 0;
}

DreamSky 2009-04-15 19:01 发表评论

hdu 1159 Common Subsequence

DreamSky — Wed, 15 Apr 2009 05:16:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=1159

#include<iostream>
#include<string>
using namespace std;
const int M = 1001;
char s[M],t[M];
int dp[M][M];//dp[i][j]表示串s的前i个与串t的前j个的最大匹配
//int track[M][M];//保存路径
/*
void LCS(int i,int j)//输出公共子串
{
    if(i==0 || j==0)
        return ;
    if(track[i][j] == 1)
    {
        LCS(i-1,j-1);
        printf("%c ",s[i-1]);
    }
    else if(track[i][j] == 2)
        LCS(i-1,j);
    else
        LCS(i,j-1);
}
*/
int main()
{
    int i,j,lens,lent;
    while(scanf("%s%s",s,t)!=EOF)
    {
        lens = strlen(s);
        lent = strlen(t);

        for(i=0;i<=lens;i++)
            dp[i][0] = 0;

        for(i=0;i<=lent;i++)
            dp[0][i] = 0;

        for(i=1;i<=lens;i++)
        {
            for(j=1;j<=lent;j++)
            {
                if(s[i-1] == t[j-1])
                {
                    dp[i][j] = dp[i-1][j-1] + 1;
                    //track[i][j] = 1;
                }
                else if(dp[i][j-1] > dp[i-1][j])
                {
                    dp[i][j] = dp[i][j-1];
                    //track[i][j] = 3;
                }
                else
                {
                    dp[i][j] = dp[i-1][j];
                    //track[i][j] = 2;
                }
            }
        }
        cout<<dp[lens][lent]<<endl;
        //LCS(lens,lent);
        //cout<
    }
    return 0;
}

DreamSky 2009-04-15 13:16 发表评论

zju 1558 Euro Efficiency

DreamSky — Fri, 10 Apr 2009 12:45:00 GMT

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=558

//1828057 2009-04-10 20:14:44 Accepted  1558 C++ 0 184 aslys
//很悬
#include<iostream>
#include<queue>
using namespace std;
int corns[7],f[101];

void fun()
{
    int i,k;
    queue<int>Q;
    for(i=1;i<=6;i++)
        Q.push(corns[i]);
    while(!Q.empty())
    {
        k = Q.front();
        Q.pop();
        for(i=1;i<k;i++)
        {
            if(f[k-i] > f[k]+f[i])
            {
                f[k-i] = f[k]+f[i];
                Q.push(k-i);
            }
        }
        for(i=k+1;i<=100;i++)
        {
            if(f[i-k] > f[k]+f[i])
            {
                f[i-k] = f[k]+f[i];
                Q.push(i-k);
            }
        }
        for(i=1;i<=100;i++)
        {
            if(i+k <=100)
            {
                if(f[i+k] > f[k]+f[i])
                {
                    f[i+k] = f[k]+f[i];
                    Q.push(i+k);
                }
            }
            else
                break;
        }
    }
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int i,sum=0,maxs=1;
        cin>>corns[1]>>corns[2]>>corns[3]>>corns[4]>>corns[5]>>corns[6];
        //六种货币
        for(i=1;i<=100;i++)
            f[i]=i;
        for(i=1;i<=6;i++)
            f[corns[i]]=1;
        fun();
        for(i=1;i<=100;i++)
        {
            if(f[i] > maxs)
                maxs = f[i];
            sum+=f[i];
        }
        printf("%.2lf %d\n",(double)sum/100.0,maxs);
    }
    return 0;
}

DreamSky 2009-04-10 20:45 发表评论

hdu 1078 FatMouse and Cheese

DreamSky — Fri, 03 Apr 2009 00:46:00 GMT

http://acm.hdu.edu.cn/showproblem.php?pid=1078

#include<iostream>
using namespace std;
int gra[101][101];
int dp[101][101];//表示从i行j列能够获得的最大收益
int n,k;
int dir[4][2]={-1,0,1,0,0,-1,0,1};

int dfs(int i,int j) //记忆法搜索
{
    int s,t,p,q,temp,m=0;
    if(dp[i][j]>0)//走过的就没必要走了
        return dp[i][j];
    for(q=0;q<4;q++)
    {
        for(p=1;p<=k;p++)//一次可以走的步数
        {
            s = i + dir[q][0]*p;
            t = j + dir[q][1]*p;
             if(s>=0 && s<n && t>=0 && t<n && gra[i][j] < gra[s][t])
            {
                temp = dfs(s,t);
                if( m < temp) //从上、下、左、右中找一个最优的
                    m = temp;
            }
        }
    }
    dp[i][j] = m + gra[i][j] ;//本身的值也得加上
    return dp[i][j];
}

int main()
{
    while(cin>>n>>k)
    {
        if(n==-1 && k==-1)
            break;
        int i,j;
        for(i=0;i<n;i++)
            for(j=0;j<n;j++)
            {
                scanf("%d",&gra[i][j]);
                dp[i][j] = 0;
            }
        dfs(0,0);
        cout<<dp[0][0]<<endl;
    }
    return 0;
}

DreamSky 2009-04-03 08:46 发表评论

zju 3049 Diablo II Items

DreamSky — Wed, 01 Apr 2009 04:22:00 GMT

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=2627

/*
3 100
1 102
2 104
100 300
输出：106
*/
#include<iostream>
#include<algorithm>
using namespace std;
struct Node
{
    int normalPrice;
    int magicPrice;
};
int main()
{
    int n,cost;
    while(cin>>n>>cost)
    {
        string str;
        Node infor[1001];
        int i,j,k,a,b,cc;
        char ch;
        for(cc=0,k=0,i=1;i<=n;i++)
        {
            cin>>a;
            ch=getchar();
            if(ch=='\n')//是换行符表示该行只有一个数
            {
                cc += a;
                continue;
            }
            else
            {
                cin>>b;
                if(b-a>cost)//只能是有挣才会卖魔法价
                {
                    infor[k].normalPrice = a;
                    infor[k++].magicPrice = b;
                }
                else
                    cc += a;
            }
        }
        if(k==0 || cc >= cost)
        {
            for(i=0;i<k;i++)
                cc += infor[i].magicPrice;
            cc -= cost*k;
        }
        else
        {
            int temp;
            int remain = cost-cc;//购买鉴定卷轴欠缺的钱
            int *dp = new int[remain+1];//dp[i]表示还差i元钱时损耗的最小价值
            for(i=0;i<=remain;i++)
                dp[i]=100000000;
            dp[0] = 0;
            for(j=0;j<k;j++)//物品
            {
                temp = infor[j].magicPrice - infor[j].normalPrice - cost;//记得减去cost
                //temp表示亏损值
                if(infor[j].normalPrice >= remain && dp[remain] > temp )//正常价格够买鉴定卷轴
                    dp[remain] = temp;
                else
                {
                    for( i=remain - infor[j].normalPrice ; i <= remain ; i++)
                    //i从remain - infor[j].normalPrice
                        if(dp[remain] > dp[i] + temp)
                            dp[remain] = dp[i] + temp;
                    for( i=remain ; i >= infor[j].normalPrice ; i-- )
                        if(dp[i] > dp[i-infor[j].normalPrice] + temp)
                            dp[i] = dp[i-infor[j].normalPrice] + temp;
                }
            }
            if(dp[remain] == 100000000)//无法以魔法价格出售
            {
                for(i=0;i<k;i++)
                    cc += infor[i].normalPrice;
            }
            else
            {
                for(i=0;i<k;i++)
                    cc += infor[i].magicPrice;
                cc -= cost*k + dp[remain];//记得要扣掉买鉴定卷轴和损耗的钱
            }
        }
        cout<<cc<<endl;
    }
    return 0;
}

DreamSky 2009-04-01 12:22 发表评论

zju 2975 Kinds of Fuwas

DreamSky — Tue, 31 Mar 2009 10:49:00 GMT

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1974

//时间复杂度不好
/*
N  B  N  B  J  B
N  B  N  B  J  B
dp[0][N]=1; dp[1][N]=1; dp[2][N]=2;
*/
#include<iostream>
using namespace std;
int dp[251][5];//列号//B J H Y N
int map[251][251];
int fun(char ch)
{
    switch(ch)
    {
    case 'B':return 0;
    case 'J':return 1;
    case 'H':return 2;
    case 'Y':return 3;
    case 'N':return 4;
    }
    return 0;
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int m,n,i,j,k,t;
        char ch;
        long long sum=0;
        cin>>m>>n;
        for(i=0;i<m;i++)
        {
            for(j=0;j<n;j++)
            {
                cin>>ch;
                map[i][j]=fun(ch);
            }
        }
        for(k=0;k<n;k++)
               for(t=0;t<5;t++)
                  dp[k][t]=0;
        for(i=0;i<m-1;i++)//行
        {
            for(j=i+1;j<m;j++)//行
            {
                if(map[i][0]==map[j][0]) //第0列的 i行与j行的福娃一样
                    dp[0][map[j][0]]=1;
                for(k=1;k<n;k++)//列
                {
                    for(t=0;t<5;t++)//获取前面一列的福娃数
                        dp[k][t]=dp[k-1][t];
                    if(map[i][k]==map[j][k])
                    {
                        if(dp[k][map[i][k]]>0)//第一对出现的不能算数，因为只有一列不是矩形
                            sum += dp[k][map[i][k]];
                        dp[k][map[i][k]]++;
                    }
                }
                for(k=0;k<n;k++)
                    for(t=0;t<5;t++)
                         dp[k][t]=0;
            }
        }
        printf("%lld\n",sum);
    }
    return 0;
}

DreamSky 2009-03-31 18:49 发表评论

El Dorado

DreamSky — Tue, 31 Mar 2009 07:57:00 GMT

http://acm.zjgsu.edu.cn/JudgeOnline/showproblem?problem_id=1260

/*
动态规划 dp[i][j]表示从1到i中递增序列长度为j且最后一个数为num[i]的序列个数
*/
#include<iostream>
using namespace std;
__int64 dp[101][101];
int main()
{
    int n,k;
    while(cin>>n>>k && n+k)
    {
        int i,j,t,num[101];
        for(i=1;i<=n;i++)
            cin>>num[i];
        for(i=0;i<=n;i++)
            for(j=0;j<=n;j++)
                dp[i][j]=0;
        for(i=1;i<=n;i++)
            dp[i][1]=1;
        for(j=2;j<=k;j++)
            dp[1][j]=0;
        for(i=2;i<=n;i++)
        {
            for(j=2;j<=k;j++)
            {
                for(t=1;t<i;t++)
                {
                    if(num[i]>num[t])
                        dp[i][j]+=dp[t][j-1];
                }
            }
        }
        __int64 sum=0;
        for(i=k;i<=n;i++)
            sum+=dp[i][k];
        printf("%I64d\n",sum);
    }
    return 0;
}

DreamSky 2009-03-31 15:57 发表评论