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. SET key value [EX seconds] [PX milliseconds] [NX|XX]

    SET key value [EX seconds] [PX milliseconds] [NX|XX] 可用版本: >= 1.0.0 时间复杂度: O(1) 将字符串值 value 关联到 k ...

  2. LInux因为缺失网关出现Name or service not known的解决方法

    笔者使用的VMware和CentOS 7.0.在安装完镜像包后,便开始配置静态ip.命令如下 vi /etc/sysconfig/network-scripts/ifcfg-ens33 将BOOTPR ...

  3. 如何申请高德地图用户Key

    打开网页https://lbs.amap.com/,进入高德开发平台. 单击箭头处[注册],打开注册页面.(如果您已注册为高德地图开发者可跳过此步骤,直接登录即可). 选择[成为个人开发者],如果您是 ...

  4. GoLang的概述

    GoLang的概述 1.什么是程序 完成某个功能的指令的集合 2.Go语言的诞生小故事 2.1. 开发团队-三个大牛 2.2.Google创造Golang的原因 2.3.Golang 的发展历程 20 ...

  5. JDK提供的并发工具类

    1.CountDownLatch await(),进入等待的状态 countDown(),计数器减一 应用场景:启动三个线程计算,需要对结果进行累加. /** * * CountDownLatch D ...

  6. Intellij IDEA最全的热键表（default keymap）

    Editing Ctrl + Space Basic code completion (the name of any class, method or variable) Ctrl + Shift ...

  7. 新安装jenkins要记录的点。

    突然很早之前部署的jenkins 由于版本低存在漏洞,被人搞了,在docker中新安装jenkins出现了各种问题,避免跳坑. 不能正确安装插件, Jenkins -> Plugin Manag ...

  8. Redis AOF持久化（二）

    1.AOF持久化的配置 AOF持久化,默认是关闭的,默认是打开RDB持久化 appendonly yes,可以打开AOF持久化机制,在生产环境里面,一般来说AOF都是要打开的,除非你说随便丢个几分钟的 ...

  9. js获取对象的属性个数

    for (var i = 0; i < dt.length; i++) { if (Object.keys(dt[i]).length <= 1) { dt.splice(i, 1); i ...

  10. window 包管理器--Chocolatey

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