\(IDA ^ * \)
正直に言うと、私はこの質問には正の解決策を考え始め、その後書いた\(IDA ^ * \)を。。。
しかし、魔法はこのことです\(IDA ^ * \)がいなくても、文字列の長さがある(2500,4000 \)\、データを高速に実行されますデータビットを下された後、しかし、私は唯一の45点を測定しました。
ただ、常に最適化(IDA ^ * \)を\の過程で、私は突然、正のソリューションの実践を思い付いこれは、最初の暴力を決定するために失敗した後と思われます。
\(DP \)最初のクエリを解決するために
検討\(DP \) 、我々は設定\(F_ {I、Jを} \) 電流が第1文字列の最初であることを示している(私は\)\ 2番目の文字列のビットである(J \を\ )工程ビットの最小数。
記録場合\(nxt1_ {X、0/ 1}、nxt2_ {X、0/1} \) 2つの文字列を表す(X \)\下部後方位置\(0/1 \)それらが現れ、その後、我々は、このシフトを取得することができます:
\ [F_ {nxt1_ {I、0}、{nxt2_ J、0} =分(F_ {nxt1_ {I、0}、{nxt2_ J、0}、F_ {I、J})\]
\ [F_ {nxt1_ {I、1}、{nxt2_ J、1}} =分(F_ {nxt1_ {I、1}、{nxt2_ J、1}}、F_ {I、J})\]
これは尋ね最初に対処します。
\(BFS \)第二の問い合わせを解決します
我々が考える場合には、\(DP \)録音\を(LSTの\は)への転送の位置を示し、あなたが出力をプログラムすることができます。
しかし、問題は辞書順最小、コモン必要です\(DP \)または\(DFS \)の形で\(DP \)は、この条件を満たすことができません。
だから我々は考えることができます\(BFS \) 。
よると、\(BFS \)の次である\(DP \) 、我々は満足して辞書的に避けられない最低限の条件を保証することができます。
コード
#include<bits/stdc++.h>
#define Tp template<typename Ty>
#define Ts template<typename Ty,typename... Ar>
#define Reg register
#define RI Reg int
#define Con const
#define CI Con int&
#define I inline
#define W while
#define N 4000
using namespace std;
int n,m;string s1,s2;
class DpSolver//BFS+DP
{
private:
string ans;short qx[(N+2)*(N+2)+5],qy[(N+2)*(N+2)+5],nxt1[N+5][2],nxt2[N+5][2];
short f[N+5][N+5],gx[N+5][N+5],gy[N+5][N+5];bool glst[N+5][N+5];
public:
I void Solve()
{
RI i,j,x,y,H=1,T=0,p[2];s1="%"+s1,s2="%"+s2;
for(p[0]=p[1]=n+1,i=n+1;~i;--i) nxt1[i][0]=p[0],nxt1[i][1]=p[1],p[s1[i]&1]=i;//初始化nxt1
for(p[0]=p[1]=m+1,i=m+1;~i;--i) nxt2[i][0]=p[0],nxt2[i][1]=p[1],p[s2[i]&1]=i;//初始化nxt2
for(i=0;i<=n+1;++i) for(j=0;j<=m+1;++j) f[i][j]=m+1;f[0][0]=0,qx[++T]=0,qy[T]=0;//初始化f数组和BFS队列
W(H<=T) i=qx[H],j=qy[H++],//取出队首
f[x=nxt1[i][0]][y=nxt2[j][0]]==m+1&&(qx[++T]=x,qy[T]=y),//未访问过就入队
f[i][j]+1<f[x][y]&&(f[x][y]=f[i][j]+1,gx[x][y]=i,gy[x][y]=j,glst[x][y]=0),//更新f和g
f[x=nxt1[i][1]][y=nxt2[j][1]]==m+1&&(qx[++T]=x,qy[T]=y),//未访问过就入队
f[i][j]+1<f[x][y]&&(f[x][y]=f[i][j]+1,gx[x][y]=i,gy[x][y]=j,glst[x][y]=1);//更新f和g
x=n+1,y=m+1;W(x||y) ans=(char)(glst[x][y]+48)+ans,i=gx[x][y],j=gy[x][y],x=i,y=j;//倒着找答案
cout<<ans<<endl;//输出答案
}
}D;
int main()
{
freopen("notme.in","r",stdin),freopen("notme.out","w",stdout);
cin>>s1>>s2,s1.length()>s2.length()&&(swap(s1,s2),0),n=s1.length(),m=s2.length();
return D.Solve(),0;
}