NOIP2017提高组 模拟赛 27（总结）

最新推荐文章于 2020-10-21 19:34:25 发布
最新推荐文章于 2020-10-21 19:34:25 发布
阅读量503 收藏
点赞数
分类专栏： noip2017模拟总结
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/kekxy/article/details/78301883
版权

NOIP2017提高组 模拟赛 27（总结）

第一题 回文数 （推公式+快速幂）

【题目描述】
【解题思路】

　　长度为1的字符串，s[1]=9。
　　长度为3的字符串，s[3]=10*9。
　　长度为5的字符串，s[5]=10*10*9。
　　长度为7的字符串，s[7]=10*10*10*9。
　　……
　　所以，长度为2k+1的字符串，s[]=9*10^k
　　先把9忽略掉。
　　那么就是他们的系数就是

k01234567……
系数x13579111315……
系数x(乘以10)013579111315

　　最后
　　 sum=9∑nk=0x[k]∗10k
　　
　　 cha=2∗(∑nk=010k)−1
　　根据等比数列求和：
　　 cha=2∗(10n+1−1)/9−1
　　
　　 10sum−sum=9sum=9∗n∗10n+1−9∗cha
　　 sum=n∗10n+1−cha
　　p=233333不是质数，用欧拉函数求出phi(p)
　　9的逆元(%p)= 9phi(p)−1mod　p
　　或者是枚举i(1~p)，9*i%p==1，那么i就是9的逆元(%p)。
　　用快速幂随便搞搞……

【代码】
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>

using namespace std;

typedef long long ll;

int T,n;
const ll mods=233333;
ll s,ni;
char ch;
int st[50];

ll Pow(ll x,ll y)
{
    ll s=1ll;
    for(;y;y>>=1,x=x*x%mods)
    if(y&1) s=(s*x)%mods;
    return s;
}

/*int read()
{
    int yu=0;
    while((ch=getchar())!='\n')
    {
        yu=(yu<<3)+(yu<<1)+ch-'0';
    }
    return yu;
}

void print()
{
    int le=0;
    while(s)
    {
        st[++le]=s%10;
        s/=10;
    }
    for(int i=le;i>=1;i--) putchar('0'+st[i]);
    putchar('\n');
}*/

int main()
{
    freopen("2255.in","r",stdin);
    freopen("2255.out","w",stdout);
    scanf("%d",&T);
    //T=read();
    ni=Pow(9,232320-1);
    for(register int h=1;h<=T;++h)
    {
        scanf("%d",&n);
        //n=read();
        n=(n+1)>>1;
        ll ty=Pow(10,n);
        ll cy=ty*((n-1)<<1|1)%mods;
        ty=2*(ty-1+mods)%mods;
        ty=(ty*ni%mods-1+mods)%mods;
        s=(cy-ty+mods)%mods;
        //print();
        printf("%lld\n",s);
    }
    return 0;
}

第二题 购物 （DP）

【题目描述】
【解题思路】

　　一开始，看错题了，额，把题目想复杂了，以为是行动的总和（第一次扫）的平方（也就是不必连续，其实是要求连续的）……
　　首先，若一段区间A，完全包含区间B，那么区间B是完全没有的。
　　那么，预处理出左端点为L，所能延伸的最长的R，即far[L]。
　　所以最多也只会有n区间而已（n=5000）。
　　预处理出前缀和sum[]
　　假设，先选区间[3,5]，再选区间[1,4]，那么就是[1,2][3,5]。
　　如果先选区间[1,4]，再选区间[3,5]，那就是[1,4][5]。
　　因为happy值是非负的，所以，若存在区间[a,far[a]]，那么就可以选区间[a,a]~[a,far[a]]。
　　为什么？
　　平方的和不大于和的平方。所以转移之后所得到的最优解中的某个断点(即这个点之前是一次购买，之后是另一次购买)一定是对应的两个区间中的某一个的端点。这样是一定可以找到一种符合要求的购买顺序的。

【代码】
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>

#define imax(a,b) ((a>b)?(a):(b))

using namespace std;

