hdu6172

模板的简单应用

先根据题中的表达式求出前几项，再上BM，注意一下n的大小关系。

#include <bits/stdc++.h>

using namespace std;

#define rep(i,a,n) for (long long i=a;i<n;i++)

#define per(i,a,n) for (long long i=n-1;i>=a;i--)

#define pb push_back

#define mp make_pair

#define all(x) (x).begin(),(x).end()

#define fi first

#define se second

#define SZ(x) ((long long)(x).size())

typedef vector<long long> VI;

typedef long long ll;

typedef pair<long long,long long> PII;

const ll mod=1e9+;

ll powmod(ll a,ll b) {ll res=;a%=mod; assert(b>=); for(;b;b>>=){if(b&)res=res*a%mod;a=a*a%mod;}return res;}

// head

long long _,n;

namespace linear_seq

{

    const long long N=;

    ll res[N],base[N],_c[N],_md[N];

    vector<long long> Md;

    void mul(ll *a,ll *b,long long k)

    {

        rep(i,,k+k) _c[i]=;

        rep(i,,k) if (a[i]) rep(j,,k)

            _c[i+j]=(_c[i+j]+a[i]*b[j])%mod;

        for (long long i=k+k-;i>=k;i--) if (_c[i])

            rep(j,,SZ(Md)) _c[i-k+Md[j]]=(_c[i-k+Md[j]]-_c[i]*_md[Md[j]])%mod;

        rep(i,,k) a[i]=_c[i];

    }

    long long solve(ll n,VI a,VI b)

    { // a 系数 b 初值 b[n+1]=a[0]*b[n]+...

//        printf("%d\n",SZ(b));

        ll ans=,pnt=;

        long long k=SZ(a);

        assert(SZ(a)==SZ(b));

        rep(i,,k) _md[k--i]=-a[i];_md[k]=;

        Md.clear();

        rep(i,,k) if (_md[i]!=) Md.push_back(i);

        rep(i,,k) res[i]=base[i]=;

        res[]=;

        while ((1ll<<pnt)<=n) pnt++;

        for (long long p=pnt;p>=;p--)

        {

            mul(res,res,k);

            if ((n>>p)&)

            {

                for (long long i=k-;i>=;i--) res[i+]=res[i];res[]=;

                rep(j,,SZ(Md)) res[Md[j]]=(res[Md[j]]-res[k]*_md[Md[j]])%mod;

            }

        }

        rep(i,,k) ans=(ans+res[i]*b[i])%mod;

        if (ans<) ans+=mod;

        return ans;

    }

    VI BM(VI s)

    {

        VI C(,),B(,);

        long long L=,m=,b=;

        rep(n,,SZ(s))

        {

            ll d=;

            rep(i,,L+) d=(d+(ll)C[i]*s[n-i])%mod;

            if (d==) ++m;

            else if (*L<=n)

            {

                VI T=C;

                ll c=mod-d*powmod(b,mod-)%mod;

                while (SZ(C)<SZ(B)+m) C.pb();

                rep(i,,SZ(B)) C[i+m]=(C[i+m]+c*B[i])%mod;

                L=n+-L; B=T; b=d; m=;

            }

            else

            {

                ll c=mod-d*powmod(b,mod-)%mod;

                while (SZ(C)<SZ(B)+m) C.pb();

                rep(i,,SZ(B)) C[i+m]=(C[i+m]+c*B[i])%mod;

                ++m;

            }

        }

        return C;

    }

    long long gao(VI a,ll n)

    {

        VI c=BM(a);

        c.erase(c.begin());

        rep(i,,SZ(c)) c[i]=(mod-c[i])%mod;

        return solve(n,c,VI(a.begin(),a.begin()+SZ(c)));

    }

};

int main()

{

    int T;

    scanf("%d", &T);

    while(T--)

    {

        scanf("%lld", &n);

        n--;

        /*求第n项*/

        printf("%I64d\n",linear_seq::gao(VI{,,,,,,,,},n-));

    }

}

hdu6185

模板的简单应用，先写个暴力程序找出前几项，可见 铺砖问题

#include <bits/stdc++.h>

using namespace std;

#define rep(i,a,n) for (long long i=a;i<n;i++)

#define per(i,a,n) for (long long i=n-1;i>=a;i--)

#define pb push_back

#define mp make_pair

#define all(x) (x).begin(),(x).end()

#define fi first

#define se second

#define SZ(x) ((long long)(x).size())

typedef vector<long long> VI;

typedef long long ll;

typedef pair<long long,long long> PII;

