hihocoder1457

http://hihocoder.com/problemset/problem/1457

找不重复子串的和

topo序搞一搞，用父亲更新儿子节点的val，记得乘上节点数

//#pragma comment(linker, "/stack:200000000")

//#pragma GCC optimize("Ofast,no-stack-protector")

//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

//#pragma GCC optimize("unroll-loops")

#include<bits/stdc++.h>

#define fi first

#define se second

#define mp make_pair

#define pb push_back

#define pi acos(-1.0)

#define ll long long

#define mod (1000000007)

#define C 0.5772156649

#define ls l,m,rt<<1

#define rs m+1,r,rt<<1|1

#define pil pair<int,ll>

#define pii pair<int,int>

#define ull unsigned long long

#define base 1000000000000000000

#define fio ios::sync_with_stdio(false);cin.tie(0)

using namespace std;

const double g=10.0,eps=1e-;

const int N=+,maxn=+,inf=0x3f3f3f3f,INF=0x3f3f3f3f3f3f3f3f;

inline void sub(ll &x){if(x>=mod)x-=mod;}

char s[N];

int c[N],a[N];

ll val[N],sum[N];

struct SAM{

    int last,cnt;

    int ch[N<<][],fa[N<<],l[N<<];

    int mx[N<<],tmp[N<<];

    void ins(int c)

    {

        int p=last,np=++cnt;last=np;l[np]=l[p]+;

        for(;p&&!ch[p][c];p=fa[p])ch[p][c]=np;

        if(!p)fa[np]=;

        else

        {

            int q=ch[p][c];

            if(l[p]+==l[q])fa[np]=q;

            else

            {

                int nq=++cnt;l[nq]=l[p]+;

                memcpy(ch[nq],ch[q],sizeof(ch[q]));

                fa[nq]=fa[q];fa[q]=fa[np]=nq;

                sum[nq]=1ll*c;

                for(;ch[p][c]==q;p=fa[p])ch[p][c]=nq;

            }

        }

        sum[np]=1ll*c;

    }

    SAM()

    {

        cnt=;

    }

    void build(){

        int len=strlen(s+);

        last=;

        for(int i=;i<=len;i++)ins(s[i]-'');

//        topo();

//        for(int i=1;i<=cnt;i++)mx[i]=l[i],tmp[i]=0;

    }

    void topo(){

        for(int i=;i<=cnt;i++)c[l[i]]++;

        for(int i=;i<=cnt;i++)c[i]+=c[i-];

        for(int i=;i<=cnt;i++)a[c[l[i]]--]=i;

    }

    void debug()

    {

        puts("");

        for(int i=;i<=cnt;i++)

        {

            printf("%d ",i);

            for(int j=;j<;j++)

                printf("%d ",ch[i][j]);

            printf("%d %d %lld %lld\n",fa[i],l[i],val[i],sum[i]);

        }

    }

    void faupdate()

    {

        for(int i=cnt;i;i--)

            tmp[fa[a[i]]]=max(tmp[fa[a[i]]],tmp[a[i]]);

    }

    void maupdate()

    {

        for(int i=;i<=cnt;i++)

            mx[i]=min(mx[i],tmp[i]),tmp[i]=;

    }

    void match()

    {

        int len=strlen(s+),now=,res=;

        for(int i=;i<=len;i++)

        {

            if(ch[now][s[i]-'a'])

            {

                now=ch[now][s[i]-'a'];

                res++;

            }

            else

            {

                while(now&&!ch[now][s[i]-'a'])now=fa[now];

                if(!now)now=,res=;

                else

                {

                    res=l[now]+;

                    now=ch[now][s[i]-'a'];

                }

            }

            tmp[now]=max(tmp[now],res);

        }

        faupdate();

        maupdate();

    }

    int maxlen()

    {

        int ans=;

        for(int i=;i<=cnt;i++)ans=max(ans,mx[i]);

        return ans;

    }

    void cal()