typedef long long ll;

const int N=6000;
int n,m,a[N],far[N];
ll sum[N],ans,f[N];

int main()
{
    freopen("2256.in","r",stdin);
    freopen("2256.out","w",stdout);
    scanf("%d%d",&n,&m); ans=0ll;
    for(int i=1;i<=n;i++) scanf("%d",&a[i]),sum[i]=sum[i-1]+a[i];
    for(int i=1;i<=m;i++)
    {
        int al,ar;
        scanf("%d%d",&al,&ar);
        far[al]=imax(far[al],ar);
    }
    for(int i=1;i<=n;i++)
    if(far[i])
    {
        for(int j=i+1;j<=far[i];j++)
        if(far[j]<=far[i]) far[j]=0;
    }
    f[0]=0ll; int mx=0ll;
    for(int i=1;i<=n;i++)
    {
        f[i]=imax(f[i],f[i-1]);
        far[i]=imax(far[i],far[i-1]);
        for(int j=i;j<=far[i];j++)
        {
            ll yu=sum[j]-sum[i-1];
            f[j]=imax(f[j],f[i-1]+yu*yu);
        }
        ans=imax(f[i],ans);
    }
    /*for(int i=1;i<=n;i++)
    {
        ll mx=0ll; int mi=0;
        for(int j=1;j<=n;j++)
        {
            ll yu=sum[far[j]]-sum[j-1];
            if(yu>mx)
            {
                mx=yu;
                mi=j;
            }
        }
        ans+=mx*mx;
        for(int j=1;j<=n;j++)
        {
            if(j<mi && far[j]>=mi) far[j]=mi-1;
        }
        int fl=0;
        for(int j=1;j<=n;j++)
        if(j>mi && j<=far[mi])
        {
            fl=imax(fl,far[j]);
            far[j]=0;
        }
        far[far[mi]+1]=fl;
        far[mi]=0;
    }*/
    printf("%lld\n",ans);
    return 0;
}

第三题 宗教 （SAM+LCA(RMQ优化)）

【题目描述】
【解题思路】

　　①将s2串反过来，那么，设st是s1的子串且长度与s2相同，st的结尾串的位置为len（下标从0开始）。
　　st与s2的相同字符数就是 ∑i<=leni=0[s1[i]==st[len−i]]
　　这显然是一个卷积的形式。
　　那么，直接枚举26个字符，把等于当前枚举的字符的位置设为1，其他为0。做一次FFT，每一个len就对应着不同st串，累加相同的次数。
　　时间复杂度：O(26*n log n)
　　②再次观察，k≤100。
　　判断st与s2时，直接找k次的最长公共前缀。
　　那么，如何快速找出最长的公共前缀呢？
　　因为SAM的parents树就是反串的后缀树
　　那么，只需要把s1的反串和s2的反串连起来，中间加个特殊符号。跑一次SAM。
　　因此，s1[now1]和s2[now2]开始匹配。
　　若s1[now1]!=s2[now2]，则now1++,now2++（直到相同为止）
　　用原串的位置所对应的在SAM的状态，找到parents树上的LCA，也就是最长的公共前缀，长度为dep[LCA]。
　　下一次从s1[now1+dep[LCA]]和s2[now2+dep[LCA]]处开始匹配。
　　至于LCA，可以用RMQ优化到O(1)。
　　
　　这里写图片描述
　　遍历整棵树的每一条边。如上图，遍历的结果就是1 2 4 2 5 6 7 6 5 2 1 3
　　边数为m，则最多遍历的点数为2m。
　　假设w[x]＜w[y]（w[i]为第一次遍历到点i的时间戳）
　　那么，点x和点y的LCA就是遍历结果中的w[x]~w[y]之间的dep最小的点。
　　例如，x=4，y=6时。w[x]=3，w[y]=6。之间的点为2 5，点2的dep最小，因此，点2是x和y的LCA。
　　因为w[x]~w[y]就是x走到y的路经+路径上的点的子树。路径上的点的子树是没用的（子树的dep必定大于路径上的点的dep），所以dep最小的点必定是x和y的LCA。
　　也就是RMQ维护区间最小值。
　　
　　时间复杂度：O(n*100+n log n)
　　Accept（额，加了register和inline才过）