const ll mod=1e9+;

ll powmod(ll a,ll b) {ll res=;a%=mod; assert(b>=); for(;b;b>>=){if(b&)res=res*a%mod;a=a*a%mod;}return res;}

// head

long long _,n;

namespace linear_seq

{

    const long long N=;

    ll res[N],base[N],_c[N],_md[N];

    vector<long long> Md;

    void mul(ll *a,ll *b,long long k)

    {

        rep(i,,k+k) _c[i]=;

        rep(i,,k) if (a[i]) rep(j,,k)

            _c[i+j]=(_c[i+j]+a[i]*b[j])%mod;

        for (long long i=k+k-;i>=k;i--) if (_c[i])

            rep(j,,SZ(Md)) _c[i-k+Md[j]]=(_c[i-k+Md[j]]-_c[i]*_md[Md[j]])%mod;

        rep(i,,k) a[i]=_c[i];

    }

    long long solve(ll n,VI a,VI b)

    { // a 系数 b 初值 b[n+1]=a[0]*b[n]+...

//        printf("%d\n",SZ(b));

        ll ans=,pnt=;

        long long k=SZ(a);

        assert(SZ(a)==SZ(b));

        rep(i,,k) _md[k--i]=-a[i];_md[k]=;

        Md.clear();

        rep(i,,k) if (_md[i]!=) Md.push_back(i);

        rep(i,,k) res[i]=base[i]=;

        res[]=;

        while ((1ll<<pnt)<=n) pnt++;

        for (long long p=pnt;p>=;p--)

        {

            mul(res,res,k);

            if ((n>>p)&)

            {

                for (long long i=k-;i>=;i--) res[i+]=res[i];res[]=;

                rep(j,,SZ(Md)) res[Md[j]]=(res[Md[j]]-res[k]*_md[Md[j]])%mod;

            }

        }

        rep(i,,k) ans=(ans+res[i]*b[i])%mod;

        if (ans<) ans+=mod;

        return ans;

    }

    VI BM(VI s)

    {

        VI C(,),B(,);

        long long L=,m=,b=;

        rep(n,,SZ(s))

        {

            ll d=;

            rep(i,,L+) d=(d+(ll)C[i]*s[n-i])%mod;

            if (d==) ++m;

            else if (*L<=n)

            {

                VI T=C;

                ll c=mod-d*powmod(b,mod-)%mod;

                while (SZ(C)<SZ(B)+m) C.pb();

                rep(i,,SZ(B)) C[i+m]=(C[i+m]+c*B[i])%mod;

                L=n+-L; B=T; b=d; m=;

            }

            else

            {

                ll c=mod-d*powmod(b,mod-)%mod;

                while (SZ(C)<SZ(B)+m) C.pb();

                rep(i,,SZ(B)) C[i+m]=(C[i+m]+c*B[i])%mod;

                ++m;

            }

        }

        return C;

    }

    long long gao(VI a,ll n)

    {

        VI c=BM(a);

        c.erase(c.begin());

        rep(i,,SZ(c)) c[i]=(mod-c[i])%mod;

        return solve(n,c,VI(a.begin(),a.begin()+SZ(c)));

    }

};

int main()

{

    while(scanf("%lld", &n) == )

    {

        /*求第n项*/

        printf("%I64d\n",linear_seq::gao(VI{,,,,,,,},n-));

    }

}

hdu6198

模板的简单应用

先写一个爆搜找出前几项，

#include<bits/stdc++.h>

using namespace std;

const int maxn =  + ;

bool vis[maxn];

int k;

int f[] = {,,,,,,,,,,,,,,,,,,,};

bool dfs(int n, int cnt)

{

    //printf("%d\n", n);

    if(n ==  && cnt == k)  return  true;

    if(cnt == k)  return false;

    for(int i = ;i <  && f[i] <= n; i++)

    {

        if(dfs(n - f[i], cnt+))  return true;  //有一个满足条件的分解即可返回

    }

    return false;

}

int main()

{

    for(k = ;k < ;k++)

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

        {

            if(!dfs(i, ))

            {

                printf("%d %d\n",k,  i);

                break;

            }

        }

    return ;

}

能很快找出前7项，这道题中够了。

#include <bits/stdc++.h>

using namespace std;

#define rep(i,a,n) for (long long i=a;i<n;i++)

#define per(i,a,n) for (long long i=n-1;i>=a;i--)

#define pb push_back

#define mp make_pair

#define all(x) (x).begin(),(x).end()

#define fi first

#define se second

#define SZ(x) ((long long)(x).size())

