C++博客-Acmiyou/Sulingr-随笔分类-ACM

有道难题...

米游 — Sun, 31 May 2009 11:46:00 GMT

Problem Statement

一个二进制序列由下面的伪代码生成：

string A = "0"
            While (A的长度小于等于n)
            创建一个和A一样长度的字符串B
            For i=0,1,...length(A)-1
            If (i 是完全平方数)
            B[i] = 1-A[i]
            Else
            B[i] = A[i]
            令A = A + B (即将B拼接在A后面)
            End While
            Return A

请注意，在上面的伪代码中，A[i]和B[i]分别表示字符串A和B中下标为i的字符（下标编号从0开始）。对“完全平方数”的定义是，对于整数i，存在整数j，使得i= j *j，则称i为完全平方数。

下面具体说明序列生成的过程：如果n=7，则在每一轮迭代中，A的取值依次为：0, 01, 0110, 01101010，所以最终产生的二进制序列就是0,1,1,0,1,0,1,0

请返回上述序列中下标为n的数字（该序列的下标从0开始）(0=

Definition

Class:	BinarySequence
Method:	getValue
Parameters:	int
Returns:	int
Method signature:	int getValue(int n)
(be sure your method is public)

//题目看错了...i是完全平方数..郁闷咯..以下代码是错的哈.

#include<iostream>
#include<string>
using namespace std;
class BinarySequence
{
public:
    int getValue(int n)
    {
        int l,i=0,k;
            //n由 n-2^(log2(n))变换而来
        for(l=n;l>0;i^=1)
        {
            for(k=1;k<l;k<<=1);
            l-=(k==l?k:(k>>1));
        }
        return i;
    }
};  //这样写对嘛？刚开始写错代码才得了一半分数..

米游 2009-05-31 19:46 发表评论

zoj 3211 Dream City

米游 — Mon, 25 May 2009 07:35:00 GMT

The 6th Zhejiang Provincial Collegiate Programming Contest 的J题比赛的时候少考虑了一个排序..DP状态都出来了居然还是少考虑一个情况..

6题与7题的差距是一个天一个地哎....

dp[i][j] 表示前i颗树第j天可达到的最大值

1#include<iostream>
2#include<algorithm>
3using namespace std;
4struct node
5{
6    int a,b;
7};
8bool operator<(const node& d,const node& e)
9{
10    return d.b<e.b;
11}
12node r[300];
13int dp[300][300];
14int main()
15{
16    int n,m,t;
17    cin>>t;
18    while(t--)
19    {
20        cin>>n>>m;
21        for(int i=1;i<=n;i++)
22            cin>>r[i].a;
23        for(int i=1;i<=n;i++)
24            cin>>r[i].b;
25        sort(r+1,r+n+1);
26        memset(dp,0,sizeof(dp));
27        //dp[i][j] 表示前i个第j天的最大值
28        for(int i=1;i<=n;i++)
29            for(int j=1;j<=m;j++)
30            {
31                if(j>i)
32                    break;
33                dp[i][j]=max(dp[i-1][j],dp[i-1][j-1]+r[i].a+r[i].b*(j-1));
34            }
35        cout<<dp[n][m]<<endl;
36    }
37}

米游 2009-05-25 15:35 发表评论

09.5.23 退役感言

米游 — Sun, 24 May 2009 11:20:00 GMT

    在左手西瓜右手泡面下，思索着如何展开这段感言。其实我真的不知道自己想说些什么.兴许对于那些并不熟悉ACM的朋友们，也许他们仅就认为这是一项竞赛而已。而对于这群深入并了解它的我们而已，并不只是表面上的那么简单。这正如你深入而了解它的同时，它也在同一时刻改变你，也渗透到你生活里。     初识程序设计是在大一的C语言。而初识ACM却是在大一下。这对于那些时候爱玩游戏的我是种既陌生又冲动的感觉。大一下那段时间我并没有学到多少真正的ACM知识。也仅仅只是初识。称之为New Acmer.而后经过院赛选拔进了08年的校集训队。其实当时我们那个队伍真的没多少人。而热爱AC并坚持下来的更没多少人。几个老师的上课也是匆匆了事。其实我觉得我自己在当时根本没学到多少的知识也曾一度想过放弃。当时认识了陆老大老张老陈。而我跟老陈后来都是一同激励。其实我当初并不认识老李和老邱。当时的集训队里只能算的上是各自做各自的。缺少很好的交流。而且老师给予的帮助也是相当少的。全是自己啃题AC题不懂的只能在网上搜罗各种资料。再AC了150道左右后发觉一直停滞不前了。陆老大和老张说过这是个坎。你必须了解更多的算法知识和坚持做更多的题目才能熬过去。而后熬过那段时间也到了省赛组队我老李老邱分到了一组。其实我们当时也是懵懵懂懂的。什么都不懂却什么都懂一点的那种感觉。只是没想到意外的摘得了铜牌那是我们对自己那段时间的一个肯定，只是遗憾的是陆老大和老张在那次比赛后退役了老邱也因考研和退役了。最后留下了我老李老陈老黄。
    于是教练也换了人。我们这几个老队员担当起了小教练的职责。选拔新队员给集训队讲解各种算法。也许对于我们也是个锻炼但那段时间里真的非常辛苦。老陈一直抱怨。后来的队员里有许多生面孔了集训队也热闹了，可是热爱ACM和努力并不多。而新集训队的模式一直没能得到很好的改观。我和老李以及池MM 肖剑是大三的其他的都是大二的。经过一次次比赛和磨练最终留下也是仅有的15个人。也是09年省赛的队伍。09年4月底去宁波参加邀请赛结果惨败而归对我们是一大打击。省赛里也只是收获了三个铜。小菜的队伍 42名。我和老李及池MM的队伍55名。老陈的队伍65名。
     其实我们并没有失败。有着许多的像我们一样的热爱ACM的朋友们。大家都在不停的AC中。AC是快乐的。我们始终相信AC带给我们的不仅仅只是一种竞赛。而是一种不同于他人的生活。
     我们收获友谊收获知识更收获了许多别处收获不到得能力。虽然某位大牛说Programming is not all your life。但我还是相信ACM依然是种不一样的生活。
     ACM是快乐的是痛苦的是百感交集的而我们却在当中成长。

仅以此文
铭记AC的两个年华
                                                                   Acmiyou 落笔

米游 2009-05-24 19:20 发表评论

RMQ ST算法 (区间最大(最小)值问题)

米游 — Mon, 18 May 2009 06:50:00 GMT

#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
#define M 100010
#define MAXN 500
#define MAXM 500
int dp[M][18];
int dp2[MAXN][MAXM][10][10]
/*
*一维RMQ ST算法
*构造RMQ数组 makermq(int n,int b[]) O(nlog(n))的算法复杂度
*dp[i]j] 表示从i到i+2^j -1中最大的一个值
*dp[i][j]=max{dp[i][j-1],dp[i+2^(j-1)][j-1]}
*查询RMQ rmq(int s,int v)
*将s ->v 分成两个2^k的区间
*即 k=(int)log2(s-v+1)
*查询结果应该为 max(dp[s][k],dp[v-2^k+1][k])
*/
int rmq(int s,int v)
{
    int k=(int)(log((v-s+1)*1.0)/log(2.0));
    return max(dp[s][k],dp[v-(1<<k)+1][k]);
}
void makermq(int n,int b[])
{
    int i,j;
    for(i=1;i<=n;i++)
        dp[i][0]=b[i];
    for(j=1;(1<<j)<=n;j++)
        for(i=1;i+(1<<j)-1<=n;i++)
            dp[i][j]=max(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
}

/*
*二维RMQ ST算法
*构造RMQ数组 makermq(int n,int m,int b[][]) O(n*m*log(n)*log(m))算法复杂度
*dp2[row][col][i][j] 表示行从row ->row +2^i-1 列从col ->col +2^j-1 二维区间里最大值
*dp2[row][col][i][j] = 下行
*max{dp2[row][col][i][j-1],dp2[row][col][i-1][j],dp2[row][col+2^(j-1)][i][j-1],dp2[row+2^(i-1)][col][i-1][j]}
*查询RMQ rmq(int sx,int ex,int sy,int ey)
*同一维的将sx->ex 分为两个2^kx区间将 sy->ey分为两个2^ky的区间
*kx=(int)log2(ex-sx+1) ky=(int)log2(ey-sy+1)
*查询结果为
*max{dp2[sx][sy][kx][ky],dp2[sx][ey-2^ky+1][kx][ky],dp2[ex-2^kx+1][sy][kx][ky],dp2[ex-2^kx+1][ey-2^ky+1][kx][ky]}
*/

void makermq(int n,int m,int b[][MAXM])
{
    int row,col,i,j;
    for(row=1;row<=n;row++)
        for(col=1;col<=m;col++)
            dp2[row][col][0][0]=b[row][col];
    for(i=0;(1<<i)<=n;i++)
        for(j=0;(1<<j)<=m;j++)
        {
            if(i==0&&j==0) continue;
            for(row=1;row+(1<<i)-1<=n;row++)
                for(col=1;col+(1<<j)-1<=m;col++)
                {
                if(i==0)
                    dp2[row][col][i][j]=max(dp2[row][col][i][j-1],dp2[row][col+(1<<(j-1))][i][j-1]);
                else
                    dp2[row][col][i][j]=max(dp2[row][col][i-1][j],dp2[row+(1<<(i-1))][col][i-1][j]);
                }
        }
}
int rmq(int sx,int ex,int sy,int ey)
{
int kx=(int)(log((ex-sx+1)*1.0)/log(2.0)),ky=(int)(log((ey-sy+1)*1.0)/log(2.0));
return max(max(dp2[sx][sy][kx][ky],dp2[sx][ey-(1<<ky)+1][kx][ky]),max(dp2[ex-(1<<kx)+1][sy][kx][ky],dp2[ex-(1<<kx)+1][ey-(1<<ky)+1][kx][ky]));
}

米游 2009-05-18 14:50 发表评论

使用后缀数组解决zoj 3199 Longest Repeated Substring

米游 — Sat, 09 May 2009 03:22:00 GMT

摘要: /*初看这题以为是传统意义上的最长重复子串.其实不然,看例子就明白*/接触这题后才开始看Suffix_array的资料.一篇论文,里面谈到如何使用O(nlogn)的方法构造后缀数组SA.并且用0(nlongn)的方法构造height数组. 点击下载后缀数组论文以下代码写的有点粗糙..排序上其实可以优化很多.我只使用sort()进行排序惭愧..... #include阅读全文

米游 2009-05-09 11:22 发表评论

线段树求矩形覆盖的周长 pku 1177

米游 — Sun, 03 May 2009 13:02:00 GMT

摘要: /**//*Author:miyouPROG:the perimeter of many rectanglesLANG:C++*/#include#include#includeusing namespace std;#define M... 阅读全文

米游 2009-05-03 21:02 发表评论

hdu 2816 即老菜鸟杯的1008题目

米游 — Sat, 02 May 2009 12:49:00 GMT

这次比赛里最简单的一道题目了...
直接模拟5次解密过程就可以过的水题..一道出自百度贴吧的真实故事的题目哈哈..

#include<iostream>
#include<cstring>
using namespace std;
char num[10][5]={"","_","ABC","DEF","GHI","JKL","MNO","PQRS","TUV","WXYZ"};
char zh1[27]={"QWERTYUIOPASDFGHJKLZXCVBNM"};
char zh[27];
char str[1005];
char c[505];
int main()
{
    for(int i=0;i<26;i++)
        zh[zh1[i]-'A']='A'+i;

    while(gets(str))
    {
        int len=strlen(str);
        int k=0;
        for(int i=0;i<len;i+=2)
            c[k++]=num[str[i]-'0'][str[i+1]-'0'-1];
        len=k;
        c[len]='\0';
        for(int i=0;i<len;i++)
            c[i]=zh[c[i]-'A'];
        int l1;
        if(len%2==0)
            l1=len/2;
        else
            l1=len/2+1;
        int i,j;
        char res[505];
        k=0;
        for(i=0,j=l1;i<l1&&j<len;i++,j++)
        {
            res[k++]=c[i];
            res[k++]=c[j];
        }
        for(;i<l1;i++)
            res[k++]=c[i];
        for(;j<len;j++)
            res[k++]=c[j];
        for(i=k-1;i>=0;i--)
            putchar(res[i]);
        putchar('\n');
    }
    return 0;
}

米游 2009-05-02 20:49 发表评论

hdu 2813 即老菜鸟杯 1005题

米游 — Sat, 02 May 2009 12:47:00 GMT

二分图的最小权值匹配问题。。。直接用浙大模板水过..

#include<iostream>
#include<string>
#include<map>
#define MAXN 310
#define inf 1000000000
#define _clr(x) memset(x,0xff,sizeof(int)*n)
using namespace std;
int kuhn_munkras(int m,int n,int mat[][MAXN],int* match1,int* match2){
    int s[MAXN],t[MAXN],l1[MAXN],l2[MAXN],p,q,ret=0,i,j,k;
    for (i=0;i<m;i++)
        for (l1[i]=-inf,j=0;j<n;j++)
            l1[i]=mat[i][j]>l1[i]?mat[i][j]:l1[i];
    for (i=0;i<n;l2[i++]=0);
    for (_clr(match1),_clr(match2),i=0;i<m;i++){
        for (_clr(t),s[p=q=0]=i;p<=q&&match1[i]<0;p++)
            for (k=s[p],j=0;j<n&&match1[i]<0;j++)
                if (l1[k]+l2[j]==mat[k][j]&&t[j]<0){
                    s[++q]=match2[j],t[j]=k;
                    if (s[q]<0)
                        for (p=j;p>=0;j=p)
                            match2[j]=k=t[j],p=match1[k],match1[k]=j;
                }
        if (match1[i]<0){
            for (i--,p=inf,k=0;k<=q;k++)
                for (j=0;j<n;j++)
                    if (t[j]<0&&l1[s[k]]+l2[j]-mat[s[k]][j]<p)
                        p=l1[s[k]]+l2[j]-mat[s[k]][j];
            for (j=0;j<n;l2[j]+=t[j]<0?0:p,j++);
            for (k=0;k<=q;l1[s[k++]]-=p);
        }
    }
    for (i=0;i<m;i++)
        ret+=mat[i][match1[i]];
    return ret;
}

int main()
{
    int m,n,k,v;
    char name1[20],name2[20];
    map<string,int>::iterator it;
    int mat[MAXN][MAXN];
    map<string,int> name,mapname;
    while(cin>>m>>n>>k)
    {
        int j=0,l1,l2,l=0;
        for(int i=0;i<m;i++)
            for(int j=0;j<n;j++)
                mat[i][j]=-10000000;
        name.clear();
        for(int i=0;i<k;i++)
        {
            cin>>name1>>name2>>v;
            string name3(name1),name4(name2);
            it=name.find(name3);
            if(it==name.end())
            {
                name[name3]=j;
                l1=j++;
            }
            else
            {
                l1=it->second;
            }
            it=mapname.find(name4);
            if(it==mapname.end())
            {
                mapname[name4]=l;
                l2=l++;
            }
            else
            {
                l2=it->second;
            }
            mat[l1][l2]=v*(-1);
        }
        int match1[MAXN],match2[MAXN];
        cout<<-1*kuhn_munkras(m,n,mat,match1,match2)<<endl;
    }
    return 0;
}

米游 2009-05-02 20:47 发表评论

hdu 2812 即老菜鸟杯 1004

米游 — Sat, 02 May 2009 12:45:00 GMT

队友写的代码..没看题目。
这是hdu大牛的解题报告:http://acm.hdu.edu.cn/forum/read.php?tid=13133

#include<iostream>
using namespace std;
int h[15];
bool f[15];
double h1[15];
int ji_suan(double t1,double t2);
int main()
{
    int n,flag,i,a,t;
    while(scanf("%d",&n)!=EOF)
    {
        flag=0;
        //memset(f,false,sizeof(f));
        for(i=0;i<14;i++)
            f[i]=false;
        for(i=0;i<n;i++)
        {
            scanf("%d%lf",&a,&h1[i]);
            h[a]=i;
            f[a]=true;
            if(a==14)
            {
                if(h1[i]>5.0||h1[i]<5.0)
                    flag=1;
            }
            if(a==7)
            {
                if(h1[i]<0||h1[i]>0.8)
                    flag=1;
            }
        }
        if(flag==1)
        {
            printf("No\n");
            continue;
        }
        printf("Yes\n");
        if(f[1]==true&&f[2]==true)
        {
            t=ji_suan(h1[h[1]],h1[h[2]]);
            if(t==0)
                printf("3 L\n");
            else
                printf("3 H\n");
            flag=1;

        }
        if(f[4]==true&&f[5]==true)
        {
            t=ji_suan(h1[h[4]],h1[h[5]]);
            if(t==0)
                printf("6 L\n");
            else
                printf("6 H\n");
            flag=1;

        }
        if(f[9]==true&&f[10]==true)
        {
            t=ji_suan(h1[h[9]],h1[h[10]]);
            if(t==0)
                printf("8 L\n");
            else
                printf("8 H\n");
            flag=1;

        }
        if(f[12]==true&&f[13]==true)
        {
            t=ji_suan(h1[h[12]],h1[h[13]]);
            if(t==0)
                printf("11 L\n");
            else
                printf("11 H\n");
            flag=1;
        }
    }
    return 0;
}
int ji_suan(double t1,double t2)
{
    if(t1>=4.5&&t1<=5.5&&t2>=4.5&&t2<=5.5)
        return 0;
    if(t1>=4.5&&t1<=5.5&&(t2<4.5||t2>5.5))
        return 1;
    if(t2>=4.5&&t2<=5.5&&(t1<4.5||t1>5.5))
        return 1;
    else
        return 1;
}

米游 2009-05-02 20:45 发表评论

hdu 2811 即老菜鸟杯 1003

米游 — Sat, 02 May 2009 12:42:00 GMT

比赛最后才和同学一起AC的一道题目...
f(n)=sum(C(n-1,j)*f(j)*f(n-1-j),0<=j<=(n-1)/2) 当然 n-1=2*j的时候组合数必须除以2.
这个时候还得考虑是否会越界..所以有个%k的处理

#include<stdio.h>
__int64 mmg[1005],n,k;
__int64 C[1005][1000],flag[1005][1000];
int main()
{
    int i,j;
    while(scanf("%I64d%I64d",&n,&k)!=EOF)
    {
        C[1][0]=C[1][1]=1%k;
        flag[1][0]=flag[1][1]=(1/k)%2;
        for(i=2;i<=n;i++)
        {
            C[i][0]=C[i][i]=1%k;
            flag[i][0]=flag[i][i]=(1/k)%2;
            for(j=1;j<i;j++)
            {
                C[i][j]=(C[i-1][j]+C[i-1][j-1])%k;
                flag[i][j]=flag[i-1][j]+flag[i-1][j-1]+(C[i-1][j]+C[i-1][j-1])/k;
                flag[i][j]=flag[i][j]%2;
            }

        }
        mmg[1]=1%k;
        mmg[2]=1%k;
        for(i=3;i<=n;i++)
        {
            mmg[i]=mmg[i-1];
            for(j=1;j<=(i-1)/2;j++)
            {
                if(i-1==2*j)
                {
                    if(flag[i-1][j]==1)
                        mmg[i]+=(((mmg[j]%k)*(mmg[i-j-1]%k))%k)*(((C[i-1][j]+k)/2)%k);
                    else
                        mmg[i]+=(((mmg[j]%k)*(mmg[i-j-1]%k))%k)*(((C[i-1][j])/2)%k);
                }
                else
                    mmg[i]+=(((mmg[j]%k)*(mmg[i-j-1]%k))%k)*(C[i-1][j]%k);
                mmg[i]=mmg[i]%k;
            }
        }
        printf("%I64d\n",mmg[n]);
    }
}

米游 2009-05-02 20:42 发表评论

hdu 2810 即老菜鸟杯1002题解题报告

米游 — Sat, 02 May 2009 12:37:00 GMT

摘要: 刚拿到题目，就想到将集合看成一个独立点，求次MST。再对每个集合内的所有点求MST..可惜比赛的时候没过这题。。错误原来在空的集合是不需要连接(这个没考虑所以出错了)我用了prim算法没什么优化..903ms过的.(可以用堆优化下) #include#includeusing namespace std;cons... 阅读全文

米游 2009-05-02 20:37 发表评论

参数搜索(求最优解问题)

米游 — Sun, 12 Apr 2009 04:41:00 GMT

:题目例子：http://acm.zjgsu.edu.cn/JudgeOnline/showproblem?problem_id=1277
题目的大意是：
有N个房子在一条街上，（街是条直线）每个房子有个离这个街的最末端得距离。现在要在这条街上安装m个路由器使得每个房子都能有个距离最近的路由器。并且使得所有房子和它距离最近的路由器之间的距离最大一个是最小的。
初一想就是一个经典DP问题..
比如考虑两个房子一个路由器.肯定是将路由器放在中间位置..
于是对房子的距离排个序就可以很快得出一个DP的表达式。
F[i,k]表示前i个房子放k个路由器的最优解..
则F[i,k]=min{max(F[j,k-1],(house[i]-house[j+1])/2),1<=j表示前j个房子放k-1路由器，第K个路由器放在第j+1和第i个房子的中点..
算法是O(m*n^2) 对于数据两n<=100000 无论如何都优化不了的DP 已经难以解决这个问题了..
后来经同实验室的朋友介绍了一个算法参数搜索的应用..
并找到一个论文.http://acm.zjgsu.edu.cn/Software/canshusousuo.doc(汪汀著) 转载至中华信息网

关于怎么解这个问题就不多说了..就是找最优解.首先找一个最小的解，最大解，然后二分判断这个某个解是否是可行解。

代码如下：

#include<iostream>
#include<algorithm>
using namespace std;
int house[100002];
int n,m;
bool isok(int value)
{
    int i=1,j,p=1;
    for(j=1;j<=n;j++)
    {
                if(i>m) return false;
        if((house[j]-house[p])*10>value*2)
        {
            p=j;
            i++;
        }
    }
    if(i<=m)
        return true;
    return false;
}
void AC()
{
    int ans;
    int low=0,hight=(house[n]-house[0])*10,mid;
    while(low<=hight)
    {
        mid=(low+hight)/2;
        if(isok(mid))
            hight=mid-1;
        else
            low=mid+1;
    }
    if(isok(low))
        ans=low;
    else
        ans=hight;
    printf("%.1lf\n",ans*1.0/10);
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&m,&n);
        for(int i=1;i<=n;i++)
            scanf("%d",&house[i]);
        house[0]=0;
        sort(house,house+n+1);
        AC();
    }
}

米游 2009-04-12 12:41 发表评论

pku 2353 DP(五)

米游 — Fri, 03 Apr 2009 09:55:00 GMT

题目的意思是..第一行走到第M行的最小消费. 可以往下走一步，可以往左走一步,可以往右走一步.
可以从第一行任意个位置出发。只要到达第M行的任意个位置就结束。。
按意思可得出一个简单dp(i,j) = Min{dp(i-1,j),dp(i,j-1),dp(i,j+1)}+v[i][j].但是你会发觉DP的时候似乎
dp(i,j+1)是在dp(i,j)之后求的..故而必须得双向DP..

现在考虑第一行的任意一个列都不会去往左右方向走的..因为它如果往左右走的话,则可以直接选择从左边或者右边开始就行.
故可以初始化dp数组
for(int i=1;i<=m;i++)
dp[1][j]=v[1][j];

在考虑dp[i,j]时候.可以这么考虑.
先比较dp[i-1,j]和dp[i-1,j]。通过这个可先求得dp[i][j+1],然后在做一次的dp[i][j+1]的比较..

代码如下：

#include<iostream>
using namespace std;
int m,n;
int dp[105][505],v[105][505],flag[105][505];
void print(int i,int j)
{
    if(i==1)
    {
        printf("%d\n",j);
        return ;
    }
    if(flag[i][j]==1)
        print(i-1,j);
    if(flag[i][j]==2)
        print(i,j-1);
    if(flag[i][j]==3)
        print(i,j+1);
    printf("%d\n",j);
}
int main()
{
    scanf("%d%d",&m,&n);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            scanf("%d",&v[i][j]);
    for(int j=1;j<=n;j++)
        dp[1][j]=v[1][j];
    for(int i=2;i<=m;i++)
    {
        dp[i][1]=dp[i-1][1]+v[i][1];//表示每行第一个房间暂时只能从上边走下来
        flag[i][1]=1; //flag标记等于1表示从上边走下
        for(int j=2;j<=n;j++) //求从上往下走和从左往右走的最小值
        {
            dp[i][j]=v[i][j];
            if(dp[i-1][j]<dp[i][j-1])
            {
                dp[i][j]+=dp[i-1][j];
                flag[i][j]=1;//flag标记等于1表示从上边走到当前位置
            }
            else
            {
                dp[i][j]+=dp[i][j-1];
                flag[i][j]=2;//flag标记等于2表示从左边走到当前位置
            }
        }
        for(int j=n-1;j>=1;j--) //再比较从右往左走的与之前的比较,取更小的.
            if(dp[i][j+1]+v[i][j]<dp[i][j])
            {
                dp[i][j]=dp[i][j+1]+v[i][j];
                flag[i][j]=3;//flag标记等于3表示从右边走到当前位置
            }
    }
    int Min=0x7fffffff,my;
    for(int i=1;i<=n;i++)
        if(Min>dp[m][i])
        {
            Min=dp[m][i];my=i;
        }
    print(m,my);
    return 0;
}

米游 2009-04-03 17:55 发表评论

pku 1631 DP(四)+二分查找

米游 — Wed, 01 Apr 2009 05:14:00 GMT

简单的概括下这题的意思就是求一个最长上升子序列.由于number的数很多.用普通的n^2的DP求最长上升子序列肯定会超时。。故采用二分查找的方法。

#include<iostream>
using namespace std;
int num[40001],dp[40001],n;
int DP()
{
    int len=1,left,right,mid;
    dp[0]=num[0];
    for(int i=1;i<n;i++)
    {
        right=len-1,left=0;
        while(left<=right)
        {
            mid=(left+right)/2;
            if(num[i]<=dp[mid])
                right=mid-1;
            else
                left=mid+1;
        }
        dp[left]=num[i];
        if(left==len) len++;
    }
    cout<<len<<endl;
    return len;
}
int main()
{
    int T;
    cin>>T;
    while(T--&&cin>>n)
    {
        for(int i=0;i<n;i++)
            cin>>num[i];
        int len=DP();
    }
    return 0;
}

米游 2009-04-01 13:14 发表评论

pku 1191 棋盘分割 (DP)(三)

米游 — Wed, 01 Apr 2009 04:15:00 GMT

摘要: 题目中文不解释了。刚开始做这题，我想到的就是递归.因为子问题与大问题完全一样..于是连公式都未化简，写了个递归程序。测试的时候发现n越大,我的程序就很慢才得出结果了..递归的函数是这么写的double getMinr(int mini,int maxi,int minj,int maxj,int k) // 表示这个矩形块分割成k块的最小均方差. double getMinr(int&... 阅读全文

米游 2009-04-01 12:15 发表评论

pku 1609解题报告(DP)(二)

米游 — Tue, 31 Mar 2009 07:58:00 GMT

题目意思：
对于给定的矩形块,宽a,长b，而对于其他的矩形块，如果宽度小于等于a,长度小于等于b,它可以被(a,b)的矩形块覆盖.题目意思就是给n个矩形块,并且已知每个矩形块的宽度和长度..让你求覆盖之后最大的块数..
比如有3块矩形块
1 1
2 3
3 2
则 (1,1)可被(2,3)覆盖 (1,1)也可被(3,2)覆盖而(2,3)不能被(3,2) 所以覆盖之后最多只能有2块。。

所以(a1,b1)<=(a2,b2)的情况下,可以被覆盖，故这题应该变成一个二维的最大上升子序列。

故可以考虑对宽度a进行从小到大排列之后.可对长度b求最大上升子序列.
最大上升子序列的求法就是.
考虑b[i]当前这个数对于i之前的数b[j](0<=j则dp[i]=max(dp[j]+1){(b[j]
代码如下：

#include<iostream>
using namespace std;
int d[10001][2],n,dp[10001];
int cmp(void const *a,void const *b)
{
    int *aa=(int *)a,*bb=(int *)b;
    if(aa[0]!=bb[0])
        return aa[0]-bb[0];
    else
        return aa[1]-bb[1];
}
int main()
{
    while(cin>>n,n)
    {
        for(int i=0;i<n;i++)
            cin>>d[i][0]>>d[i][1];
        qsort(d,n,sizeof(d[0]),cmp);
        for(int i=0;i<n;i++)
            dp[i]=1;
        int max=0;
        for(int i=1;i<n;i++)
        {
            for(int j=0;j<i;j++)
            {
                if(d[j][1]<=d[i][1]&&dp[i]<dp[j]+1)
                    dp[i]=dp[j]+1;
            }
            if(max<dp[i])
                max=dp[i];
        }
        cout<<max<<endl;
    }
    cout<<'*'<<endl;
}

米游 2009-03-31 15:58 发表评论

pku 1179 解题报告（DP）

米游 — Mon, 30 Mar 2009 08:57:00 GMT

题目意思如下:
对于给定多边形，N个顶点，N条边。每条边表示一个‘+’或者‘*’,每个顶点是个数值.
在删除一条边的情况下可进行以下操作：
.选择一条边和该边的两个顶点。通过边上的符号计算合并这两个顶点为一个顶点。直到没有边的可选择时.游戏结束。则最后会得到一个顶点的值
题目的意思让你求出再删除哪一条边的情况，并且通过合理选择边得到一个最大值。

假设 R(i,j)表示从第i个结点开始，顺时针方向连续j个顶点的某个子链,可以考虑该链最后通过中间某条边合成为一个顶点.令 i<=s<=i+j 从R(i,s) 和R(s+1,j)这两个子链通过边(s,s+1)合成为一个顶点
令
m1为R(i,s)这个子链的结果 a为R(i,s)这个子链的最小结果，b为为R(i,s)这个子链的最大结果，则a<=m1<=b.
m2为R(s+1,j)这个子链的结果.c为R(s+1,j)的这个子链的最小结果,d为R(s+1,j)的这个子链的最大结果,则c<=m2<=d
故
如果边(s,s+1)为符号‘+’时候 R(i,j)的最大值应该为(b+d)
如果边(s,s+1)为符号‘*’时候R(i,j)的最大值应该为max(ac,ad,bc,bd)
//需要一正一负的情况..(a,b)=(-8,-5),(c,d)=(5,8)

1#include<iostream>
2#include<algorithm>
3using namespace std;
4int v[52],R[52][52][2],n;
5char ed[52];
6int getMaxr(int vec,int len,int flag) //flag=0 表示求最小值,flag=1表示求最大值
7{
8    int Lmax,Rmax,Lmin,Rmin;
9    if(len==0)
10        return v[vec];
11    if(len==1)
12    {
13        if(ed[vec%n+1]=='t')
14            R[vec][len][flag]=v[vec]+v[vec%n+1];
15        else
16            R[vec][len][flag]=v[vec]*v[vec%n+1];
17        return R[vec][len][flag];
18    }
19    if(R[vec][len][flag]>-32769&&R[vec][len][flag]<32769)
20        return R[vec][len][flag];
21    for(int i=1;i<=len;i++)
22    {
23        Lmax=getMaxr(vec,i-1,1);
24        Rmax=getMaxr((vec+i-1)%n+1,len-i,1);
25        Lmin=getMaxr(vec,i-1,0);
26        Rmin=getMaxr((vec+i-1)%n+1,len-i,0);
27    //    cout<
28        if(ed[((vec+i-1)%n+1)]=='t')
29        {
30            if(flag&&R[vec][len][flag]<Lmax+Rmax)
31                R[vec][len][flag]=Lmax+Rmax;
32            if(!flag&&R[vec][len][flag]>Lmin+Rmin)
33                R[vec][len][flag]=Lmin+Rmin;
34        }
35        else
36        {
37            int temp1=max(max(Lmax*Rmin,Lmax*Rmax),max(Lmin*Rmin,Lmin*Rmax));
38            int temp2=min(min(Lmax*Rmin,Lmax*Rmax),min(Lmin*Rmin,Lmin*Rmax));
39            if(flag&&R[vec][len][flag]<temp1)
40                R[vec][len][flag]=temp1;
41            if(!flag&&R[vec][len][flag]>temp2)
42                R[vec][len][flag]=temp2;
43        }
44        //cout<
45    }
46    return R[vec][len][flag];
47}
48int main()
49{
50    int flag,re[52],Maxx;
51    while(cin>>n)
52    {
53        Maxx=-32769;
54        flag=0;
55        for(int i=1;i<=n;i++)
56        {
57            getchar();
58            cin>>ed[i]>>v[i];
59        }
60        for(int i=0;i<=n;i++)
61            for(int j=0;j<=n;j++)
62            {
63                R[i][j][0]=32769;
64                R[i][j][1]=-32769;
65            }
66        for(int i=1;i<=n;i++)
67        {
68            re[i]=getMaxr(i,n-1,1);
69            if(re[i]>Maxx)
70                Maxx=re[i];
71        }
72        printf("%d\n",Maxx);
73        for(int i=1;i<=n;i++)
74        {
75            if(Maxx==re[i])
76                if(flag==0)
77                    flag++,cout<<i;
78                else
79                    cout<<" "<<i;
80        }
81        cout<<endl;
82    }
83}

米游 2009-03-30 16:57 发表评论