@ysner
2018-05-25T22:30:03.000000Z
字数 1556
阅读 2161
DP
我们要求出第种方案,莫过于看其前面什么时候有种方案。
于是,我们要求出每种情况的方案数。
设表示第个字母中,已分段,第个字母为s()字母的序号 的方案数。
状态转移方程易得:(其中是下一个字母)
为了下面运算方便,要对求前缀和
然后就可以从前往后推了,若碰到的字母已知,看是否影响段数就可以了;
如果碰到的字母未知,枚举的情况,如,说明第个还在当前枚举字母情况的后面,继续枚举;否则,这一位就应该填当前枚举到的字母。
#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#define ll long long
#define re register
#define il inline
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int N=1e5+100;
char s[N];
int m,k,a[N];
ll r,n,dp[5][15][N],sum[5][15][N];
il ll gi()
{
re ll x=0,t=1;
re char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if(ch=='-') t=-1,ch=getchar();
while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
return x*t;
}
il void wri(re int x)
{
if(x<0) putchar('-'),x=-x;
if(x>9) wri(x/10);
putchar(x%10+'0');
}
int main()
{
n=gi();m=gi();r=gi();
scanf("%s",s+1);
fp(i,1,n)
if(s[i]=='A') a[i]=1;else if(s[i]=='C') a[i]=2;else if(s[i]=='G') a[i]=3;else if(s[i]=='T') a[i]=4;
if(a[n]) dp[a[n]][1][n]=1;
else fp(i,1,4) dp[i][1][n]=1;
fq(i,n-1,1)
fp(j,1,4)
if(!a[i]||a[i]==j)
fp(k,1,m)
fp(l,1,4)
dp[j][k][i]+=dp[l][k-(j>l)][i+1];
fp(i,1,4)
fp(j,1,n)
fp(k,1,m)
sum[i][k][j]=sum[i][k-1][j]+dp[i][k][j];
re int las=0;
fp(i,1,n)
if(a[i])
{
if(a[i]<las) --m;
las=a[i];
putchar(s[i]);
}
else
{
re int j;
for(j=1;j<=4&&r>sum[j][m-(j<las)][i];j++) r-=sum[j][m-(j<las)][i];
if(j==1) putchar('A');if(j==2) putchar('C');if(j==3) putchar('G');if(j==4) putchar('T');
if(j<las) --m;
las=j;
}
puts("");
return 0;
}