topcode srm SRM 557
topcode srm SRM 557 div2
这是今年长春赛之前最后一次写博客了，心理还是很紧张的，毕竟这还是我第一次参加regional ，看题目心不在焉，不思考的不仔细，代码出错率搞，code时容易犯各种脑残的错误……可能这就是什么赛前综合症吧，归根结底还是水平太菜呀
250pt if else 的基础编程题就不说了
500pt 点数才10，用二进制枚举所有可能取变成魔法少女的点，然后顺着题意看这种局面是否存在，然后取最优值即可
1000pt
很明显当n为奇数数，答案为0
n为偶数的时候， U 和 D 各占一半 ,然后用点技巧就是用ac自动机来构建状态（虽然有点大财小用），但是还是比较方便的编程复查度都很低，来表示当前这个串的结尾情况
很容易想到的一种状态就是
dp1[ i ] [ j ][ state ] U的个数为i, D 的个数为 j ，结尾的状态为state 并且已经匹配了 history 串了
dp2[ i ] [ j ][ state ] U的个数为i, D 的个数为 j ，结尾的状态为state 并且没有匹配了 history 串
仔细想想其实dp1没有必要搞成三维 ，表示成dp [ i ] [ j ] 就够了
贴上代码：
500pt
bool vis[15],g[15][15],ok;
int bit[1024],n;
class IncubatorEasy
{
public:
void dfs(int u,int state)
{
//cout<<u<<" "<<state<<endl;
if(!ok) return;
for(int i=0;i<n;i++)
if(g[u][i]&&!vis[i])
{
vis[i]=1;
if(state&(1<<i)) { ok=0;break; }
dfs(i,state);
}
}
int judge(int state)
{
ok=1;
memset(vis,0,sizeof(vis));
for(int i=0;i<n;i++)
if((1<<i)&state)
{
dfs(i,state);
if(!ok) break;
}
return ok;
}
int maxMagicalGirls(vector <string> love)
{
for(int i=1;i<1024;i++)
bit[i]=bit[i>>1]+(i&1);
n=love.size();
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
g[i][j]=(love[i][j]=='Y');
int ret=0;
for(int i=(1<<n)-1;i>0;i--)
if(judge(i)) ret=max(ret,bit[i]);
return ret;
}
1000pt
const int mod=1000000009;
int chd[55][2],fail[55],sz,sw[127],word[55];
char ss[55];
int dp1[30][30][55],dp2[30][30][55];
void ins(char *s)
{
int p=0;
for(;*s;s++)
{
int id=sw[*s];
if(!chd[p][id])
{
memset(chd[sz],0,sizeof(chd[sz]));
word[sz]=0;
chd[p][id]=sz++;
}
p=chd[p][id];
}
word[p]=1;
}
int Que[55];
void ac()
{
int *s=Que,*e=Que;
if(chd[0][0]) fail[chd[0][0]]=0,*e++=chd[0][0];
if(chd[0][1]) fail[chd[0][1]]=0,*e++=chd[0][1];
while(s!=e)
{
int p=*s++;
for(int i=0;i<2;i++)
if(chd[p][i])
{
*e++=chd[p][i];
fail[chd[p][i]]=chd[fail[p]][i];
word[chd[p][i]]|=word[fail[chd[p][i]]];
}else chd[p][i]=chd[fail[p]][i];
}
}
class FoxAndMountain
{
public:
int count(int n, string h)
{
if(n&1) return 0;
//cout<<n<<" "<<h<<endl;
sw['U']=0;
sw['D']=1;
for(int i=0;i<h.length();i++) ss[i]=h[i];
ss[h.length()]='\0';
memset(chd[0],0,sizeof(chd[0]));
sz=1;
ins(ss);
ac();
//cout<<sz<<endl;
memset(dp1,0,sizeof(dp1));
memset(dp2,0,sizeof(dp2));
dp1[0][0][0]=1;
for(int i=0;i<=n/2;i++)
for(int j=0;j<=i;j++)
for(int r=0;r<sz;r++)
if(dp1[i][j][r])
{
int next=chd[r][0];
if(word[next]) dp2[i+1][j][next]=(dp2[i+1][j][next]+dp1[i][j][r])%mod;
else dp1[i+1][j][next]=(dp1[i+1][j][next]+dp1[i][j][r])%mod;
if(i!=j){
next=chd[r][1];
if(word[next]) dp2[i][j+1][next]=(dp2[i][j+1][next]+dp1[i][j][r])%mod;
else dp1[i][j+1][next]=(dp1[i][j+1][next]+dp1[i][j][r])%mod;
}
}
for(int i=0;i<=n/2;i++)
for(int j=0;j<=i;j++)
for(int r=0;r<sz;r++)
if(dp2[i][j][r])
{
int next=chd[r][0];
dp2[i+1][j][next]=(dp2[i+1][j][next]+dp2[i][j][r])%mod;
if(i!=j){
next=chd[r][1];
dp2[i][j+1][next]=(dp2[i][j+1][next]+dp2[i][j][r])%mod;
}
}
//cout<<"bug"<<endl;
int ret=0;
for(int i=0;i<sz;i++)
ret=(ret+dp2[n/2][n/2][i])%mod;
return ret;
}
我搞acm的时间其实不长了，希望我能放下紧张，好好享受剩下的每场acm比赛吧