题目链接:bzoj2956
题目大意:
求
∑ni=1∑mj=1(nmodi)×(mmodj),i≠j
题解:
数论+分块
看清题啊
i≠j
!
所以式子就变成:
∑i=1n∑j=1m(nmodi)×(mmodj)−∑i=1min(n,m)(nmodi)×(mmodj)
像 bzoj1257那样处理
于是就变成了:
∑i=1n∑j=1m(n−i×⌊ni⌋)×(m−j×⌊mj⌋)−∑i=1min(n,m)(n−i×⌊ni⌋)×(m−j×⌊mj⌋)
傻傻的我把两边都拆开直接搞了。
于是对于 ∑ni=1∑mj=1i×j×⌊ni⌋×⌊mj⌋ 表示打了个 O(n√m−−√) 的…诶好蠢啊ww
应该把式子化成:
∑i=1n(n−i×⌊ni⌋)∑j=1m(m−j×⌊mj⌋)−∑i=1min(n,m)(n−i×⌊ni⌋)×(m−j×⌊mj⌋))
这样前面的就可以把 n 、
总结:
有时可以考虑内外层求和对调啊
还有分块的套路要熟悉。要有意识看到
⌊ni⌋
、
⌊ni⌋×⌊mj⌋
之后用分块啊
话说,这题mod得我真恶心,因为mod得不够多WA了两三次!
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
const LL ny=3323403;
const LL mod=19940417;
LL mymin(LL x,LL y){return (x<y)?x:y;}
LL get(LL n,LL lim)
{
LL ret=0,r;
for (LL i=1;i<=lim;i=r+1)
{
LL x=n/i;r=mymin(lim,n/x);
ret=(ret+x%mod*((r-i+1)*(r+i)/2%mod))%mod;
}return ret;
}
int main()
{
//freopen("a.in","r",stdin);
//freopen("a.out","w",stdout);
LL n,m,ans;
scanf("%lld%lld",&n,&m);
ans=(((n*n)%mod-get(n,n)+mod)%mod)*(((m*m)%mod-get(m,m)+mod)%mod)%mod;
LL lim=mymin(n,m),x,y,r,lx,ly;
ans=(ans-(lim*n%mod*m%mod)+mod)%mod;
for (LL i=1;i<=lim;i=r+1)
{
x=n/i;y=m/i;
r=mymin(n/x,m/y);
ans=(ans+m*x%mod*((r-i+1)*(r+i)/2%mod))%mod;
ans=(ans+n*y%mod*((r-i+1)*(r+i)/2%mod))%mod;
lx=((r*(r+1)%mod*(2*r+1)%mod)*ny)%mod;
ly=(((i-1)*i%mod*(2*i-1)%mod)*ny)%mod;
ans=(ans-(x*y%mod*((lx-ly+mod)%mod))%mod)%mod;
ans=(ans+mod)%mod;
}
while (ans<0) ans+=mod;
printf("%lld\n",ans%mod);
return 0;
}