【代码】
//SAM+LCA(RMQ优化)
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>

#define imax(a,b) ((a>b)?(a):(b))
#define imin(a,b) ((a<b)?(a):(b))

using namespace std;

typedef long long ll;

const int M=600600;
const int N=201000;
char s1[N],s2[N];
int n,m,kk,ans;
int Td,root,last;
int to[M],ne[M],h[M],tt;
int wei[M];
int tis[55],lg2[M<<1];
int rmq[M<<1][22],tim[M],dfstime;

int op[M],tail;

struct state { int par,val,go[27]; } T[M];

void Insert(int w,int yi)
{
    int p=last;
    T[++Td].val=T[p].val+1;
    int np=Td;
    last=wei[yi]=np;
    for(;p && T[p].go[w]==0;p=T[p].par) T[p].go[w]=np;
    if(p==0) T[np].par=root; else
    {
        int q=T[p].go[w];
        if(T[p].val+1==T[q].val) T[np].par=q; else
        {
            T[++Td].val=T[p].val+1;
            int nq=Td;
            memcpy(T[nq].go,T[q].go,sizeof(T[q].go));
            T[nq].par=T[q].par;
            T[q].par=nq;
            T[np].par=nq;
            for(;p && T[p].go[w]==q;p=T[p].par) T[p].go[w]=nq;
        }
    }
}

void addedge(int a,int b) { to[++tt]=b; ne[tt]=h[a]; h[a]=tt; }

/*void dfs(int x)
{
    tim[x]=++dfstime;
    rmq[dfstime][0]=x;
    for(int p=h[x];p;p=ne[p])
    {
        dfs(to[p]);
        rmq[++dfstime][0]=x;
    }
}*/

void dfs(int x)
{
    tail=1; op[1]=x;
    while(tail)
    {
        rmq[++dfstime][0]=op[tail];
        if(h[op[tail]])
        {
            int yt=op[tail];
            op[++tail]=to[h[yt]];
            h[yt]=ne[h[yt]];
        } else tail--;
    }
    for(int i=1;i<=dfstime;i++) 
    if(!tim[rmq[i][0]]) tim[rmq[i][0]]=i;
}

inline int ance(int x,int y) { return ((T[x].val<T[y].val)?(x):(y)); }

inline int len(int x,int y)
{
    x=tim[wei[n-x]];
    y=tim[wei[n+m-y]];
    if(x>y) swap(x,y);
    int ly=lg2[y-x];
    return (T[ance(rmq[x][ly],rmq[y-tis[ly]+1][ly])].val);
}

void ready()
{
    scanf("%s",s1); Td=1; tt=1;
    T[1].val=0; root=1; last=1;
    n=strlen(s1);
    for(int i=n-1;i>=0;i--) Insert(s1[i]-'a',n-i);
    scanf("%s",s2);
    m=strlen(s2);
    Insert(26,0);
    for(int i=m-1;i>=0;i--) Insert(s2[i]-'a',n+m-i);
    scanf("%d",&kk);
    for(int i=2;i<=Td;i++) addedge(T[i].par,i);
    dfstime=0; dfs(1);
    tis[0]=1;
    for(int i=1;i<=25;i++) tis[i]=tis[i-1]<<1;
    lg2[1]=0;
    for(int i=2;i<=dfstime;i++)
    if(tis[lg2[i-1]+1]<=i) lg2[i]=lg2[i-1]+1; else lg2[i]=lg2[i-1];
    for(int j=1;j<=21;j++)
    for(int i=1;i<=dfstime;i++)
    {
        rmq[i][j]=rmq[i][j-1];
        if(i+tis[j-1]<=dfstime) rmq[i][j]=ance(rmq[i][j-1],rmq[i+tis[j-1]][j-1]);
    }
}

