首页 > 其他分享 >[国家集训队] Crash的数字表格 / JZPTAB

[国家集训队] Crash的数字表格 / JZPTAB

时间:2024-05-06 21:24:09浏览次数:29  
标签:frac gcd min sum JZPTAB ij Crash 国家集训队 dp

题目所求即

\[\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{gcd(i,j)} \]

这里没有出现\([gcd(x,y)=1]\),所以我们枚举\(gcd\)的值来硬凑,原式就等于

\[\sum_{d=1}^{min(n,m)}\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{gcd(i,j)}[gcd(i,j)=d] \]

为了出现\([gcd(i,j)=1]\),直接将\(i,j\)变成\(d\)的倍数,原式就等于

\[\sum_{d=1}^{min(n,m)}\sum_{d|i}\sum_{d|j}\frac{ij}{d}[gcd(i,j)=d]=\sum_{d=1}^{min(n,m)}\sum_{k=1}^{\frac{n}{d}}\sum_{t=1}^{\frac{m}{d}}ktd[gcd(k,t)=1]=\sum_{d=1}^{min(n,m)}d\sum_{k=1}^{\frac{n}{d}}\sum_{t=1}^{\frac{m}{d}}kt\sum_{p|gcd(k,t)}u(p)=\sum_{d=1}^{min(n,m)}d\sum_{p=1}^{min(\frac{n}{d},\frac{m}{d})}u(p)\sum_{p|k}\sum_{p|t}kt=\sum_{d=1}^{min(n,m)}d\sum_{p=1}^{min(\frac{n}{d},\frac{m}{d})}u(p)\sum_{i=1}^{\frac{n}{dp}}ip\sum_{j=1}^{\frac{m}{dp}}jp=\sum_{d=1}^{min(n,m)}d\sum_{p=1}^{min(\frac{n}{d},\frac{m}{d})}u(p)p^2(\sum_{i=1}^{\frac{n}{dp}}i)(\sum_{j=1}^{\frac{m}{dp}}j) \]

最后两项用等差数列求和公式就好了,然后分块套分块即可

标签:frac,gcd,min,sum,JZPTAB,ij,Crash,国家集训队,dp
From: https://www.cnblogs.com/dingxingdi/p/18175963

相关文章

  • [国家集训队] 矩阵乘法 题解
    发现实际上就是二维静态区间最大值,可以用整体二分维护。时间复杂度\(O((q+n^2)\log\max(a_{i,j})\logn^2)\)。#include<bits/stdc++.h>#definelllonglongusingnamespacestd;constintW=310005;constintQ=6e4+5;intn,q,w,ans[Q];intc[505][505],m;voidadd(i......
  • P2839 [国家集训队] middle
    Sol:首先注意到答案是具有单调性的,考虑二分答案\(x\)解决。令$up(l,r,x)/down(l,r,x)$是\([l,r]\)中大于等于/小于\(x\)的数。那么对于一个区间\([l,r]\),显然中位数\(\gex\)的条件为\(up(l,r,x)\gedown(l,r,x)\).变形得到\(up(l,r,x)-down(l,r,......
  • 洛谷 P4451 [国家集训队] 整数的lqp拆分
    洛谷传送门设\(G\)为斐波那契数列的生成函数,显然有\(F=F\timesG+1\)。那么\(F=\frac{1}{1-G}=1+\frac{x}{1-2x-x^2}\)。问题是如何展开\(\frac{x}{1-2x-x^2}\)。因为\(\frac{1}{1-ax}=\sum\limits_{i\ge0}(ax)^i\),所以考虑设\(\frac{x}{1-......
  • C113 带修莫队 P1903 [国家集训队] 数颜色/维护队列
    视频链接:   LuoguP1903[国家集训队]数颜色/维护队列//带修莫队O(n^(5/3))#include<iostream>#include<cstring>#include<algorithm>#include<cmath>usingnamespacestd;constintN=1000005;intn,m,B,mq,mr,a[N];intsum,cnt[N],ans[N];st......
  • C112 莫队算法 P1494 [国家集训队] 小 Z 的袜子
    视频链接:  LuoguP1494[国家集训队]小Z的袜子//普通莫队O(n*sqrt(n))#include<iostream>#include<cstring>#include<algorithm>#include<cmath>usingnamespacestd;constintN=50005;intn,m,B,a[N];intsum,cnt[N],ans1[N],ans2[N];str......
  • 【QT+QGIS跨平台编译】之九十:【QGIS_Crashhandler+Qt跨平台编译】(一套代码、一套框架,
    文章目录一、QGIS_Crashhandler介绍二、QGIS下载三、文件分析四、pro文件五、编译实践一、QGIS_Crashhandler介绍  QGIS_Crashhandler模块是QGIS中的一个重要组成部分,它提供了QGIS程序的错误崩溃处理与跟踪。二、QGIS下载QGIS网址:QGISSourceDownload......
  • P1829 / SP8099 Crash的数字表格 题解
    P1829/SP8099Crash的数字表格题解莫比乌斯反演、数论分块、杜教筛。题目传送门题意简述求以下式子的值,对\(20101009\)(一个质数)取模:\[\sum\limits_{i=1}^n\sum\limits_{j=1}^m\operatorname{lcm}(i,j)\]\(n,m\le10^7\)。加强:\(n,m\le10^{10}\)。解法莫比乌斯反......
  • openGauss 由于RemoveIPC未关闭导致数据库crash
    openGauss由于RemoveIPC未关闭导致数据库crashsemop引发的数据库crash--主库FATAL:semop(id=xxxxx)failed:IdentifierremovedFATAL:semctl(xxxxxx,11,SETVAL,0)failed:Invalidargument--备库FATAL:semctl(xxxxxx,11,SETVAL,0)failed:InvalidargumentLOG......
  • P1975 [国家集训队] 排队 题解
    题目链接:排队水紫,\(n\)不大,树套树或者分块都能做。分块的话,最优序列分块套套值域分块最优。观察到是可差性问题维护,即权值数量维护,那我们就树状数组套权值线段树即可。由于\(n\)不大,我们可以不用回收标记,直接数组空间开大点就行。我们预处理出初始逆序对,每一次操作都是基于......
  • P4555 [国家集训队] 最长双回文串
    题目链接:https://www.luogu.com.cn/problem/P4555题解:要找一个由两个回文组成字符串,一定有一个分割的位置,将两个回文串分开,设分割x,x+1,可能成为最后答案的值一定是左边的最长回文串和右边的最长的回文串长度之和。   回文中心一个<x,一个>x+1且一定包含x和x+1。如果最......