typedef vector<long long> VI;

typedef long long ll;

typedef pair<long long,long long> PII;

const ll mod=;

ll powmod(ll a,ll b) {ll res=;a%=mod; assert(b>=); for(;b;b>>=){if(b&)res=res*a%mod;a=a*a%mod;}return res;}

// head

long long _,n;

namespace linear_seq

{

    const long long N=;

    ll res[N],base[N],_c[N],_md[N];

    vector<long long> Md;

    void mul(ll *a,ll *b,long long k)

    {

        rep(i,,k+k) _c[i]=;

        rep(i,,k) if (a[i]) rep(j,,k)

            _c[i+j]=(_c[i+j]+a[i]*b[j])%mod;

        for (long long i=k+k-;i>=k;i--) if (_c[i])

            rep(j,,SZ(Md)) _c[i-k+Md[j]]=(_c[i-k+Md[j]]-_c[i]*_md[Md[j]])%mod;

        rep(i,,k) a[i]=_c[i];

    }

    long long solve(ll n,VI a,VI b)

    { // a 系数 b 初值 b[n+1]=a[0]*b[n]+...

//        printf("%d\n",SZ(b));

        ll ans=,pnt=;

        long long k=SZ(a);

        assert(SZ(a)==SZ(b));

        rep(i,,k) _md[k--i]=-a[i];_md[k]=;

        Md.clear();

        rep(i,,k) if (_md[i]!=) Md.push_back(i);

        rep(i,,k) res[i]=base[i]=;

        res[]=;

        while ((1ll<<pnt)<=n) pnt++;

        for (long long p=pnt;p>=;p--)

        {

            mul(res,res,k);

            if ((n>>p)&)

            {

                for (long long i=k-;i>=;i--) res[i+]=res[i];res[]=;

                rep(j,,SZ(Md)) res[Md[j]]=(res[Md[j]]-res[k]*_md[Md[j]])%mod;

            }

        }

        rep(i,,k) ans=(ans+res[i]*b[i])%mod;

        if (ans<) ans+=mod;

        return ans;

    }

    VI BM(VI s)

    {

        VI C(,),B(,);

        long long L=,m=,b=;

        rep(n,,SZ(s))

        {

            ll d=;

            rep(i,,L+) d=(d+(ll)C[i]*s[n-i])%mod;

            if (d==) ++m;

            else if (*L<=n)

            {

                VI T=C;

                ll c=mod-d*powmod(b,mod-)%mod;

                while (SZ(C)<SZ(B)+m) C.pb();

                rep(i,,SZ(B)) C[i+m]=(C[i+m]+c*B[i])%mod;

                L=n+-L; B=T; b=d; m=;

            }

            else

            {

                ll c=mod-d*powmod(b,mod-)%mod;

                while (SZ(C)<SZ(B)+m) C.pb();

                rep(i,,SZ(B)) C[i+m]=(C[i+m]+c*B[i])%mod;

                ++m;

            }

        }

        return C;

    }

    long long gao(VI a,ll n)

    {

        VI c=BM(a);

        c.erase(c.begin());

        rep(i,,SZ(c)) c[i]=(mod-c[i])%mod;

        return solve(n,c,VI(a.begin(),a.begin()+SZ(c)));

    }

};

int main()

{

    while(scanf("%lld", &n) == )

    {

        /*求第n项*/

        printf("%I64d\n",linear_seq::gao(VI{, , , , , , },n-));

//         VI res = linear_seq::BM(VI{4, 12, 33, 88, 232, 609, 1596});

//        for(int i = 1;i < res.size();i++)  printf("%lld\n", (mod-res[i]) % mod);

    }

}

P5487

要求输出递推式和第 $m$ 项，

可能用来验板子（单组，$ 1 \leq n \leq 1000, m\leq 10^9 $，时限3.5s，不开O2会T，应该是vector的原因吧）

#include <bits/stdc++.h>

using namespace std;

#define rep(i,a,n) for (long long i=a;i<n;i++)

#define per(i,a,n) for (long long i=n-1;i>=a;i--)

#define pb push_back

#define mp make_pair

#define all(x) (x).begin(),(x).end()

#define fi first

#define se second

#define SZ(x) ((long long)(x).size())

typedef vector<long long> VI;

typedef long long ll;

typedef pair<long long,long long> PII;

const ll mod=;

ll powmod(ll a,ll b) {ll res=;a%=mod; assert(b>=); for(;b;b>>=){if(b&)res=res*a%mod;a=a*a%mod;}return res;}

// head

