BZOJ5294 BJOI2018二进制(线段树)

  二进制数能被3整除相当于奇数、偶数位上1的个数模3同余。那么如果有偶数个1,一定存在重排方案使其合法;否则则要求至少有两个0且至少有3个1,这样可以给奇数位单独安排3个1。

  考虑线段树维护区间内的一堆东西,合并两节点时计算跨过区间中点的答案。可以对每个节点记录f[0/1][0/1][0/1][0/1/2]表示前/后缀,异或和为0/1,是否至少出现了两个1,出现了0/1/超过2个0。大力讨论即可。

  成功写了一晚上才不是因为要补十几面数学作业

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 100010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
    int x=0,f=1;char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
    while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
    return x*f;
}
int n,m,a[N];
struct data{int f[2][2][2][3],l,r,sum;ll ans;
}tree[N<<2];
data merge(data x,data y)
{
    data p;memset(p.f,0,sizeof(p.f));
    p.l=x.l,p.r=y.r;
    p.sum=x.sum+y.sum;
    p.ans=x.ans+y.ans;
    for (int i=0;i<2;i++)
    {
        int s1=0,s2=0;
        for (int j=0;j<2;j++)
            for (int k=0;k<3;k++)
            s1+=x.f[1][i][j][k],s2+=y.f[0][i][j][k];
        p.ans+=1ll*s1*s2;
    }
    for (int i=0;i<3;i++) 
        for (int j=2-i;j<3;j++)
            for (int u=0;u<2;u++)
                for (int v=0;v<2;v++)
                if (u|v) p.ans+=1ll*x.f[1][0][u][i]*y.f[0][1][v][j]+1ll*x.f[1][1][u][i]*y.f[0][0][v][j];
    int lone=x.sum,lzero=x.r-x.l+1-lone,rone=y.sum,rzero=y.r-y.l+1-rone;
    for (int i=0;i<2;i++) 
        for (int j=0;j<2;j++)
            for (int k=0;k<3;k++)
            p.f[0][i][j][k]+=x.f[0][i][j][k],p.f[0][i^(lone&1)][lone+i>=2||j][min(2,lzero+k)]+=y.f[0][i][j][k],
            p.f[1][i][j][k]+=y.f[1][i][j][k],p.f[1][i^(rone&1)][rone+i>=2||j][min(2,rzero+k)]+=x.f[1][i][j][k];
    return p;
}
void newpoint(int k,int x)
{
    memset(tree[k].f,0,sizeof(tree[k].f));
    tree[k].f[0][x][0][x^1]=tree[k].f[1][x][0][x^1]=1;tree[k].ans=x^1;tree[k].sum=x;
}
void build(int k,int l,int r)
{
    tree[k].l=l,tree[k].r=r;
    if (l==r){newpoint(k,a[l]);return;}
    int mid=l+r>>1;
    build(k<<1,l,mid);
    build(k<<1|1,mid+1,r);
    tree[k]=merge(tree[k<<1],tree[k<<1|1]);
}
void modify(int k,int p,int x)
{
    if (tree[k].l==tree[k].r) {newpoint(k,x);return;}
    int mid=tree[k].l+tree[k].r>>1;
    if (p<=mid) modify(k<<1,p,x);
    else modify(k<<1|1,p,x);
    tree[k]=merge(tree[k<<1],tree[k<<1|1]);
}
data query(int k,int l,int r)
{
    if (tree[k].l==l&&tree[k].r==r) return tree[k];
    int mid=tree[k].l+tree[k].r>>1;
    if (r<=mid) return query(k<<1,l,r);
    else if (l>mid) return query(k<<1|1,l,r);
    else return merge(query(k<<1,l,mid),query(k<<1|1,mid+1,r));
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj5294.in","r",stdin);
    freopen("bzoj5294.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read();
    for (int i=1;i<=n;i++) a[i]=read();
    build(1,0,n);
    m=read();
    while (m--)
    {
        int op=read();
        if (op==1)
        {
            int x=read();
            modify(1,x,a[x]^=1);
        }
        else
        {
            int l=read(),r=read();
            printf(LL,query(1,l,r).ans);
        }
    }
    return 0;
}

 

转载于:https://www.cnblogs.com/Gloid/p/10095788.html

内容概要:本文详细探讨了基于樽海鞘算法(SSA)优化的极限学习机(ELM)在回归预测任务中的应用,并与传统的BP神经网络、广义回归神经网络(GRNN)以及未优化的ELM进行了性能对比。首先介绍了ELM的基本原理,即通过随机生成输入层与隐藏层之间的连接权重及阈值,仅需计算输出权重即可快速完成训练。接着阐述了SSA的工作机制,利用樽海鞘群体觅食行为优化ELM的输入权重和隐藏层阈值,从而提高模型性能。随后分别给出了BP、GRNN、ELM和SSA-ELM的具体实现代码,并通过波士顿房价数据集和其他工业数据集验证了各模型的表现。结果显示,SSA-ELM在预测精度方面显著优于其他三种方法,尽管其训练时间较长,但在实际应用中仍具有明显优势。 适合人群:对机器学习尤其是回归预测感兴趣的科研人员和技术开发者,特别是那些希望深入了解ELM及其优化方法的人。 使用场景及目标:适用于需要高效、高精度回归预测的应用场景,如金融建模、工业数据分析等。主要目标是提供一种更为有效的回归预测解决方案,尤其是在处理大规模数据集时能够保持较高的预测精度。 其他说明:文中提供了详细的代码示例和性能对比图表,帮助读者更好地理解和复现实验结果。同时提醒使用者注意SSA参数的选择对模型性能的影响,建议进行参数敏感性分析以获得最佳效果。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值