【牛客练习赛53】:B, C, E 题解

B:美味果冻

分析:

更换一下枚举顺序得到:

\sum_{j=1}^{n}\sum_{i=j}^{n}i*\left \lfloor \frac{i}{j} \right \rfloor^j

容易发现 i 为 [j,2j),[2j,3j)... [(t-1)*j,t*j)时后面向下取整的结果是一样的,分块计算即可,j 次幂通过(j-1)次幂 O(1) 得到,时间复杂度为 O(NlogN)

代码:

#include <bits/stdc++.h>
 
using namespace std;
typedef long long LL;
const int mod = 1e9+7;
const int maxn = 3e6+36;
LL Pow[maxn];
inline LL qpow(int a,int x){
    if(x == 1) return Pow[a];
    return Pow[a] = Pow[a]*a%mod;
}
LL solve(int n){
    LL res = 0;
    for(int j = 1;j <= n; ++j){
        int t = n/j; LL tep = 0;
        for(int i = 1;i < t; ++i){
            tep += 1ll*j*(2ll*i*j+j-1)/2%mod*qpow(i,j)%mod;
            if(tep>=mod) tep -= mod;
        }
        tep += 1ll*(n-t*j+1)*(n+t*j)/2%mod*qpow(t,j)%mod;     //单独计算最后非完整的块
        res = (res + tep) % mod;
    }
    return res;
}
int main(){
    int n;  cin >> n;
    for(int i = 1;i<=n; ++i) Pow[i] = i;
    cout << solve(n) << '\n';
    return 0;
}

C:富豪凯匹配串

分析:

AC自动机是不可能的

数据范围是支持暴力枚举文本串的,那就枚举每一个文本串,只用把匹配串下划线对应的位置搞成 0,然后比较二者是否一样即可,字符串只包含0,1,bitset 优化上述两种操作即可

代码:

#include <bits/stdc++.h>

using namespace std;
const int maxn = 1e3+13;
bitset<maxn> b[maxn];
int n,m,Q;
char s[maxn];
int main(){
    scanf("%d %d",&n,&m);
    for(int i = 0;i < n; ++i){
        scanf("%s",s);
        for(int j = 0;j < m; ++j) if(s[j]=='1') b[i].set(j); 
    }
    scanf("%d",&Q);
    while(Q--){
        scanf("%s",s);
        bitset<maxn> bb,aa;
        for(int i = 0;i < m; ++i){
            if(s[i] != '_') bb.set(i);
            if(s[i] == '1') aa.set(i);
        }
        int cnt = 0;
        for(int i = 0;i < n; ++i){
            if((bb&b[i])==aa) cnt++;
        }
        printf("%d\n",cnt);
    }
    return 0;
}

 E:老瞎眼 pk 小鲜肉

分析:

区间异或值为 0,套路的算出每个位置的前缀异或和,只要两个位置的前缀异或和一样,那么这段区间的异或值就为 0;问题就转换为求区间内两个最近的点对,先将询问离线按右端点排序,依次移动扫描指针,求出当前位置的异或前缀和,只需要更新这个值上一次出现的位置的最小值即可,因为这两个点对一定是最近的(仅对于这个异或值作为R来说),那么扫描到询问直接查询区间最小值即可

代码:

#include <bits/stdc++.h>

using namespace std;
const int maxn = 5e5+55;
const int inf = 0x3f3f3f3f;
int n,m,a[maxn];
struct query{
    int l,r,id;
}q[maxn];
bool cmp(query a,query b){
    return a.r < b.r;
}
int tr[maxn<<2],res[maxn];
void build(int l,int r,int x){
    tr[x] = inf;
    if(l == r){ tr[x] = inf; return ;}
    int mid = (l+r)>>1;
    build(l,mid,x<<1);
    build(mid+1,r,x<<1|1);
}
void updata(int l,int r,int pos,int x,int v){
    if(l == r){tr[x] = v; return ;}
    int mid = (l+r) >> 1;
    if(pos > mid) updata(mid+1,r,pos,x<<1|1,v);
    else updata(l,mid,pos,x<<1,v);
    tr[x] = min(tr[x<<1],tr[x<<1|1]);
}
void query(int l,int r,int L,int R,int x,int &ans){
    if(l > R || r < L) return ;
    if(l >= L && r <= R){ans = min(ans,tr[x]); return ;}
    int mid = (l+r) >> 1;
    query(l,mid,L,R,x<<1,ans);
    query(mid+1,r,L,R,x<<1|1,ans);
}
map<int,int> mp;
int main(){
    scanf("%d %d",&n,&m);
    for(int i = 1;i <= n; ++i) scanf("%d",a+i);
    for(int i = 0;i < m; ++i) scanf("%d %d",&q[i].l,&q[i].r),q[i].id=i;
    sort(q,q+m,cmp); build(1,n,1);
    int R = 0,sum = 0; mp[0] = 0;
    for(int i = 0;i < m; ++i){
        while(R < q[i].r){
            sum = sum^a[++R];
            if(mp.count(sum)) updata(1,n,mp[sum]+1,1,R-mp[sum]);
            mp[sum] = R;
        }
        int ans = inf; query(1,n,q[i].l,q[i].r,1,ans);
        res[q[i].id] = ans > n ? -1 : ans;
    }
    for(int i = 0;i < m; ++i) printf("%d\n",res[i]);
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值