1. 程式人生 > >BZOJ2423 HAOI2010最長公共子序列(動態規劃)

BZOJ2423 HAOI2010最長公共子序列(動態規劃)

out color tchar efi turn name 子序列 lld 動態

  大討論。註意去重。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=0,f=1;char c=getchar();
    while (c<0||c>9) {if (c==-) f=-1;c=getchar();}
    while (c>=
0&&c<=9) x=(x<<1)+(x<<3)+(c^48),c=getchar(); return x*f; } #define N 5010 #define P 100000000 int n,m,f[2][N],g[2][N]; char a[N],b[N]; int main() { #ifndef ONLINE_JUDGE freopen("bzoj2423.in","r",stdin); freopen("bzoj2423.out","w",stdout); const char LL[]="%I64d\n"; #else
const char LL[]="%lld\n"; #endif char c=getchar(); while (c>=A&&c<=Z) a[++n]=c,c=getchar(); a[++n]=c,c=getchar(); while (c<A||c>Z) c=getchar(); while (c>=A&&c<=Z) b[++m]=c,c=getchar(); b[++m]=c,c=getchar(); for (int j=0;j<=m;j++) g[0
][j]=1; for (int i=1;i<=n;i++) { memset(f[i&1],0,sizeof(f[i&1])); memset(g[i&1],0,sizeof(g[i&1])); g[i&1][0]=1; for (int j=1;j<=m;j++) if (a[i]!=b[j]) { f[i&1][j]=max(f[i&1^1][j],f[i&1][j-1]); if (f[i&1^1][j]==f[i&1][j-1]) { if (f[i&1^1][j]==f[i&1^1][j-1]) g[i&1][j]=((g[i&1^1][j]+g[i&1][j-1])%P-g[i&1^1][j-1]+P)%P; else g[i&1][j]=(g[i&1^1][j]+g[i&1][j-1])%P; } else if (f[i&1^1][j]>f[i&1][j-1]) g[i&1][j]=g[i&1^1][j]; else g[i&1][j]=g[i&1][j-1]; } else { f[i&1][j]=f[i&1^1][j-1]+1; g[i&1][j]=g[i&1^1][j-1]; if (f[i&1][j]==f[i&1^1][j]) g[i&1][j]=(g[i&1][j]+g[i&1^1][j])%P; if (f[i&1][j]==f[i&1][j-1]) g[i&1][j]=(g[i&1][j]+g[i&1][j-1])%P; } } cout<<f[n&1][m]-1<<endl<<g[n&1][m]; return 0; }

BZOJ2423 HAOI2010最長公共子序列(動態規劃)