title : 暨南大学2023东软教育杯ACM校赛 题解
tags : ACM,练习记录
date : 2023-3-26
author : Linno
题目链接:https://ac.nowcoder.com/acm/contest/47948
出题数量:12/12 ,AK了
出题顺序:A->B->C->D->F->G->J->E->I->H->L->K
简评:首先感谢出题组手下留情,再来点中等水平思维题估计就能搞死我了。因为题出的都比较板所以顺理成章地AK了,打搅了大家的做题兴致非常抱歉!这一场如果板子带够了并且刷题量达到一定水平,打得都不会差的。如果觉得自己打得不好,那希望这篇题解能给你带来帮助(不懂也可以继续问),喜欢ACM的话就请继续加油吧!
选某一个数,它的因数都会被选中,因此选它本身肯定是最优的,那么我们直观觉得从大到小枚举的时候把因数都打上标签,然后没打标签的数统计到答案里面就必然是正确的。但可惜n有1e18这样做肯定超时,因此我们继续考虑更直观的结论:
void solve(){
int n;
cin>>n;
cout<<(n+1)/2<<"\n";
}
注意读题,苏神是可以第一的,但是成绩不能并列。因此一般情况下苏神的成绩是第二名成绩-1秒,但是如果那是第一名的成绩,就要跑得比这个更快一秒才行了。
void solve(){
cin>>n;
for(int i=1;i<=n;++i) cin>>a[i];
stable_sort(a+1,a+1+n);
if(a[1]+1==a[2]) ans=a[1]-1;
else ans=a[2]-1;
cout<<ans<<"\n";
}
前置知识:中国剩余定理
原题意也就是让我们解同余方程组:
{ x ≡ b 1 m o d a 1 x ≡ b 2 m o d a 2 . . . x ≡ b n m o d a n \begin{cases} x\equiv b_1 \mod a_1 \\ x\equiv b_2 \mod a_2 \\ ...\\ x\equiv b_n \mod a_n \\ \end{cases} ⎩ ⎨ ⎧x≡b1moda1x≡b2moda2...x≡bnmodan
由于 a a a两两互质,令 x = ( N / a i ) ∗ y x=(N/a_i)*y x=(N/ai)∗y,方程组等同于解同余方程 ( N / a i ) y ≡ 1 m o d a i (N/a_i)y\equiv 1 \mod a_i (N/ai)y≡1modai,得到特解 y i y_i yi,则方程组的解为: x 0 = b 1 x 1 + b 2 x 2 + . . . + b r x r m o d N x_0=b_1x_1+b_2x_2+...+b_rx_r \mod N x0=b1x1+b2x2+...+brxrmodN,在模N意义下唯一。
注意到准备的箱子数 a i a_i ai肯定是不同的质数,也就是说不需要用到扩展CRT,直接套板子即可。
#include
#define int long long
using namespace std;
const int N=20;
int n,a[N],b[N],m[N],t[N],M=1;
inline int exgcd(int a,int b,int &x,int &y){
if(!b){
x=1;y=0;
return a;
}
int d=exgcd(b,a%b,y,x);
y-=a/b*x;
return d;
}
inline int inv(int a,int b){
int d,x,y;
d=exgcd(a,b,x,y);
return (x<0)?(x+b):x;
}
void solve(){
cin>>n;
for(int i=1;i<=n;++i){
cin>>a[i]>>b[i];
M*=a[i];
}
int ans=0;
for(int i=1;i<=n;++i){
m[i]=M/a[i];
t[i]=inv(m[i],a[i]);
ans+=b[i]*m[i]%M*t[i]%M;
ans%=M;
}
cout<<ans<<"\n";
}
signed main(){
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
int T=1;
while(T--){
solve();
}
return 0;
}
题目给定 x , b x,b x,b,我们可以转化为每个节点的做右端点 [ l , r ] [l,r] [l,r],并且按照 r r r从小到大排序,每次贪心地选最右点建基站并跳过已被覆盖掉的节点,即可保证最优。
#include
#define int long long
using namespace std;
const int N=1e5+7;
struct E{
int l,r;
bool operator <(const E &B)const{
return r<B.r;
}
}s[N];
void solve(){
int n;
cin>>n;
for(int i=1,x,b;i<=n;++i){
cin>>x>>b;
s[i].l=x-b;s[i].r=x+b;
}
stable_sort(s+1,s+1+n);
int lst=-0x3f3f3f3f,ans=0;
for(int i=1;i<=n;++i){
if(lst<s[i].l){
lst=s[i].r;
++ans;
}
}
cout<<ans<<"\n";
}
signed main(){
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
int T=1;
while(T--){
solve();
}
return 0;
}
前置知识:数位DP
数位DP板题,形式一般都如下: 给定 L , R , 问 [ L , R ] 范围内有多少 X X X 的数,数据范围可出到 1 e 18 给定L,R,问[L,R]范围内有多少XXX的数,数据范围可出到1e18 给定L,R,问[L,R]范围内有多少XXX的数,数据范围可出到1e18。
流程大概就是按位拆分给的数,然后再按位记忆化搜索统计答案。代码细节各位自己看吧,dfs(stp,zero,lim,mx,now)表示搜到第stp位数,最多连续mx位6,当前连续了now位6时的状态,zero是前置0标签,lim是最高位标签。
#include
#define int long long
using namespace std;
const int N=66;
int len=0,num[N],vis[N][N][N];
int dfs(int stp,bool zero,bool lim,int mx,int now){
if(!stp) return (mx>=3);
if(!zero&&!lim&&vis[stp][mx][now]) return vis[stp][mx][now];
int j=lim?num[stp]:9,ans=0;
for(int i=0;i<=j;++i){
int tmp=(i==6)?now+1:0;
ans+=dfs(stp-1,zero&(i==0),lim&(i==num[stp]),max(mx,tmp),tmp);
}
if(!zero&&!lim) vis[stp][mx][now]=ans;
return ans;
}
int solve(int x){
len=0;
memset(vis,0,sizeof(vis));
while(x){
num[++len]=x%10;
x/=10;
}
return dfs(len,1,1,0,0);
}
signed main(){
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
int L,R;
cin>>L>>R;
cout<<solve(R)-solve(L-1)<<"\n";
return 0;
}
/*
1 1000000000000000000
*/
看大家的代码写的都好短,我开数组记了前缀和和后缀和。然后枚举‘z’的位置,每个位置对答案的贡献就是前面‘j’的数量乘上后面‘h’的数量,记得取模。
#include
#define int long long
using namespace std;
const int mod=998244353;
const int N=2e6+7;
int len;
string str;
int numj[N],numh[N];
void solve(){
cin>>len;
cin>>str;
int ans=0;
numj[0]=(str[0]=='j');numh[len]=0;
for(int i=1;i<len;++i){
numj[i]=numj[i-1]+(str[i]=='j');
}
for(int i=len-1;i>=0;--i){
numh[i]=numh[i+1]+(str[i]=='h');
}
// for(int i=0;i
for(int i=1;i<len;++i){
if(str[i]=='z') ans+=numj[i-1]*numh[i+1]%mod,ans%=mod;
// cout<
}
cout<<ans<<"\n";
}
signed main(){
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
int T=1;
while(T--){
solve();
}
return 0;
}
/*
14
mjbajzhmuaxing
9
jzhjzhjzh
*/
爆搜就行了,首先把最外面一圈全变成1,然后剩下里面的就直接全变成6,最后统计答案即可。
#include
#define int long long
using namespace std;
const int N=55;
int n,vis[N][N],mp[N][N];
int xx[]={0,0,-1,1},yy[]={-1,1,0,0};
bool check(int x,int y){return (x>=1&&x<=n&&y>=1&&y<=n&&!vis[x][y]&&mp[x][y]==0);}
inline void dfs(int x,int y,int id){
mp[x][y]=id;vis[x][y]=1;
for(int d=0;d<4;++d){
int nx=x+xx[d],ny=y+yy[d];
if(check(nx,ny)) dfs(nx,ny,id);
}
}
void solve(){
cin>>n;
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j) cin>>mp[i][j];
}
for(int i=1;i<=n;++i){
if(!vis[i][1]&&mp[i][1]==0) dfs(i,1,1);
if(!vis[i][n]&&mp[i][n]==0) dfs(i,n,1);
if(!vis[n][i]&&mp[n][i]==0) dfs(n,i,1);
if(!vis[1][i]&&mp[1][i]==0) dfs(1,i,1);
}
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j){
if(!vis[i][j]&&mp[i][j]==0) dfs(i,j,6);
}
}
int ans=0;
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j){
if(mp[i][j]==6) ++ans;
}
}
cout<<ans<<"\n";
}
signed main(){
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
int T=1;
while(T--){
solve();
}
return 0;
}
没想到这题居然是暴力……大家可以忽略我的做法去看别人的。题目如果不保证所有i累加小于1e7的话,可以考虑建平衡树,我用的是FHQ Treap。这种数据结构可以很方便地进行序列的拆分和合并操作。然后对于拆分之后的两颗平衡树每个节点都有加法标签和乘法标签,仿照线段树的pushdown操作进行先乘后加的处理即可。可能场内没有想我一样写平衡树的同学吧……pushdown真的很坑,让我wa了两发。
#include
#define int long long
using namespace std;
const int N=2e5+7;
const int mod=1000000007;
inline int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){
if(ch=='-') f=f*-1;
ch=getchar();
}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x*f;
}
int n,q,a[N];
int ch[N][2],val[N],pri[N],sz[N],tg1[N],tg2[N],cnt,rt;
void update(int x){
sz[x]=1+sz[ch[x][0]]+sz[ch[x][1]];
}
#define mid ((l+r)>>1)
#define ls (ch[x][0])
#define rs (ch[x][1])
void pushdown(int x){
if(x&&tg2[x]!=1){
if(ls){
val[ls]=val[ls]*tg2[x]%mod;
tg1[ls]=tg1[ls]*tg2[x]%mod;
tg2[ls]=tg2[ls]*tg2[x]%mod;
}
if(rs){
val[rs]=val[rs]*tg2[x]%mod;
tg1[rs]=tg1[rs]*tg2[x]%mod;
tg2[rs]=tg2[rs]*tg2[x]%mod;
}
tg2[x]=1;
}
if(x&&tg1[x]){
if(ls) val[ls]=(val[ls]+tg1[x])%mod,tg1[ls]=(tg1[ls]+tg1[x])%mod;
if(rs) val[rs]=(val[rs]+tg1[x])%mod,tg1[rs]=(tg1[rs]+tg1[x])%mod;
tg1[x]=0;
}
}
int new_node(int v){
sz[++cnt]=1;
tg1[cnt]=0;tg2[cnt]=1;
val[cnt]=v;
pri[cnt]=rand();
return cnt;
}
int merge(int x,int y){
if(!x||!y) return x+y;
pushdown(x);pushdown(y);
if(pri[x]<pri[y]){
ch[x][1]=merge(ch[x][1],y);
update(x);
return x;
}else{
ch[y][0]=merge(x,ch[y][0]);
update(y);
return y;
}
}
void split(int now,int k,int &x,int &y){
if(!now) x=y=0;
else{
pushdown(now);
if(k<=sz[ch[now][0]]) y=now,split(ch[now][0],k,x,ch[now][0]);
else x=now,split(ch[now][1],k-sz[ch[now][0]]-1,ch[now][1],y);
update(now);
}
}
int build(int l,int r){
if(l>r) return 0;
int now=new_node(a[mid]);
// cout<
ch[now][0]=build(l,mid-1);
ch[now][1]=build(mid+1,r);
update(now);
return now;
}
void dfs(int x){
if(!x) return;
pushdown(x);
dfs(ls);
printf("%lld ",val[x]%mod);
dfs(rs);
}
signed main(){
srand(time(0));
n=read();
for(int i=1;i<=n;++i) a[i]=read()%mod;
rt=build(1,n);
q=read();
for(int i=1,x,op,y,a,b;i<=q;++i){
x=read();op=read();y=read();
split(rt,x,a,b);
if(op==1) val[a]=(val[a]+y)%mod,tg1[a]=(tg1[a]+y)%mod;
else val[a]=val[a]*y%mod,tg2[a]=tg2[a]*y%mod;
rt=merge(b,a);
}
dfs(rt);
puts("");
}
/*
7
1 2 3 4 5 6 7
3
3 1 7
2 1 1
1 2 2
4
2 5 4 3
0
4
2 5 4 3
3
1 2 2
3 1 3
1 2 5
——8 7 6 20
3
0 5 4
3
2 2 5
2 2 3
3 1 7
*/
前置知识:线段树
首先要想到一个很直观的结论就是:我们每次都只选择长度为2的区间可以保证答案是最大的。因此我们只需要维护相邻两项的乘积即可。题目转化为了区间查询最大值和单点修改,学过数据结构的同学可以快速通过。
#include
#define int long long
using namespace std;
const int N=2e5+7;
#define ls (p<<1)
#define rs (p<<1|1)
#define mid ((l+r)>>1)
int n,k,a[N],b[N],mx[N<<2];
void build(int p,int l,int r){
if(l==r){
mx[p]=b[l];
return;
}
build(ls,l,mid);
build(rs,mid+1,r);
mx[p]=max(mx[ls],mx[rs]);
}
void update(int p,int l,int r,int ql,int qr,int k){
if(ql<=l&&r<=qr){
mx[p]=k;
return;
}
if(ql<=mid) update(ls,l,mid,ql,qr,k);
if(qr>mid) update(rs,mid+1,r,ql,qr,k);
mx[p]=max(mx[ls],mx[rs]);
}
int query(int p,int l,int r,int ql,int qr){
if(ql<=l&&r<=qr) return mx[p];
int res=0;
if(ql<=mid) res=query(ls,l,mid,ql,qr);
if(qr>mid) res=max(res,query(rs,mid+1,r,ql,qr));
return res;
}
void solve(){
cin>>n;
for(int i=1;i<=n;++i) cin>>a[i];
for(int i=1;i<n;++i) b[i]=a[i]*a[i+1];b[n]=b[n-1];
build(1,1,n);
cin>>k;
for(int i=1,op,l,r;i<=k;++i){
cin>>op>>l>>r;
if(op==1){
b[l]=b[l]/a[l]*r;
if(l>1) b[l-1]=b[l-1]/a[l]*r;
if(l==n) b[n-1]=b[n];
else if(l==n-1) b[n]=b[n-1];
a[l]=r;
update(1,1,n,n-1,n-1,b[n-1]);
update(1,1,n,n,n,b[n]);
update(1,1,n,l,l,b[l]);
if(l>1) update(1,1,n,l-1,l-1,b[l-1]);
}else{
cout<<query(1,1,n,l,r-1)<<"\n";
}
}
}
signed main(){
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
int T=1;
while(T--){
solve();
}
return 0;
}
/*
5
1 2 3 4 5
4
2 1 2
2 1 3
2 4 5
2 3 5
5
1 2 5 3 4
4
2 1 2
2 1 3
2 4 5
2 3 5
6
2 8 9 1 11 3
6
2 1 5
2 1 3
1 1 10
2 1 2
2 2 3
2 1 3
6
2 8 9 1 11 3
7
2 1 5
2 1 3
1 6 10
2 1 6
2 5 6
2 4 5
2 4 6
*/
前置知识:图论基础知识,最小费用最大流
最大流是在一个有向图中从S到T,中间有许多节点,每条边权限制了流量,问单位时间能流过的最大流量是多少的模型。而最小费用最大流流是在这基础上给出每条边每流过一单位水流的费用,问在最大流量的基础上最小费用的模型。显然题目可以套入这个模型,从源点S到每个地点之间的最大流量便是单车最开始的摆放量 a i a_i ai,不造成费用;每两个地点之间的流量无限,即可以不限量的搬运,但是需要造成的费用便是两地的距离 d i s i j dis_{ij} disij;每个地点到汇点的流量便是最终要摆放的车辆 b i b_i bi,也不造成费用。这样,在满足摆放好车辆的前提下(最大流),所需要的努力最少(最小费用)。模型是怎么做到这一点的?去学前置知识。
#include
#define int long long
#define inf 0x7f7f7f7f
using namespace std;
const int N=305,M=2e5+7;
inline int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){
if(ch=='-') f=f*-1;
ch=getchar();
}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x*f;
}
inline void write(int x){
if(x>9) write(x/10);
putchar(x%10+'0');
}
int v[M],nxt[M],f[M],cnt=1;
int flow[N],w[M],head[N];
inline void addedge(int x,int y,int z,int k){
++cnt;v[cnt]=y;w[cnt]=z;f[cnt]=k;nxt[cnt]=head[x];head[x]=cnt;
++cnt;v[cnt]=x;w[cnt]=0;f[cnt]=-k;nxt[cnt]=head[y];head[y]=cnt;
}
int n,m,S,T;
int dis[N],mcost,mflow;
int inq[N],pre[N];
int mp[N][N],tot=0,q[M*10];
inline bool spfa(){
for(int i=S;i<=T;++i) dis[i]=inf,inq[i]=0;
int L=1,R=1;q[L]=S;
inq[S]=1;dis[S]=0;flow[S]=inf;
while(L<=R){
int fro=q[L];++L;
inq[fro]=0;
for(int i=head[fro];i;i=nxt[i]){
if(w[i]&&dis[v[i]]>dis[fro]+f[i]){
dis[v[i]]=dis[fro]+f[i];
flow[v[i]]=min(flow[fro],w[i]);
pre[v[i]]=i;
if(!inq[v[i]]){
q[++R]=v[i];
inq[v[i]]=1;
}
}
}
}
return dis[T]!=inf;
}
inline void update(){
int x=T;
while(x!=S){
int i=pre[x];
w[i]-=flow[T];
w[i^1]+=flow[T];
x=v[i^1];
}
mflow+=flow[T];
mcost+=flow[T]*dis[T];
}
void EK(){
while(spfa()) update();
}
void solve(){
n=read();
S=0;T=n+1;
for(int i=1;i<=n;++i) addedge(S,i,read(),0);
for(int i=1;i<=n;++i) addedge(i,T,read(),0);
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j){
mp[i][j]=read();
addedge(i,j,inf,mp[i][j]);
}
}
EK();
write(mcost);
}
signed main(){
int T=1;
while(T--){
solve();
}
return 0;
}
前置知识:SAM(后缀自动机)
因为被选中的位置不可用,所对于两个字符串S,T来说都可以分段处理。也就是说S的每一段与T的每一段两两匹配,求最长公共子串的最大长度。显然两两匹配会超时,我们不妨对S的每一段字符串建广义SAM,然后枚举T的每一段在SAM上跑最长匹配长度。别问我SAM是怎么做到的(其实就是在DAG上一个一个字符去匹配,如果匹配失败就像KMP一样跳转到Fail的位置),去学前置知识。听说大佬用后缀数组被卡了,可惜。
#include
using namespace std;
const int N=2e6+7;
int a[N],b[N],lst,tot=1,rt=1,num,ans=0;
int fa[N],len[N],ch[N][30],sz[N][12];
char s[N];
vector<int>G[N];
vector<string>s1,s2;
string str1,str2;
void insert(int c,int id){
int p=lst,np=lst=++tot;
len[np]=len[p]+1,sz[np][id]++;
while(p&&!ch[p][c]) ch[p][c]=np,p=fa[p];
if(!p) fa[np]=rt;
else{
int q=ch[p][c];
if(len[q]==len[p]+1) fa[np]=q;
else{
int nq=++tot;
memcpy(ch[nq],ch[q],sizeof(ch[q]));
len[nq]=len[p]+1,fa[nq]=fa[q],fa[q]=fa[np]=nq;
while(ch[p][c]==q) ch[p][c]=nq,p=fa[p];
}
}
}
void dfs(int x){
int len=G[x].size();
for(int i=0;i<len;++i){
int y=G[x][i];
dfs(y);
for(int id=1;id<=num;++id) sz[x][id]+=sz[y][id];
}
}
bool check(int p){
if(!p) return false;
for(int i=1;i<=num;++i){
if(sz[p][i]) return true;
}
return false;
}
void work(string s){
int p=rt,l=0;
for(int i=0;i<s.length();++i){
int c=s[i]-'a';
if(check(ch[p][c])) l++,p=ch[p][c];
else{
while(p&&!check(ch[p][c])) p=fa[p];
if(p) l=len[p]+1,p=ch[p][c];
else l=0,p=rt;
}
ans=max(ans,l);
}
}
signed main(){
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
cin>>str1>>str2;
int l1=str1.length(),l2=str2.length(),n,m;
memset(a,-1,sizeof(a));
memset(b,-1,sizeof(b));
cin>>n;for(int i=1;i<=n;++i) cin>>a[i];
cin>>m;for(int i=1;i<=m;++i) cin>>b[i];
stable_sort(a+1,a+1+n);
stable_sort(b+1,b+1+m);
//
string tmp="";
for(int i=0,j=1;i<l1;++i){
if(i==a[j]){
if(tmp!="") s1.emplace_back(tmp);
++j;tmp="";
}else tmp.push_back(str1[i]);
}
if(tmp!="") s1.emplace_back(tmp);
tmp="";
for(int i=0,j=1;i<l2;++i){
if(i==b[j]){
if(tmp!="") s2.emplace_back(tmp);
++j;tmp="";
}else tmp.push_back(str2[i]);
}
if(tmp!="") s2.emplace_back(tmp);
//
for(auto s:s1){
lst=rt;num++;
// cout<
for(int i=0;i<s.length();++i) insert(s[i]-'a',num);
}
for(int i=2;i<=tot;++i) G[fa[i]].emplace_back(i);
dfs(rt);
for(auto s:s2){
work(s);
// cout<
}
cout<<ans<<"\n";
return 0;
}
/*
aabaabaab
aa
3
3 6 9
0
abcdabcdbacbd
bcd
1
4
1
1
*/
前置知识:数论分块,迪利克雷卷积,莫比乌斯反演
下面是推导过程,如果有不懂的私信我或者去看前置知识。
题意转化为求 ∑ i = 1 S 1 ∑ j = 1 S 2 l c m ( i , j ) \sum_{i=1}^{S1}\sum_{j=1}^{S2}lcm(i,j) ∑i=1S1∑j=1S2lcm(i,j),最小公倍数 l c m ( i , j ) = i ∗ j g c d ( i , j ) lcm(i,j)=\frac{i*j}{gcd(i,j)} lcm(i,j)=gcd(i,j)i∗j
原式等于 ∑ i = 1 n ∑ j = 1 m i ⋅ j gcd ( i , j ) = ∑ d = 1 n d ⋅ ∑ i = 1 ⌊ n d ⌋ ∑ j = 1 ⌊ m d ⌋ [ gcd ( i , j ) = 1 ] i ⋅ j = ∑ d = 1 n ∑ d ∣ i n ∑ d ∣ j m μ ( d ) ⋅ i ⋅ j 令 g ( n , m ) = ∑ i = 1 n ∑ j = 1 m i ⋅ j = n ⋅ ( n + 1 ) 2 × m ⋅ ( m + 1 ) 2 sum ( n , m ) = ∑ d = 1 n μ ( d ) ⋅ d 2 ⋅ g ( ⌊ n d ⌋ , ⌊ m d ⌋ ) 原式 = ∑ d = 1 n d ⋅ sum ( ⌊ n d ⌋ , ⌊ m d ⌋ ) , 可用数论分块和线性筛解决。 原式等于 \sum_{i=1}^n\sum_{j=1}^m\frac{i\cdot j}{\gcd(i,j)} \\ = \sum_{d=1}^n d\cdot\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}[\gcd(i,j)=1]\ i\cdot j \\ = \sum_{d=1}^n\sum_{d\mid i}^n\sum_{d\mid j}^m\mu(d)\cdot i\cdot j\\ 令 g(n,m)=\sum_{i=1}^n\sum_{j=1}^m i\cdot j=\frac{n\cdot(n+1)}{2}\times\frac{m\cdot(m+1)}{2} \\ \operatorname{sum}(n,m)=\sum_{d=1}^n\mu(d)\cdot d^2\cdot g(\lfloor\frac{n}{d}\rfloor,\lfloor\frac{m}{d}\rfloor) \\ 原式=\sum_{d=1}^n d\cdot\operatorname{sum}(\lfloor\frac{n}{d}\rfloor,\lfloor\frac{m}{d}\rfloor) ,可用数论分块和线性筛解决。 原式等于i=1∑nj=1∑mgcd(i,j)i⋅j=d=1∑nd⋅i=1∑⌊dn⌋j=1∑⌊dm⌋[gcd(i,j)=1] i⋅j=d=1∑nd∣i∑nd∣j∑mμ(d)⋅i⋅j令g(n,m)=i=1∑nj=1∑mi⋅j=2n⋅(n+1)×2m⋅(m+1)sum(n,m)=d=1∑nμ(d)⋅d2⋅g(⌊dn⌋,⌊dm⌋)原式=d=1∑nd⋅sum(⌊dn⌋,⌊dm⌋),可用数论分块和线性筛解决。
#include
#define int long long
using namespace std;
const int mod=1e9+7;
const int N=1e7+7;
int np[N],pri[N],mu[N],sum[N],cnt=0;
void init(){
mu[1]=np[1]=1;
for(int i=2;i<N;++i){
if(!np[i]) pri[++cnt]=i,mu[i]=-1;
for(int j=1;j<=cnt&&pri[j]*i<N;++j){
np[pri[j]*i]=1;
if(i%pri[j]) mu[i*pri[j]]=-mu[i];
else break;
}
}
for(int i=1;i<N;++i) sum[i]=(sum[i-1]+i*i%mod*(mu[i]+mod)%mod)%mod;
}
int Sum(int x,int y){
return (x*(x+1)/2%mod)*(y*(y+1)/2%mod)%mod;
}
int func(int x,int y){
int res=0;
for(int l=1,r;l<=min(x,y);l=r+1){
r=min(x/(x/l),y/(y/l));
res=(res+(sum[r]-sum[l-1]+mod)*Sum(x/l,y/l)%mod)%mod;
}
return res;
}
int Solve(int x,int y){
int res=0;
for(int l=1,r;l<=min(x,y);l=r+1){
r=min(x/(x/l),y/(y/l));
res=(res+(r-l+1)*(l+r)/2%mod*func(x/l,y/l)%mod)%mod;
}
return res;
}
signed main(){
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
init();
int t,n,m;
cin>>t;
while(t--){
cin>>n>>m;
cout<<Solve(n,m)<<"\n";
}
return 0;
}
/*
2
4 5
4 5
*/