BZOJ 3107 [cqoi2013]二进制a+b (DP)

3107: [cqoi2013]二进制a+b

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 995 Solved: 444
[Submit][Status][Discuss]

Description

输入三个整数a, b, c，把它们写成无前导0的二进制整数。比如a=7, b=6, c=9，写成二进制为a=111, b=110, c=1001。接下来以位数最多的为基准，其他整数在前面添加前导0，使得a, b, c拥有相同的位数。比如在刚才的例子中，添加完前导0后为a=0111, b=0110, c=1001。最后，把a, b, c的各位进行重排，得到a’, b’, c’，使得a’+b’=c’。比如在刚才的例子中，可以这样重排：a’=0111, b’=0011, c’=1010。

你的任务是让c’最小。如果无解，输出-1。

Input

输入仅一行，包含三个整数a, b, c。

Output

输出仅一行，为c’的最小值。

Sample Input

7 6 9

Sample Output

HINT

a,b,c<=2^30

题解：

参考代码：dp[i][j][k][l][m]表示：第I位，第1,2,3个数分别用了几个一，m表示是否有进位：

分类讨论：

　　对于m=0的情况:　

　　 dp[i+1][j+1][k+1][l+1][1]=min(dp[i+1][j+1][k+1][l+1][1],tmp+(1<<i+1));

    dp[i+1][j+1][k][l+1][0]=min(dp[i+1][j+1][k][l+1][0],tmp+(1<<i));

    dp[i+1][j][k+1][l+1][0]=min(dp[i+1][j][k+1][l+1][0],tmp+(1<<i));

    dp[i+1][j][k][l][0]=min(dp[i+1][j][k][l][0],tmp);
　对于m=1的情况：

    dp[i+1][j+1][k+1][l+1][1]=min(dp[i+1][j+1][k+1][l+1][1],tmp+(1<<i+1));

    dp[i+1][j][k+1][l][1]=min(dp[i+1][j][k+1][l][1],tmp+(1<<i));

    dp[i+1][j+1][k][l][1]=min(dp[i+1][j+1][k][l][1],tmp+(1<<i));

    dp[i+1][j][k][l][0]=min(dp[i+1][j][k][l][0],tmp);

 #include<bits/stdc++.h>

 using namespace std;

 #define clr(a,b) memset(a,b,sizeof(a))

 typedef long long ll;

 const ll inf=0x3f3f3f3f3f3f3f3fll;

 ll a,b,c;

 ll la,lb,lc,len;

 ll dp[][][][][];

 void getlen(){len=max(log2(a)+,max(log2(b)+,log2(c)+));}

 ll getnum(ll x)

 {

     ll sum=;

     for(int i=;i<;++i) {if( x & (1ll<<i) ) sum++; }

     return sum;

 }

 void work()

 {

     dp[][][][][]=;

     for(ll i=;i<len;++i)

     {

         for(ll j=;j<=la;++j)

         {

             for(ll k=;k<=lb;++k)

             {

                 for(ll l=;l<=lc;++l)

                 {

                      long long tmp=dp[i][j][k][l][];//枚举最后一位不进位的情况

                     dp[i+][j+][k+][l+][]=min(dp[i+][j+][k+][l+][],tmp+(<<i+));

                     dp[i+][j+][k][l+][]=min(dp[i+][j+][k][l+][],tmp+(<<i));

                     dp[i+][j][k+][l+][]=min(dp[i+][j][k+][l+][],tmp+(<<i));

                     dp[i+][j][k][l][]=min(dp[i+][j][k][l][],tmp);

                     tmp=dp[i][j][k][l][];//枚举最后一位进位的情况

                     dp[i+][j+][k+][l+][]=min(dp[i+][j+][k+][l+][],tmp+(<<i+));

                     dp[i+][j][k+][l][]=min(dp[i+][j][k+][l][],tmp+(<<i));

                     dp[i+][j+][k][l][]=min(dp[i+][j+][k][l][],tmp+(<<i));

                     dp[i+][j][k][l][]=min(dp[i+][j][k][l][],tmp);

                 }

             }

         }

     }

     //cout<<len<<" "<<la<<" "<<lb<<" "<<lc<<endl;

     if(dp[len][la][lb][lc][]>=inf) printf("-1\n");

     else printf("%d\n",dp[len][la][lb][lc][]);

 }

 int main()

 {

     scanf("%lld%lld%lld",&a,&b,&c);

     getlen(); memset(dp,inf,sizeof dp);

     la=getnum(a);lb=getnum(b);lc=getnum(c);

     work();

     return ;

 }

巴特西