int main()
{
    freopen("2257.in","r",stdin);
    freopen("2257.out","w",stdout);
    ready();
    for(register int i=0;i<=n-m;i++)
    {
        int now1=i,now2=0,df=0;
        while(now2<m)
        {
            while(s1[now1]!=s2[now2] && now2<m)
            {
                df++;
                if(df>kk) break;
                now1++; now2++;
            }
            if(df>kk || now2==m) break;
            int go=len(now1,now2);
            now1+=go; now2+=go;
        }
        if(df<=kk) ans++;
    }
    printf("%d\n",ans);
    return 0;
}
//FFT 常数太大 TLE
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>

#define imax(a,b) ((a>b)?(a):(b))
#define imin(a,b) ((a<b)?(a):(b))

using namespace std;

typedef long long ll;

const ll mods=1000000007;
const int N=201000;
const double pi=acos(-1.0);
char st[N];
int d[N],w[N],n,m,kk,k,th[N<<1],sum;

struct Complex
{
    double real,image;
    Complex() {}
    Complex(double _real,double _image)
    {
        real=_real; image=_image;
    }
    friend Complex operator + (Complex A,Complex B) { return Complex(A.real+B.real,A.image+B.image); }
    friend Complex operator - (Complex A,Complex B) { return Complex(A.real-B.real,A.image-B.image); }
    friend Complex operator * (Complex A,Complex B) { return Complex(A.real*B.real-A.image*B.image,A.image*B.real+A.real*B.image); }
} A[N<<2],B[N<<2];
int rev[N<<2];

void FFT(Complex *A,int n,int DFT)
{
    for(int i=0;i<n;i++) if(i<rev[i]) swap(A[i],A[rev[i]]);
    for(int s=1;(1<<s)<=n;s++)
    {
        int mi=(1<<s);
        Complex wn=Complex(cos(2*pi/mi),DFT*sin(2*pi/mi));
        for(int t=0;t<n;t+=mi)
        {
            Complex w=Complex(1,0);
            for(int j=0;j<(mi>>1);j++)
            {
                Complex u=A[t+j],v=w*A[t+j+(mi>>1)];
                A[t+j]=u+v; A[t+j+(mi>>1)]=u-v;
                w=w*wn;
            }
        }
    }
    if(DFT==-1) for(int i=0;i<n;i++) A[i].real/=n,A[i].image/=n;
}

