P5495 Dirichlet 前缀和
题目描述
给定一个长度为 nnn 的数列 a1,a2,a3,…,ana_1,a_2,a_3,\dots,a_na1,a2,a3,…,an。
现在你要求出一个长度为 nnn 的数列
b1,b2,b3,…,bnb_1,b_2,b_3,\dots,b_nb1,b2,b3,…,bn,满足bk=∑i∣kaib_k=\sum_{i|k}a_i bk=i∣k∑ai
由于某些神秘原因,这里的 bkb_kbk 要对 2322^{32}232 取模。
输入格式
为了避免过大的输入,本题的输入使用随机数生成器。输入中只有一行两个整数 n,seedn,seedn,seed。其中 seedseedseed 为 323232 位无符号整数,用来生成数据。接下来,你要调用 nnn 次随机数生成器,分别生成 a1∼ana_1\sim a_na1∼an。
对于C/C++选手,生成器模板如下:
#define uint unsigned int
uint seed;
inline uint getnext(){
seed^=seed<<13;
seed^=seed>>17;
seed^=seed<<5;
return seed;
}
注意:所有 nnn 个数均为 323232 位无符号整数。
输出格式
为了避免过大的输出,你只需输出一个 323232 位无符号整数,表示所有 bib_ibi 的异或和。
输入
5 1477
输出
2608816472
题解
#include<bits/stdc++.h>
#define uint unsigned int
using namespace std;
uint seed;
inline uint getnext(){
seed^=seed<<13;
seed^=seed>>17;
seed^=seed<<5;
return seed;
}
uint num[30000000],t[30000000],prim[30000000],pri[30000000],cnt,ans;
int main()
{
uint n;
cin>>n>>seed;
for(uint i=2;i<=n;i++)
if(!prim[i])
for(uint j=2*i;j<=n;j+=i)
prim[j]=1;
for(uint i=2;i<=n;i++)
if(!prim[i]) pri[++cnt]=i;
for(uint i=1;i<=n;i++)
num[i]=getnext(),t[i]=num[i];
for(uint i=1;i<=cnt;i++){
for(uint j=1;pri[i]*j<=n;j++){
t[j*pri[i]]+=t[j];
}
}
for(uint i=1;i<=n;i++){
ans^=t[i];
}
cout<<ans;
}
其中,inline用于声明内联函数,内联函数作用可以省去函数调用的开销,编译器会将内联函数复制展开。
- inline定义与性质
定义
当函数被声明为内联函数之后,编译器会将其内联展开,而不是按通常函数调用。
性质
只有当函数只有10行甚至更少时才将其定义为内联函数
内联以代码膨胀(拷贝)为代价,仅仅省去了函数调用的开销。
函数调用开销:参数压栈、跳转、退栈和返回等操作,不包括执行函数体的开销关键字 inline 必须与函数定义体放在一起才能使函数成为内联,仅将 inline 放在函数声明前面不起任何作用。
// 如下风格的函数 Foo 不能成为内联函数:
inline void Foo(int x, int y); // inline 仅与函数声明放在一起
void Foo(int x, int y)
{
...
}
// 而如下风格的函数 Foo 则成为内联函数:
void Foo(int x, int y);
inline void Foo(int x, int y) // inline 与函数定义体放在一起
{
...
}