第十周 11.1-11.7

11.1-11.3

什么都没干。

11.4

HNU 10104 病毒

又是卡了很久的题目。其实已经改出了可以AC的代码。但是不能理解差别。

后来试了一种可以理解的简单做法发现能过就先这样吧- -。

题意：给n个二进制串，问能否构造出不包含这n个串的无限长二进制串。

最后的做法是简单的丢到AC自动机里。然后dfs找一个没有单词结尾的环。

没有单词结尾是这个节点和它顺着fail到root的路上都没有单词结尾。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
const int maxnode=3e4+10,sigma_size=2;
const int maxn=3e4+10;
int n,p[maxn];
char s[maxnode];

struct Ac_auto
{
    int Next[maxnode][sigma_size];
    int fail[maxnode];
    int val[maxnode];
    int vis[maxnode];
    int sz;
    Ac_auto(){sz=1;memset(Next[0],0,sizeof(Next[0]));}
    void init()
    {
        vis[0]=val[0]=0;
        sz=1;memset(Next[0],0,sizeof(Next[0]));
    }
    int idx(char c) {return c-'0';}

    void insert(char *s)
    {
        int u=0,n=strlen(s);
        for(int i=0;i<n;i++)
        {
            int c=idx(s[i]);
            if(!Next[u][c])
            {
                memset(Next[sz],0,sizeof(Next[sz]));
                val[sz]=vis[sz]=0;
                Next[u][c]=sz++;
            }
            u=Next[u][c];
        }
        val[u]=1;
    }

    void build()
    {
        queue<int> q;
        fail[0]=0;
        for(int i=0;i<sigma_size;i++) if(Next[0][i])
        {
            fail[Next[0][i]]=0;
            q.push(Next[0][i]);
        }
        while(!q.empty())
        {
            int pos=q.front(); q.pop();
            for(int i=0;i<sigma_size;i++)
            {
                if(!Next[pos][i]) Next[pos][i]=Next[fail[pos]][i];
                else
                {
                    fail[Next[pos][i]]=Next[fail[pos]][i];
                    q.push(Next[pos][i]);
                }
            }
        }
    }

    bool dfs(int pos)
    {
        vis[pos]=1;
        for(int i=0;i<sigma_size;i++)
        {
            int u=Next[pos][i],flag=0;
            while(u)
            {
                if(val[u]) {flag=1; break;}
                u=fail[u];
            }
            if(flag) continue;
            if(vis[Next[pos][i]]||dfs(Next[pos][i])) return true;
        }
        vis[pos]=0;
        return false;
    }

}ACA;

int main(void)
{
    while(~scanf("%d",&n))
    {
        ACA.init();
        for(int i=1;i<=n;i++)
        {
            scanf("%s",s+p[i-1]);
            ACA.insert(s+p[i-1]);
            p[i]=p[i-1]+strlen(s+p[i-1]);
        }
        ACA.build();
        puts(ACA.dfs(0)?"TAK":"NIE");
    }
    return 0;
}

11.5

补BC。

HDU 5524 Subtrees

我的做法是把子树分成满二叉树和非满二叉树两种统计。

深度记为dep。如果本身就是满。那么满二叉子树有dep种。

如果自身不是满。那么左右子树必有一颗满。

左满有dep-2种。右满有dep-1种。

再统计非满的二叉子树个数。

只要找到最小的非满二叉子树。那么所有包含这棵树的都是。

#include <iostream>
#include <cstdio>
using namespace std;
typedef long long LL;
LL pow[61];

int main(void)
{
    pow[0]=1LL;
    for(int i=1;i<61;i++) pow[i]=pow[i-1]*2LL;
    LL N;
    while(~scanf("%I64d",&N))
    {
        int ans,dep=1,cnt=0;
        while(N>pow[dep]-1LL) dep++;
        LL t=N-pow[dep-1]+1;
        if(t==pow[dep-1]) ans=2*dep-1;
        else if(t<pow[dep-2]) ans=2*dep-3;
        else ans=2*dep-2;
        while(t%2LL==0) ans--,t/=2LL;
        printf("%d
",ans);
    }
    return 0;
}

HDU 5525 Product

司老大题解

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
const LL MOD=1e9+7;
const int maxn=1e5;
LL pr[maxn],cnt[maxn];
int pos[maxn];

void GetPrime()
{
    memset(pr,0,sizeof(pr));
    for(int i=2;i<=maxn;i++)
    {
        if(!pr[i]) pr[++pr[0]]=i;
        for(int j=1;j<=pr[0]&&pr[j]*i<=maxn;j++)
        {
            pr[i*pr[j]]=1;
            if(i%pr[j]==0) break;
        }
    }
    for(int i=1;i<=pr[0];i++) pos[pr[i]]=i;
    return;
}

LL qpow(LL a,LL b)
{
    LL d=1,t=a;
    while(b)
    {
        if(b%2) d=(d*t)%MOD;
        b/=2;
        t=(t*t)%MOD;
    }
    return d;
}

int main(void)
{
    GetPrime();
    int n;
    while(~scanf("%d",&n))
    {
        memset(cnt,0,sizeof(cnt));
        for(int i=1;i<=n;i++)
        {
            LL x; scanf("%I64d",&x);
            int cur=i;
            for(int j=1;j<=pr[0];j++)
            {
                if(cur<(LL)pr[j]*pr[j]) break;
                if(cur%pr[j]==0)
                {
                    LL t=0;
                    while(cur%pr[j]==0) {t++; cur/=pr[j];}
                    cnt[j]=(cnt[j]+x*t)%(2*MOD-2);
                }
            }
            if(cur>1) cnt[pos[cur]]=(cnt[pos[cur]]+x)%(2*MOD-2);
        }
        LL Q=1LL,ans=1LL;
        for(int i=1;i<=pr[0];i++) Q=Q*(cnt[i]+1)%(2*MOD-2);
        for(int i=1;i<=pr[0];i++)
        {
            LL b=(cnt[i]*Q)%(2*MOD-2)/2;
            ans=ans*qpow(pr[i],b)%MOD;
        }
        printf("%I64d
",ans);
    }
    return 0;
}

HDU 5526 Lie

按照官方题解敲的。（好像并没有多少人去补这题。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int num[1010][1010];
int val[1010][1010];
int dp[1010];

int main(void)
{
    int N;
    while(~scanf("%d",&N))
    {
        memset(num,0,sizeof(num));
        memset(val,0,sizeof(val));
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=N;i++)
        {
            int A,B;
            scanf("%d%d",&A,&B);
            if(A+B+1>N) continue;
            num[A+B+1][0]++;
            num[A+B+1][A+1]++;
        }
        for(int i=1;i<=N;i++)
        {
            for(int j=1;i*j<=N;j++)
            {
                int tmp=0;
                for(int k=1;k<=i*j;k++) tmp+=min(j,num[i][k]);
                val[i][j]=tmp;
            }
            for(int v=N;v>=0;v--)
                for(int j=1;i*j<=N;j++)
                    if(v>=i*j) dp[v]=max(dp[v],dp[v-i*j]+val[i][j]);
        }
        printf("%d
",dp[N]);
    }
    return 0;
}

补完一套BC。

11.6-11.7

什么都没干。