long long _,n;

namespace linear_seq

{

    const long long N=;

    ll res[N],base[N],_c[N],_md[N];

    vector<long long> Md;

    void mul(ll *a,ll *b,long long k)

    {

        rep(i,,k+k) _c[i]=;

        rep(i,,k) if (a[i]) rep(j,,k)

            _c[i+j]=(_c[i+j]+a[i]*b[j])%mod;

        for (long long i=k+k-;i>=k;i--) if (_c[i])

            rep(j,,SZ(Md)) _c[i-k+Md[j]]=(_c[i-k+Md[j]]-_c[i]*_md[Md[j]])%mod;

        rep(i,,k) a[i]=_c[i];

    }

    long long solve(ll n,VI a,VI b)

    { // a 系数 b 初值 b[n+1]=a[0]*b[n]+...

//        printf("%d\n",SZ(b));

        ll ans=,pnt=;

        long long k=SZ(a);

        assert(SZ(a)==SZ(b));

        rep(i,,k) _md[k--i]=-a[i];_md[k]=;

        Md.clear();

        rep(i,,k) if (_md[i]!=) Md.push_back(i);

        rep(i,,k) res[i]=base[i]=;

        res[]=;

        while ((1ll<<pnt)<=n) pnt++;

        for (long long p=pnt;p>=;p--)

        {

            mul(res,res,k);

            if ((n>>p)&)

            {

                for (long long i=k-;i>=;i--) res[i+]=res[i];res[]=;

                rep(j,,SZ(Md)) res[Md[j]]=(res[Md[j]]-res[k]*_md[Md[j]])%mod;

            }

        }

        rep(i,,k) ans=(ans+res[i]*b[i])%mod;

        if (ans<) ans+=mod;

        return ans;

    }

    VI BM(VI s)

    {

        VI C(,),B(,);

        long long L=,m=,b=;

        rep(n,,SZ(s))

        {

            ll d=;

            rep(i,,L+) d=(d+(ll)C[i]*s[n-i])%mod;

            if (d==) ++m;

            else if (*L<=n)

            {

                VI T=C;

                ll c=mod-d*powmod(b,mod-)%mod;

                while (SZ(C)<SZ(B)+m) C.pb();

                rep(i,,SZ(B)) C[i+m]=(C[i+m]+c*B[i])%mod;

                L=n+-L; B=T; b=d; m=;

            }

            else

            {

                ll c=mod-d*powmod(b,mod-)%mod;

                while (SZ(C)<SZ(B)+m) C.pb();

                rep(i,,SZ(B)) C[i+m]=(C[i+m]+c*B[i])%mod;

                ++m;

            }

        }

        return C;

    }

    long long gao(VI a,ll n)

    {

        VI c=BM(a);

        c.erase(c.begin());

        rep(i,,SZ(c)) c[i]=(mod-c[i])%mod;

        return solve(n,c,VI(a.begin(),a.begin()+SZ(c)));

    }

};

int main()

{

    int k;

    vector<ll>tmp;

    scanf("%d%lld", &k, &n);

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

    {

        ll t;

        scanf("%lld", &t);

        tmp.push_back(t);

    }

     /*输出系数*/

    /*前k项递推，需要2*k项能确定*/

    VI res = linear_seq::BM(tmp);

    for(int i = ;i < res.size();i++)  printf("%lld%c", (mod-res[i]) % mod, i == res.size()- ? '\n' : ' ');

    /*输出第n项*/

    printf("%lld\n",linear_seq::gao(tmp, n));

}

参考链接：

1. https://blog.csdn.net/Anxdada/article/details/77817850

2. https://blog.csdn.net/WilliamSun0122/article/details/77926806