int main()
{
    freopen("2257.in","r",stdin);
    freopen("2257.out","w",stdout);
    scanf("%s",st);
    n=strlen(st);
    for(int i=0;i<n;i++) d[i]=st[i]-'a'+1;
    scanf("%s",st);
    m=strlen(st);
    for(int i=0;i<m;i++) w[i]=st[i]-'a'+1;
    scanf("%d",&kk);
    for(int i=1;i<=(m>>1);i++) swap(w[i-1],w[m-i]);
    n--; m--;
    int yn=n+m,nn; k=0;
    for(nn=1;nn<=yn;nn<<=1) k++;
    for(int i=0;i<nn;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<(k-1));

    for(int h=1;h<=26;h++)
    {
        for(int i=0;i<=nn;i++)
        {
            A[i].real=(d[i]==h)?1:0;
            B[i].real=(w[i]==h)?1:0;
            A[i].image=B[i].image=0;
        }
        FFT(A,nn,1); FFT(B,nn,1);
        for(int i=0;i<=nn;i++) A[i]=A[i]*B[i];
        FFT(A,nn,-1);
        for(int i=m;i<=n;i++)
        th[i]+=(int)(A[i].real+0.4);
    }
    for(int i=m;i<=n;i++)
    if(th[i]>=(m+1-kk)) sum++;
    printf("%d\n",sum);
    return 0;
}
kekxy
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
NOIP提高组模拟赛---好文章
聂泽宇的blog
07-13 175
NOIP 好文章 题解 这道题呢看起来就貌似是一个纯的HASH，只需要枚举起点然后呢再来一个set自动去重，然后输出set的size就OK，理论上是这样的。不过呢鉴于各种各样的hash冲突，这道题哪怕是用双hash也只有七十分 ...
NOIP2017提高组 模拟赛24（总结）
kxy的博客
10-11 358
##NOIP2017提高组 模拟赛24（总结） ###第一题 Karen与游戏 （贪心） ###第二题 puzzle （期望树形DP） ###第三题 Karen与测试 （找规律，快速幂）
NOIP2017提高组 模拟赛13（总结）
kxy的博客
09-12 835
NOIP2017提高组 模拟赛13（总结）第一题 函数 第二题 部落 部落的基地里已经有N个建筑设施受到了严重的损伤，如果不尽快修复的话，这些建筑设施将会完全毁坏。第三题 sunny图 codeforces 603E
纪中 2020.06.20【NOIP提高组】模拟 总结反思
CJZ的博客
06-22 201
文章目录T1 【NOIP2013模拟联考5】小麦亩产一千八(kela)比赛时之后T2 【NOIP2013模拟联考6】秀姿势(sugata)比赛时之后【NOIP2013模拟联考5】休息(rest) [题解](https://blog.csdn.net/qq_39005980/article/details/106910787)比赛时之后总结 T1 【NOIP2013模拟联考5】小麦亩产一千八(kela) 比赛时 水题一道，不过居然花了我30分钟～～～ 预估：100 之后 实际：100 居然有人没做出来（这也
[NOIP2017模拟]小店购物
g19zwk的博客
11-08 426
2017.11.7 T3 2046
noip提高组模拟题5
07-21
【压缩包子文件的文件名称列表】"BZOJ NoipRound 5"可能表示这是BZOJ（Blue Zoo Online Judge）平台上的第五次NOIP模拟赛。BZOJ是一个在线的编程题解平台，它提供各种编程题目供用户练习，并且支持代码提交和自动...
noip提高组模拟题6
07-21
"noip模拟题"表明这是针对NOIP比赛的模拟练习题目，而“提高组”则表示这些题目难度相对较高，适合已经有一定基础的信息学竞赛选手。 【压缩包子文件的文件名称列表】"BZOJ NoipRound 6"可能代表这个压缩包内的文件...
noip提高组模拟题
08-14
【标题】"noip提高组模拟题"涉及的是全国青少年信息学奥林匹克竞赛（NOIP）的提高组模拟试题。NOIP是一项针对中学生的信息技术竞赛，旨在培养学生的计算机编程、算法设计和问题解决能力。提高组是NOIP中的一个级别，...
noip提高组最新模拟题
08-14
【标题】"noip提高组最新模拟题"指的是全国青少年信息学奥林匹克联赛（NOIP）提高组的模拟试题集。NOIP是中国计算机学会举办的一项针对中学生的信息技术竞赛，旨在提高学生的计算机编程能力和算法设计能力。提高组是...
noip提高组模拟题7
07-21
【标签】"提高组"和"noip模拟题"进一步明确了这个资源的定位。提高组是NOIP竞赛中的一个级别，通常针对有一定基础和经验的参赛者，题目难度相对较高。"noip模拟题"则表明这些题目是模仿NOIP竞赛的真实题型，对于想要...
逆元详细解释（扩展欧几里得 ， 费马小定理 ， 欧拉定理）未完成
短发
09-03 589
逆元： 方程  ax≡1(mod  p)  的解称为a关于模p的逆，当gcd(a,p)==1（即a，p互质）时，方程有唯一解，否则无解。 也就是说一个数有逆元的充分必要条件是gcd(a,n)=1，此时逆元唯一存在 。 有 (a+b)%c=(a%c+b%c)%c (a*b)%c=(a%c*b%c)%c (a-b)%c=(a%c-b%c)%c 但是却没有   (a/b)%c=(a%c/b%...
【NOIP2017提高A组冲刺11.8】购物
weixin_34402090的博客
08-07 121
这个范围对DP不友好，和CF的一道C题非常像，贪心+后悔。 先使用k个优惠券购买k个q最小的（钱不购买则退出），同时把这k个p[i]-q[i]放入小根堆，然后将剩下的n-k个按p升序排序，记小根堆堆顶为top，每次比较p[i]和q[i]+top(相当于反悔，用top的代价把之前的一个优惠券用到这里) 感谢wwb的hack，虽然这题没出数据：当小根堆弹完之后就会出问题了，一个粗暴的解决方法是取...
JZOJ 5461. 【NOIP2017提高A组冲刺11.8】购物
qq_35786326的博客
09-15 178
.
JZOJ 5459 购物
sdfzsyq的博客
08-07 470
Time Limits: 1000 ms Memory Limits: 524288 KB Detailed Limits Description 小X 正困在一个密室里，他希望尽快逃出密室。 密室中有N 个房间，初始时，小X 在1 号房间，而出口在N 号房间。 密室的每一个房间中可能有着一些钥匙和一些传送门，一个传送门会单向地创造一条从房间X 到房间Y 的通道。另外，想要通过某个...
乘法逆元的三种求解方法
weixin_48842132的博客
09-02 6259
目录乘法逆元小结逆元的定义求解逆元的方法1. 快速幂测试代码2. 乘法逆元小结 乘法逆元，一般用于求（a / b）(modp)的值（p 通常为质数），是解决模意义下分数数值的必要手段。 关于求余，有以下三种规则： 加法：(a+b)%m=(a%m+b%m)%m(a+b)%m=(a%m+b%m)%m 减法：(a−b)%m=(a%m−b%m)%m(a−b)%m=(a%m−b%m)%m 乘法：(a∗b)%m=(a%m∗b%m)%m(a∗b)%m=(a%m∗b%m)%m 但是这个规则在除法不适用，所以就要用到乘法逆
[AHOI2005]洗牌
aoanping0730的博客
08-20 183
题目描述 为了表彰小联为Samuel星球的探险所做出的贡献，小联被邀请参加Samuel星球近距离载人探险活动。 由于Samuel星球相当遥远，科学家们要在飞船中度过相当长的一段时间，小联提议用扑克牌打发长途旅行中的无聊时间。玩了几局之后，大家觉得单纯玩扑克牌对于像他们这样的高智商人才来说太简单了。有人提出了扑克牌的一种新的玩法。 对于扑克牌的一次洗牌是这样定义的，将一叠N（N为偶数...
2016年省赛H题目改编--H +模逆元
转角的守望
11-17 972
Swap Time Limit: 1 Sec   Memory Limit: 128 MB Submit: 13   Solved: 2 Description Bobo has a n digits decimal number D=d1 d2…dn (It may have leading zeros). Let R(i,j)
拓展 欧几里得算法 求逆元_算法学习笔记(9)：逆元
weixin_39549312的博客
10-21 285
https://zhuanlan.zhihu.com/p/105467597在数论中，如果 ，我们就说 和 在模 意义下互为乘法逆元，记作 。逆元有什么用呢？我们常常遇到一些题目要求结果对一个大质数 取模，这是因为答案很大，出题人为了不麻烦大家写高精，就采取这样的方法。加减法和乘法对取模运算都是封闭的，所以你可以处处取模来避免溢出。一次取模和步步取模的结果是一样的但遇到除法时，就麻烦了：...
逆元的定义与计算
nowting的博客
11-18 4399
参考:博客 比如说要你求A/B%C等于多少，但是存在除法取模问题（因为(A/B）%C ！= (A%C)/(B%C)，而对于乘法却有(AB）%C ！= (A%C)(B%C) 所以要把A/B%C==&amp;amp;amp;amp;amp;gt;A*(1/B)%C 转换成A*X%C的形式，现在就是如何求X。(X在乘法上是B的逆元，意思是我们用B的逆元取代1/B就行了）这样就解决了乘法取模问题 现在讲讲逆元的定义： 在这里就是： 如果有B...
noip模拟赛 测试数据
最新发布
09-13
总之，测试数据在NOI模拟赛中起到了至关重要的作用，它既考察了选手对问题的理解和解决能力，又提高了选手编程的技巧和效率。同时，合理和恰当的测试数据也是公平竞赛的保证，确保每个参赛选手有相同的机会和条件...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值