B:美味果冻
分析:
更换一下枚举顺序得到:
容易发现 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;
}
分析:
区间异或值为 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;
}