hdu6172&&hdu6185&&P5487——BM算法的更多相关文章

  1. BM算法　　Boyer-Moore高质量实现代码详解与算法详解

    Boyer-Moore高质量实现代码详解与算法详解 鉴于我见到对算法本身分析非常透彻的文章以及实现的非常精巧的文章,所以就转载了,本文的贡献在于将两者结合起来,方便大家了解代码实现! 算法详解转自:h ...

  2. hrbustoj 1551:基础数据结构——字符串2 病毒II（字符串匹配，BM算法练习）

    基础数据结构——字符串2 病毒IITime Limit: 1000 MS Memory Limit: 10240 KTotal Submit: 284(138 users) Total Accepte ...

  3. BM算法详解

    http://www-igm.univ-mlv.fr/~lecroq/string/node14.html http://www.cs.utexas.edu/users/moore/publicati ...

  4. Boyer-Moore(BM)算法，文本查找，字符串匹配问题

    KMP算法的时间复杂度是O(m + n),而Boyer-Moore算法的时间复杂度是O(n/m).文本查找中“ctrl + f”一般就是采用的BM算法. Boyer-Moore算法的关键点: 从右遍历 ...

  5. Berlekamp_Massey 算法 (BM算法) 学习笔记

    原文链接www.cnblogs.com/zhouzhendong/p/Berlekamp-Massey.html 前言 BM算法用于求解常系数线性递推式. 它可以在 $O(n^2)$ 的时间复杂度内解 ...

  6. 【算法】BM算法

    目录 BM算法 一. 字符串比较的分析 二.BM算法的思想 三.算法实现 BM算法 @ 一. 字符串比较的分析 如果要判定长度为\(n\)两个字符串相等,比较中要进行\(n\)比较,但是如果要判定两个 ...

  7. BM算法

    BM算法 用来求解一个数列的递推式. 即给定\(\{x_i\}\)求解一个\(\{a_i\}\),满足\(|a|=m,x_n=\sum_{i=1}^ma_i*x_{n-i}\). 考虑增量法构造. 假 ...

  8. 数据结构 BM算法

    BM算法是比KMP算法更快的字符串模式匹配算法.BM算法最好情况下的时间复杂度是O(n),KMP算法最好情况下的时间复杂度是O(n+m),两者最坏情况下的时间复杂度均是O(m*n).其中,n指目标串长 ...

  9. BM算法--串匹配

    BM(Boyer-Moore)算法,后缀匹配,是指模式串的比较从右到左,模式串的移动也是从左到右的匹配过程,一般情况比KMP算法要快.时间复杂度O(m/n) C++描述(教师版) int BM(cha ...

随机推荐

  1. 手撕面试官系列（六）：并发+Netty+JVM+Linux面试专题

    并发面试专题 (面试题+答案领取方式见侧边栏) 现在有 T1.T2.T3 三个线程,你怎样保证 T2 在 T1 执行完后执行,T3 在 T2 执行完后执行? 在 Java 中 Lock 接口比 syn ...

  2. Haystack搜索框架

    1.什么是Haystack Haystack是django的开源全文搜索框架(全文检索不同于特定字段的模糊查询,使用全文检索的效率更高 ),该框架支持Solr,Elasticsearch,Whoosh ...

  3. c，使用lib,dll

    lib使用: #include "xxx.h" // lib的头文件 #pragma comment(lib, "xxx.lib") 这样会将lib里的数据编译 ...

  4. Docker之数据卷(Data Volumes)操作

    目的: 前言 Docker宿主机和容器之间文件拷贝 数据卷 数据卷容器 前言 Docker 数据管理 在生产环境中使用 Docker ,往往需要对数据进行持久化,或者需要在多个容器之间进行 数据共享, ...

  5. Django开发之登陆和登出

    使用django自带的验证模块 1.首先使用python manage.py startapp models为当前项目添加一个应用. 2.在setting.py中INSTALLED_APPS后面添加' ...

  6. python_并行与并发、多线程

    问题一: 计算机是如何执行程序指令的? 问题二: 计算机如何实现并发的? 轮询调度实现并发执行 程序1-8轮询完成,才再CPU上运行 问题三: 真正的并行需要依赖什么? 并行需要的核心条件 多进程实现 ...

  7. Sqlserver 总结(2) 存储过程

    目录 写在前面 什么是存储过程 存储过程的优点 存储过程的分类 1.只返回单一记录集的存储过程 2.没有输入输出的存储过程 3.有返回值的存储过程 4.有输入参数和输出参数的存储过程 5.同时具有返回 ...

  8. docker build 错误 /usr/share/dotnet/sdk/2.1.801/Microsoft.Common.CurrentVersion.targets(2106,5): warning MSB3245: Could not resolve this reference

    docker dotnet Restore 的时候报错, 一度怀疑是linux的dotnet core sdk没有装好, 卸了装, 装了卸, 试了好几遍还是无效(Microsoft.Common.Cu ...

  9. window 包管理器--Chocolatey

    Chocolatey 介绍 在 Linux 下,大家喜欢用 apt-get 来安装应用程序,如今在 windows 下,大家可以使用 Chocolatey 来快速下载搭建一个开发环境. Chocola ...

  10. 约束布局ConstraintLayout

    Android新特性介绍,ConstraintLayout完全解析 约束布局ConstraintLayout用法全解析 约束布局ConstraintLayout看这一篇就够了