首先对于答案

ΣΣ(n mod i)*(m mod j) i<>j

也就是Σ(n mod i)Σ(m mod j)-Σ(n mod i)(m mod i)

将mod展开，我们可以得到有floor的式子，对于这种式子，我们可以

利用分段的思想，将O(N)的简化为sqrt(n)的

/**************************************************************    Problem: 2956    User: BLADEVIL    Language: Pascal    Result: Accepted    Time:388 ms    Memory:224 kb****************************************************************/ //By BLADEVILconst    d39                             =19940417;     var    n, m                            :int64;    ans, ans2                       :int64; function min(a,b:int64):int64;begin    if a>b then min:=b else min:=a;end;     function calc(x,y:int64):int64;var    i, j                            :int64;    z                               :int64;begin    calc:=0;    i:=1;    while i<=y do    begin        j:=x div (x div i);        if j>y then j:=y;        z:=((i+j)*(j-i+1) div 2) mod d39;        calc:=(calc+(z*(x div i) mod d39) mod d39)mod d39;        i:=j+1;    end;end; function sum(x:int64):int64;var    a, b, c                         :int64;begin    if x=0 then exit(0);    a:=x; b:=x+1; c:=2*x+1;    if a mod 3=0 then a:=a div 3 else    if b mod 3=0 then b:=b div 3 else    if c mod 3=0 then c:=c div 3;    if a mod 2=0 then a:=a div 2 else    if b mod 2=0 then b:=b div 2 else    if c mod 2=0 then c:=c div 2;    sum:=a mod d39;    sum:=sum*b mod d39;    sum:=sum*c mod d39;end;  function fuck:int64;var    i, j                            :int64;    t1, t2                          :int64;    z                               :int64;begin    i:=1;    fuck:=0;    while i<=min(n,m) do    begin        t1:=n div (n div i);        t2:=m div (m div i);        j:=min(t1,t2);        z:=(((sum(j)-sum(i-1)) mod d39+d39) mod d39);        z:=(z*(n div i)) mod d39;        z:=(z*(m div i)) mod d39;        fuck:=(fuck+z) mod d39;        i:=j+1;    end;end; begin    read(n,m);    ans2:=calc(m,m) mod d39;    ans2:=((m*m-ans2) mod d39+d39) mod d39;    ans:=((n*n-calc(n,n)) mod d39*ans2) mod d39;    ans2:=(n*m mod d39)*min(n,m) mod d39;    ans2:=(ans2+fuck) mod d39;    ans2:=((ans2-m*calc(n,min(n,m)))mod d39+d39) mod d39;    ans2:=((ans2-n*calc(m,min(n,m)))mod d39+d39) mod d39;    ans:=((ans-ans2) mod d39+d39) mod d39;    writeln(ans);end.