UVA 11645 Bits

       这题写的真是无比蛋疼啊,刚开始看,感觉像是数位dp,想了想,应该能做,那就搞吧,费了半天劲终于写完了,数据看起来也没啥问题了,交上去发现wa了,后来发现数据非常大的话就会超出long long的范围,这意思是要写高精啊,当时就郁闷了,没办法,又写了半天大整数加法和乘法,为了把之前的代码改成大整数的又弄了一会儿,写完真是累死人了。。。唉,这代码写的真是不忍直视,凑合看吧……

        对于一个数,先将其转换为二进制。每个数可以用若干个模板和前缀表示,用f(m,b,d)表示这个数的结果,则m表示它的前缀中有多少个Adjacent Bits,b表示它前一位是1还是0,d表示这个数有多少个*,其中*可以表示0和1的任意一个数,比如110就可以用**、10*、110表示,分别对应f(0,0,2)、f(0,0,1)、f(1,0,0),这样就可以写出递推公式:

                 f(m,b,d)=f(m+(b==1?1:0),1,d-1)+f(m,0,d-1);

       用dp[m][b][d]表示状态,将搜过的状态直接返回。

 

附上代码:

 

#include <iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
struct BigInt             //自己写的大整数,比较挫,可以无视
{
    int num[50];
    int len;
    void clear()
    {
        for(int i=0;i<50;++i) num[i]=0;
        len=1;
    }
    void getBigInt(ll n)
    {
        clear();
        len=0;
        while(n!=0)
        {
            num[len++]=n%10;
            n/=10;
        }
        if(len==0) len=1;
    }
    void getBigInt(int n)
    {
        ll m=(ll)n;
        getBigInt(m);
    }
    BigInt operator + (const BigInt a)const
    {
        int n=max(len,a.len);
        BigInt c;
        c.clear();
        c.len=0;
        int jinwei=0;
        int tp;
        for(int i=0;i<n;++i)
        {
            tp=a.num[i]+num[i]+jinwei;
            c.num[c.len++]=tp%10;
            jinwei=tp/10;
        }
        while(jinwei!=0)
        {
            c.num[c.len++]=jinwei%10;
            jinwei/=10;
        }
        return c;
    }
    BigInt operator * (const BigInt a) const
    {
        int m=a.len;
        int n=len;
        BigInt c;
        c.clear();
        c.len=0;
        int bs=0,w=0,jinwei=0,tp;
        for(int i=0;i<m;++i)
        {
            jinwei=0;
            w=0;
            for(int j=0;j<n;++j)
            {
                tp=a.num[i]*num[j]+jinwei+c.num[bs+w];
                c.num[bs+w]=tp%10;
                jinwei=tp/10;
                w++;
            }
            while(jinwei!=0)
            {
                c.num[bs+w]=jinwei%10;
                jinwei/=10;
                w++;
            }
            c.len=max(c.len,bs+w);
            bs++;
        }
        while(c.num[c.len-1]==0&&c.len>1)
           c.len--;
        return c;
    }
};
BigInt dp[70][2][70];
BigInt pp[70];
int nu[100];
BigInt f(int m,int b,int d)
{
    BigInt tp,td;
    tp.clear();
    td.clear();
    if(d<0) return tp;
    tp.getBigInt(m);
    if(d==0) return tp;
    if(d==1)
    {
        td.getBigInt(b==1?1:0);
        return td+pp[d]*tp;
    }
    if(d==2)
    {
        td.getBigInt(b==1?3:1);
        return td+pp[d]*tp;
    }
    if(dp[m][b][d].len>0) return dp[m][b][d];
    dp[m][b][d].clear();
    dp[m][b][d]=dp[m][b][d]+f(m+(b==1?1:0),1,d-1);
    dp[m][b][d]=dp[m][b][d]+f(m,0,d-1);
    return dp[m][b][d];
}
BigInt slove(int n)
{
    int m=0,b=0,d=n;
    BigInt res;
    res.clear();
    //res+=f(0,0,n-1);
    for(int i=0;i<n;++i)
    {
        d--;
        if(nu[i]==1)
        {
            res=res+f(m,0,d);
            m+=(b==1)?1:0;
            b=1;
        }
        else
        {
            b=0;
        }
    }
    BigInt tp;
    tp.getBigInt(m);
    res=res+tp;
    return res;
}
void Init()
{
    ll z=1;
    pp[0].getBigInt(z);
    for(int i=1;i<=63;++i)
    {
        z=z<<1;
        pp[i].getBigInt(z);
    }
    for(int i=0;i<70;++i)
    {
        for(int j=0;j<2;++j)
        {
            for(int k=0;k<70;++k)
            {
                dp[i][j][k].len=-1;
            }
        }
    }
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    Init();
    int tcase=0;
    ll x;
    while(cin>>x)
    {
        if(x<0) break;
        tcase++;
        int len=0;
        while(x!=0)
        {
            nu[len++]=x%2;
            x=x>>1;
        }
        int temp;
        for(int i=0;i<len/2;++i)
        {
            temp=nu[i];
            nu[i]=nu[len-i-1];
            nu[len-i-1]=temp;
        }
        BigInt res=slove(len);
        cout<<"Case "<<tcase<<": ";
        for(int i=res.len-1;i>=0;--i)
           cout<<res.num[i];
        cout<<endl;
    }
    return 0;
}


 


 

你可能感兴趣的:(动态规划,ACM,uva)