    {

        topo();

        for(int i=;i<=cnt;i++)

        {

            int p=a[i];

            val[p] =  val[p] + 1ll*sum[p] * (l[p] - l[fa[p]])%mod;

            sub(val[p]);

            for(int j=;j<;j++)

            {

                if(ch[p][j])

                {

                    val[ch[p][j]]=val[ch[p][j]]+val[p]*10ll%mod;

                    sub(val[ch[p][j]]);

                }

            }

         //   cout << val[p] << endl;

        }

       // debug();

        ll ans=;

        for(int i=cnt;i;i--)

        {

            ans=ans+val[i];

            sub(ans);

        }

        printf("%lld\n",ans);

    }

}sam;

int main()

{

    int n;

    scanf("%d",&n);

    for(int i=;i<=n;i++)

    {

        scanf("%s",s+);

        sam.build();

    }

    sam.cal();

    return ;

}

/********************

1

101

********************/

版本1

//#pragma comment(linker, "/stack:200000000")

//#pragma GCC optimize("Ofast,no-stack-protector")

//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

//#pragma GCC optimize("unroll-loops")

#include<bits/stdc++.h>

#define fi first

#define se second

#define mp make_pair

#define pb push_back

#define pi acos(-1.0)

#define ll long long

#define mod (1000000007)

#define C 0.5772156649

#define ls l,m,rt<<1

#define rs m+1,r,rt<<1|1

#define pil pair<int,ll>

#define pii pair<int,int>

#define ull unsigned long long

#define base 1000000000000000000

#define fio ios::sync_with_stdio(false);cin.tie(0)

using namespace std;

const double g=10.0,eps=1e-;

const int N=+,maxn=+,inf=0x3f3f3f3f,INF=0x3f3f3f3f3f3f3f3f;

inline void sub(ll &x){if(x>=mod)x-=mod;}

char s[N];

int c[N],a[N];

ll val[N],sum[N],sz[N];

struct SAM{

    int last,cnt;

    int ch[N<<][],fa[N<<],l[N<<];

    int mx[N<<],tmp[N<<];

    void ins(int c)

    {

        int p=last,np=++cnt;last=np;l[np]=l[p]+;

        for(;p&&!ch[p][c];p=fa[p])ch[p][c]=np;

        if(!p)fa[np]=;

        else

        {

            int q=ch[p][c];

            if(l[p]+==l[q])fa[np]=q;

            else

            {

                int nq=++cnt;l[nq]=l[p]+;

                memcpy(ch[nq],ch[q],sizeof(ch[q]));

                fa[nq]=fa[q];fa[q]=fa[np]=nq;

                for(;ch[p][c]==q;p=fa[p])ch[p][c]=nq;

            }

        }

        sum[np]=1ll*c;

    }

    SAM()

    {

        cnt=;

    }

    void build(){

        int len=strlen(s+);

        last=;

        for(int i=;i<=len;i++)ins(s[i]-'');

//        topo();

//        for(int i=1;i<=cnt;i++)mx[i]=l[i],tmp[i]=0;

    }

    void topo(){

        for(int i=;i<=cnt;i++)c[l[i]]++;

        for(int i=;i<=cnt;i++)c[i]+=c[i-];

        for(int i=;i<=cnt;i++)a[c[l[i]]--]=i;

    }

    void debug()

    {

        puts("");

        for(int i=;i<=cnt;i++)

        {

            printf("%d ",i);

            for(int j=;j<;j++)

                printf("%d ",ch[i][j]);

            printf("%d %d %lld %lld\n",fa[i],l[i],val[i],sum[i]);

        }

    }

    void faupdate()

    {

        for(int i=cnt;i;i--)

            tmp[fa[a[i]]]=max(tmp[fa[a[i]]],tmp[a[i]]);

    }

    void maupdate()

    {

        for(int i=;i<=cnt;i++)

            mx[i]=min(mx[i],tmp[i]),tmp[i]=;

    }

    void match()

    {

        int len=strlen(s+),now=,res=;

        for(int i=;i<=len;i++)

        {

            if(ch[now][s[i]-'a'])

            {

                now=ch[now][s[i]-'a'];

                res++;

            }

            else

            {

                while(now&&!ch[now][s[i]-'a'])now=fa[now];

                if(!now)now=,res=;

                else

                {

                    res=l[now]+;

                    now=ch[now][s[i]-'a'];

                }

            }

            tmp[now]=max(tmp[now],res);

        }

        faupdate();

        maupdate();

    }

    int maxlen()

