Lưu ý: Các code mẫu dưới đây chỉ mang tính tham khảo và có thể không AC được bài tập này

Code mẫu của RR

{$R+,Q+}
var
  n,u,sum:longint;
  kq:int64;
function gcd(a,b:longint):longint;
begin
  if (a=0) or (b=0) then exit(a+b)
  else if (a<b) then exit(gcd(b,a))
  else exit(gcd(b,a mod b));
end;
procedure count;
var
  i:longint;
begin
  for i:=1 to sum div 2 do
    if gcd(i,sum-i)=1 then inc(kq);
end;
begin
  read(n);
  kq:=0;
  for u:=trunc(sqrt(n)) downto 1 do
  if n mod u=0 then
    begin
      sum:=n div u+1;
      count;
      if u*u=n then continue;
      sum:=u+1;
      count;
    end;
  writeln(kq);
end.

Code mẫu của hieult

#include <stdio.h>
//#include <conio.h>
#include <math.h>
long n[200];
long Ole(long N)
{
long k=1,KQ=1;
if(N==2)
  return 2;
else
  {                    
  for(long i=1;i<=180;i++)
      {
      if(N==1)
        break;       
      else if(n[i]>sqrt(N))
          {
          KQ*=(N-1);              
          break;
          }     
      else if(N%n[i]==0)
          {
          long t=0;                
          do
            {
            N=N/n[i];
            t++;
            }while(N%n[i]==0);
          KQ=KQ*(n[i]-1)*pow(n[i],t-1);
          }                      
      }
  return KQ;
  }
}               
main()
 {
 long t=2,T,KQ=0;
 n[1]=2;
 n[2]=3;
 for(long i=4;i<1100;i++)
   {
   int k=0;        
   for(long j=1;j<=t;j++)
     {
     if(i%n[j]==0)
       {
       k=1;
       break;
       }
     else if(n[j]>sqrt(i))     
       break;
     }
   if(k==0)
     {
     t++;
     n[t]=i;
     }        
   }                 
 scanf("%ld",&T);
 for(long i=1;i<=sqrt(T);i++)
   {
   if(i*i==T)
     KQ=KQ+Ole(i+1)/2;
   else if(T%i==0)
     KQ=KQ+Ole(i+1)/2+Ole(T/i+1)/2;  
   }
 printf("%ld",KQ);   
 //getch();
 }

Code mẫu của ll931110

Program SQUARES;
  Const
    input  = '';
    output = '';
    maxn = 1000000;
  Var
    a,F: array[1..maxn] of longint;
    t,n,i: longint;
    fi,fo: text;

Procedure openfile;
  Begin
    Assign(fi, input);
      Reset(fi);

    Assign(fo, output);
      Rewrite(fo);
  End;

Procedure solve;
  Var
    i,j,k,tmp,res: longint;
  Begin
    Fillchar(a, sizeof(a), 0);
    For i:= 2 to trunc(sqrt(maxn)) do
      if a[i] = 0 then
      Begin
        k:= 2 * i;
        While k <= maxn do
          Begin
            a[k]:= i;
            k:= k + i;
          End;
      End;

    F[1]:= 1;
    For i:= 2 to maxn do
      if a[i] = 0 then F[i]:= i - 1 else
        Begin
          k:= 0;
          tmp:= i;

          While tmp mod a[i] = 0 do
            Begin
               inc(k);
               tmp:= tmp div a[i];
            End;

          If tmp = 1 then
            Begin
              F[i]:= a[i] - 1;
              For j:= 1 to k - 1 do F[i]:= F[i] * a[i];
            End
          else F[i]:= F[tmp] * F[i div tmp];
        End;

      Readln(fi, n);
      res:= 0;
      For i:= 1 to n - 1 do
        if n mod i = 0 then res:= res + F[n div i + 1] div 2;

      Writeln(fo, res + 1);
  End;

Procedure closefile;
  Begin
    Close(fo);
    Close(fi);
  End;

Begin
  openfile;
  solve;
  closefile;
End.