2017-10-19 NOIP模拟赛

Count（哈格朗日插值）

题解: 有个定理，另sum(x)表示小于等于x的数中与x互质的数的和

sum(x)=φ(x)*x/2    最后可知f(x)=x  （f(1)=2）  当然打表能知道。

然后就转化为了求Σi^k

然后就是拉格朗日插值法了，不在我理解范畴........

但这个博客介绍挺好哒 http://www.cnblogs.com/ECJTUACM-873284962/p/6833391.html

#include<iostream>
#include<cstdio>
#include<cstring>
#define mod 998244353
#ifdef WIN32
#define PLL "%I64d"
#else
#define PLL "%lld"
#endif
using namespace std;
int l,r,k;
long long ans;
long long mul(long long a,long long b){
    long long res=;
    while(b){
        if(b&)res=res+a;
        a=a+a;
        if(res>=mod)res-=mod;
        if(a>=mod)a-=mod;
        b>>=;
    }
    return res;
}
long long Pow(long long a,int b){
    long long res=;
    while(b){
        if(b&)res=mul(res,a);
        a=mul(a,a);
        b>>=;
    }
    return res;
}
int main(){
    freopen("count.in","r",stdin);freopen("count.out","w",stdout);
    //freopen("Cola.txt","r",stdin);
    scanf("%d%d%d",&l,&r,&k);
    for(int i=l;i<=r;i++){
        if(i==)ans+=Pow(,k);
        else ans+=Pow(i,k);
        if(ans>=mod)ans-=mod;
    }
    printf(PLL,ans);
    //cout<<ans;
}

15分 暴力

Block

/*
    题意转化为求最大的区间长度使得这段区间和减k>=0
    首先做前缀和，可知若当前到了k,i<j<k && sum[i]<sum[j]则j一定不可能比i更优
    用单调栈维护这个过程。然后倒序更新答案即可。
*/ 

#include<iostream>
#include<cstdio>
#include<cstring>

#define N 1000007

using namespace std;
long long sum[N];
int a[N],st[N],top,n,m,cnt;

inline int read()
{
    int x=,f=;char c=getchar();
    while(c>''||c<''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
    return x*f;
}

void solve(int k)
{
    top=;int res=;
    for(int i=;i<=n;i++)
    {
        sum[i]=sum[i-]+a[i]-k;
        if(!top || sum[st[top]]>sum[i]) st[++top]=i;
    }
    for(int i=n;i>=;i--)
    {
        while(top && sum[i]>=sum[st[top]]) top--;
        res=max(res,i-st[top+]);
    }
    printf("%d ",res);
}

int main()
{
    freopen("blocks.in","r",stdin);
    freopen("blocks.out","w",stdout);
    n=read();m=read();
    for(int i=;i<=n;i++) a[i]=read();
    while(m--) solve(read());
    return ;
}

100分

Biology（Trie树+lca）

#include<iostream>
#include<cstring>
#include<cstdio>
#define maxn 1000010
#define root 1
using namespace std;
int n,m,tr[maxn][],num,id[],tc=,fa[maxn],dep[maxn];
char s[];
void Insert(){
    int now=root;
    for(int i=strlen(s+);i>=;i--){
        int p=s[i]-'a';
        if(tr[now][p]==){
            tr[now][p]=++tc;
            fa[tc]=now;
            dep[tc]=dep[now]+;
        }
        now=tr[now][p];
    }
    id[num]=now;
}
int lca(int a,int b){
    while(a!=b){
        if(dep[a]<dep[b])swap(a,b);
        a=fa[a];
    }
    return a;
}
int main(){
    freopen("biology.in","r",stdin);freopen("biology.out","w",stdout);
    //freopen("Cola.in","r",stdin);
    scanf("%d%d",&n,&m);
    fa[]=;dep[]=;
    for(int i=;i<=n;i++){
        scanf("%s",s+);
        num=i,Insert();
    }
    int op,nn,x;
    for(int i=;i<=m;i++){
        scanf("%d",&op);
        if(op==){
            scanf("%s",s+);
            num++;
            Insert();
        }
        else{
            scanf("%d",&nn);
            int now;
            for(int j=;j<=nn;j++){
                scanf("%d",&x);
                if(j==)now=id[x];
                else now=lca(now,id[x]);
            }
            printf("%d\n",dep[now]);
        }
    }
    /*int sz=sizeof(tr)+sizeof(id)+sizeof(fa)+sizeof(dep);
    cout<<sz/1048576;*/
    return ;
}
/*
5 5
zzj
pri
prime
ime
owaski
2 3 1 3 5
2 2 2 3
1 actri
2 2 3 4
2 3 2 6 5
*/

80分 Trie树+暴力lca

#include<iostream>
#include<cstring>
#include<cstdio>
#define maxn 1000010
#define root 1
using namespace std;
int n,m,tr[maxn][],num,id[],tc=,fa[maxn][],dep[maxn];
char s[];
void Insert(){
    int now=root;
    for(int i=strlen(s+);i>=;i--){
        int p=s[i]-'a';
        if(tr[now][p]==){
            tr[now][p]=++tc;
            dep[tc]=dep[now]+;
            fa[tc][]=now;
            for(int j=;j<=;j++)fa[tc][j]=fa[fa[tc][j-]][j-];
        }
        now=tr[now][p];
    }
    id[num]=now;
    //for(int i=1;i<=18;i++)fa[now][i]=fa[fa[now][i-1]][i-1];
}
int lca(int a,int b){
    if(dep[a]<dep[b])swap(a,b);
    for(int i=;i>=;i--)
        if(dep[fa[a][i]]>=dep[b])a=fa[a][i];
    for(int i=;i>=;i--)
        if(fa[a][i]!=fa[b][i])
            a=fa[a][i],b=fa[b][i];
    if(a==b)return a;
    return fa[a][];
}
int main(){
    freopen("biology.in","r",stdin);freopen("biology.out","w",stdout);
    //freopen("Cola.txt","r",stdin);
    scanf("%d%d",&n,&m);
    fa[][]=;dep[]=;
    for(int i=;i<=n;i++){
        scanf("%s",s+);
        num=i,Insert();
    }
    int op,nn,x;
    for(int i=;i<=m;i++){
        scanf("%d",&op);
        if(op==){
            scanf("%s",s+);
            num++;
            Insert();
        }
        else{
            scanf("%d",&nn);
            int now;
            for(int j=;j<=nn;j++){
                scanf("%d",&x);
                if(j==)now=id[x];
                else now=lca(now,id[x]);
            }
            printf("%d\n",dep[now]);
        }
    }
    /*int sz=sizeof(tr)+sizeof(id)+sizeof(fa)+sizeof(dep);
    cout<<sz/1048576;*/
    return ;
}
/*
5 5
zzj
pri
prime
ime
owaski
2 3 1 3 5
2 2 2 3
1 actri
2 2 3 4
2 3 2 6 5
*/

100分 Trie树+倍增lca