挖掘机技术哪家强

题目描述
有人问现实中为什么总是男生追求女生，反过来很少。实际上女生也是想主动追求男生的，但是世俗中对于主动追求男生的女生有种歧视，这样就使得女生不大敢主动追求男生。但是面对喜欢的男生，难道就不出手么？女生只能步步为营，挖坑来引诱男生往里跳。这时候问题就来了，挖掘机技术到底哪家强？

被热血沸腾的广告洗脑了若干天后，Matt终于下定决心，毅然登上了开往泉城的列车，决心寻找生活的希望。

来到布鲁谢特学院后，Matt逐渐地了解了各种型号的挖掘机。在这里我们可以认为有大挖掘机和小挖掘机两种。

今天Matt的任务很简单：首先他要用大挖掘机挖出恰好N单位体积的砂土。由于小挖掘机比较笨拙，它每次挖的砂土体积是固定的。也就是说，设每次挖x单位体积砂土，那么N需要被x整除。在挖出若干堆体积为x的砂土后，Matt需要计算x的“难挖指数”。体积x的“难挖指数”定义如下：对于某个不超过x的体积y，如果x与y的最大公约数为1，那我们认为体积y是“难挖的”，x的“难挖指数”就要加上y。

由于Matt之后还需要用小挖掘机处理被大挖掘机挖出的砂土，他希望知道所有可能的x的难挖指数的和，这样他好估算今天要干多久，不然作为布鲁谢特的高才生，他出门要被笑话的。

输入 第一行一个整数T，表示数据组数。

接下来T行每行一个整数表示N。

输出 对于每个数据输出一行一个整数表示难挖指数的和。

样例输入
3
2
3
4
样例输出
2
4
6
提示
对于30%的数据有T<=20，N<=10^4。
对于60%的数据有T<=100，N<=10^7。 对于100%的数据 有1<=T<=1000，1<=N<=10^9。

重要结论：ans=phi(n)*n/2;
phi()为欧拉函数
若a与n互质，则（n-a）也与n互质，而这样总共有phi(n)/2对
时间复杂度O(sqrt(n))

var
cnt,s,a:array[0..33] of int64;
ans,t:qword;
q,n,i,m,tt:longint;
procedure add(k:longint);
var
sum,ol:int64;
i:longint;
begin
  sum:=1;
  for i:=1 to k do
    sum:=sum*s[i];
  if sum=1 then begin ans:=ans+1; exit; end;
  ol:=sum;
  for i:=1 to k do
    if s[i]<>1 then ol:=ol*(a[i]-1) div a[i];
  ans:=ans+ol*sum div 2;
end;

procedure dfs(k:longint);
var
i:longint;
begin
  s[k]:=1;
  if k>m then begin add(k-1); exit; end;
  dfs(k+1);
  for i:=1 to cnt[k] do
  begin
    s[k]:=s[k]*a[k];
    dfs(k+1);
  end;
  s[k]:=1;
end;

begin
  readln(tt);
  for q:=1 to tt do
  begin
    readln(n);
    ans:=0;
    t:=n;
    m:=0;
    for i:=2 to trunc(sqrt(n)) do
    begin
      if t mod i=0 then
      begin
        m:=m+1;
        a[m]:=i;
        cnt[m]:=0;
        while t mod i=0 do
        begin
          cnt[m]:=cnt[m]+1;
          t:=t div i;
        end;
      if t=1 then break;
      end;
    end;
    if t<>1 then begin m:=m+1; a[m]:=t; cnt[m]:=1; end;
    if m=0 then begin m:=1; a[m]:=n; cnt[m]:=1; end;
    dfs(1);
    writeln(ans);
  end;
end.