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. Matlab 线性规划问题模型代码

    线性规划问题的基本内容 线性规划解决的是自变量在一定的线性约束条件下,使得线性目标函数求得最大值或者最小值的问题. \[ \min z=\sum_{j=1}^{n} f_{j} x_{j} \] \[ ...

  2. mysql 初级练习题

    1.题目 第一题: tb_user: User_id User_name User_phone 1 张三 13800138000 2 李四 13800138001 tb_customer: Custo ...

  3. c++项目经验分享

    1.C++的const比C语言#define更好的原因? 首先,它能够明确指定类型,有类型检查功能. 其次,可以使用C++的作用域规则将定义限制在特定的函数[常函数]或文件中. 第三,可以将const ...

  4. Django框架深入了解_05 (Django中的缓存、Django解决跨域流程(非简单请求，简单请求)、自动生成接口文档)

    一.Django中的缓存: 前戏: 在动态网站中,用户所有的请求,服务器都会去数据库中进行相应的增,删,查,改,渲染模板,执行业务逻辑,最后生成用户看到的页面. 当一个网站的用户访问量很大的时候,每一 ...

  5. 深入玩转K8S之外网如何访问业务应用

    有一个问题就是现在我的业务分配在多个Pod上,那么如果我某个Pod死掉岂不是业务完蛋了,当然也会有人说Pod死掉没问题啊,K8S自身机制Deployment和Controller会动态的创建和销毁Po ...

  6. python_封装redis_hash方法

    xshell 进入 虚拟环境 安装 redis workon py3env # 进入虚拟环境 pip install redis # 安装redis deactivate # 退出虚拟环境 简单的封装 ...

  7. aop 打印请求信息

    项目中使用 AOP 打印请求信息,打印响应信息.package com.example.aspect; import com.alibaba.fastjson.JSON;import com.goog ...

  8. ansible debugger 模块

    在搞TF(tungstenfabric)时遇到了一些错误,TF通过ansible playbook 来部署的.通常情况下遇到错误都是通过ansibale xxxx –vvv 来详细输出一下.出错的类型 ...

  9. k8s与监控--k8s部署grafana6.0

    原文参考:https://segmentfault.com/a/1190000018335241

  10. 【开发工具】- myeclipse安装主题

    你想用IDEA那样炫酷的符合90后气质的主题吗?废话不多说,按照下边步骤就可以安装像IDEA一样超级炫酷的主题. 下载主题 1.进入插件官网(http://eclipsecolorthemes.org ...