    {

        int ans=;

        for(int i=;i<=cnt;i++)ans=max(ans,mx[i]);

        return ans;

    }

    void cal()

    {

        topo();

        sz[]=;

        for(int i=;i<=cnt;i++)

        {

            int p=a[i];

//            val[p] =  val[p] + 1ll*sum[p] * (l[p]-l[fa[p]])%mod;

//            sub(val[p]);

            for(int j=;j<;j++)

            {

                if(ch[p][j])

                {

                    val[ch[p][j]]=val[ch[p][j]]+val[p]*10ll%mod+1ll*j*sz[p]%mod;

                    sub(val[ch[p][j]]);

                    sz[ch[p][j]]+=sz[p];

                }

            }

         //   cout << val[p] << endl;

        }

       // debug();

        ll ans=;

        for(int i=cnt;i;i--)

        {

            ans=ans+val[i];

            sub(ans);

        }

        printf("%lld\n",ans);

    }

}sam;

int main()

{

    int n;

    scanf("%d",&n);

    for(int i=;i<=n;i++)

    {

        scanf("%s",s+);

        sam.build();

//        if(i!=n)sam.ins(10);

    }

    sam.cal();

    return ;

}

/********************

1

987654321123456789

********************/

版本2

hihocoder1457的更多相关文章

【后缀自动机】【拓扑排序】【动态规划】hihocoder1457 后缀自动机四·重复旋律7
解题方法提示小Hi:我们已经学习了后缀自动机,今天我们再来看这道有意思的题. 小Ho:好!这道题目让我们求的是若干的数字串所有不同子串的和. 小Hi:你能不能结合后缀自动机的性质来思考如何解决本题? ...

随机推荐

[golang note] 函数定义
普通函数定义 √ golang函数基本组成:关键字func.函数名.参数列表.返回值.函数体和返回语句. • 语法如下 func 函数名(参数列表) (返回值列表) { // 函数体 } • 示例如下 ...
Mysql中int和varchar类型
int类型: int bigint smallint 和 tinyint 类型,如果创建新表时没有指定 int(M) 中的M时,默认分别是 : int ------- in ...
端口状态说明 LISTENING、ESTABLISHED、TIME_WAIT及CLOSE_WAIT
TCP状态转移要点 TCP协议规定,对于已经建立的连接,网络双方要进行四次握手才能成功断开连接,如果缺少了其中某个步骤,将会使连接处于假死状态,连接本身占用的资源不会被释放.网络服务器程序要同时 ...
DevStore分享：月薪3万的程序员都避开了哪些坑
程序员薪水有高有低,有的人一个月可能拿30K.50K,有的人可能只有2K.3K.同样有五年工作经验的程序员,可能一个人每月拿20K,一个拿5K.是什么因素导致了这种差异?我特意总结了容易导致薪水低的九 ...
JSP 与 Servlet 的关系
以下摘自维基百科: Java服务器页面(JSP)是HttpServlet的扩展.由于HttpServlet大多是用来响应HTTP请求,并返回Web页面(例如HTML.XML),所以不可避免地,在编写s ...
分布式系统 SOA与中间件
在分布式系统中,有一个基础的理论 CAP,Consistency一致性 Availability可用性 Partition Tolerance分区容忍性,任何一个系统都不可能同时满足这三个条件(高富帅 ...
centos7防火墙的简单配置介绍
centos7版本 1.查看已开放的端口(默认不开放任何端口) firewall-cmd --list-ports 2.开启80端口 firewall-cmd --zone=public(作用域) - ...
centos配置jdk的环境变量
1.首先呢,centos下的JDK环境配置分两种情况,一直是root用户级别的jdk配置,另一种是其他用户组级别的配置.这里讲解的是root用户级别的配置. 我们已经下载解压好了jdk的目录.如下 2 ...
C#打开文件资源管理器
需求: 需要打开windows的文件资源管理器进行浏览文件. 方法: 利用命令提示符(cmd)中输入explorer.exe命令即可打开文件资源管理器代码实现: 1.引用 using System. ...
C teaching
L1 Preparation 1. environment codeblocks 2.variable, const 3. data type int "%d" char &quo ...

hihocoder1457

hihocoder1457的更多相关文章

随机推荐

热门专题