hdu 4733 G(x) (dp+枚举)

题意：定义G(x)=x xor (x>>1)，给两个二进制数，求有多少个满足条件的数，其中有些为不确定，用？表示。

思路：写一下x和x+1就会发现二者的区别:

x:    ????01111

x+1:????10000

那么就可以进行dp，dp[i][j]表示前i位满足条件并且第i为是j有多少种情况，那么接下来只需要枚举一下x从右到左第一位的0就行了，由于再这一位后面异或后两个串都是1000...的形式，因此可以判断一下时候能构成这样的形式就行了。这样前面还算好写，至于只有一组时输出方案的问题，最开始想根据dp数组来求，但是发现这个有可能会有多种情况，行不通，由于这个运算的特殊性，如果确定了某一位，那么这一位异或前一位等于G(x)相对应的位，因此，可以根据这个算出前一位。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<bitset>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-6
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=100000+10;
const int mod=1000000000+7;
char A[maxn],B[maxn];
int n,lsts[maxn][2],zero[maxn];
ll dp[maxn][2];
bool check(int i)
{
    return (i+1>n)||(A[i+1]!='0'&&B[i+1]!='0'&&zero[i+2]);
}
void solve()
{
    lsts[0][0]=dp[0][0]=1;
    lsts[0][1]=dp[0][1]=0;
    for(int i=1;i<=n;++i)
    {
        if(A[i]=='?'&&B[i]=='?')
        {
            dp[i][0]=(dp[i-1][0]+dp[i-1][1])%mod;
            dp[i][1]=(dp[i-1][0]+dp[i-1][1])%mod;
            lsts[i][0]=lsts[i-1][0]+lsts[i-1][1];
            lsts[i][1]=lsts[i-1][0]+lsts[i-1][1];
            if(lsts[i][0]>2) lsts[i][0]=2;
            if(lsts[i][1]>2) lsts[i][1]=2;
        }
        else if(A[i]=='?'||B[i]=='?'||(A[i]==B[i]))
        {
            if(A[i]=='0'||B[i]=='0')
            {
                dp[i][0]=dp[i-1][0];
                dp[i][1]=dp[i-1][1];
                lsts[i][0]=lsts[i-1][0];
                lsts[i][1]=lsts[i-1][1];
            }
            else
            {
                dp[i][0]=dp[i-1][1];
                dp[i][1]=dp[i-1][0];
                lsts[i][0]=lsts[i-1][1];
                lsts[i][1]=lsts[i-1][0];
            }
        }
        else
        {
            dp[i][0]=dp[i][1]=0;
            lsts[i][0]=lsts[i][1]=0;
        }
    }
    zero[n+1]=1;
    for(int i=n;i>=1;--i)
    {
        if(A[i]=='1'||B[i]=='1')
            zero[i]=0;
        else zero[i]=zero[i+1];
    }
    ll sum=0,tmp;
    int dt=0,pos=-1,dtmp;
    for(int i=1;i<=n;++i)
    {
        if(A[i]=='?'&&B[i]=='?')
        {
            tmp=(dp[i-1][0]+dp[i-1][1])%mod;
            dtmp=lsts[i-1][0]+lsts[i-1][1];
            if(dtmp&&check(i))
            {
                pos=i;
                sum=(sum+tmp)%mod;
                dt+=dtmp;
                if(dt>2) dt=2;
            }
        }
        else if(A[i]=='?'||B[i]=='?'||(A[i]!=B[i]))
        {
            if(A[i]=='0'||B[i]=='1')
            {
                tmp=dp[i-1][0];
                dtmp=lsts[i-1][0];
            }
            else
            {
                tmp=dp[i-1][1];
                dtmp=lsts[i-1][1];
            }
            if(dtmp&&check(i))
            {
                pos=i;
                sum=(sum+tmp)%mod;
                dt+=dtmp;
                if(dt>2) dt=2;
            }
        }
    }
    if(dt==0) printf("Impossible\n");
    else if(dt>1) printf("Ambiguous %I64d\n",sum);
    else
    {
        for(int i=pos+1;i<=n;++i)
            A[i]='1',B[i]='0';
        char ch,last;
        if(A[pos]=='0'||B[pos]=='1') ch='0';
        else ch='1';
        A[pos]='0';B[pos]='1';
        for(int i=pos-1;i>=1;--i)
        {
            last=ch;
            if(A[i]=='1'||B[i]=='1') ch=(ch=='0')?'1':'0';
            A[i]=B[i]=last;
        }
        for(int i=n;i>1;--i)
        {
            if(A[i-1]==A[i]) A[i]='0';
            else A[i]='1';
            if(B[i-1]==B[i]) B[i]='0';
            else B[i]='1';
        }
        printf("%s\n",A+1);
        printf("%s\n",B+1);
    }
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int t,tcase=0;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%s",A+1);
        scanf("%s",B+1);
        n=strlen(A+1);
        printf("Case #%d:\n",++tcase);
        solve();
    }
    return 0;
}