piezoelectric energy harvesting

[Noi2010] Energy harvesting Time Limit: ten Sec Memory Limit:552 MBSubmit: 2324 Solved: 1387[Submit] [Status] [Discuss] DescriptionThe building has a rectangular land, he planted an energy plant on the ground, which can collect the energy of sunlight. After the plants collected

machine is large and inconvenient to move, the building puts it on a corner with coordinates exactly (0, 0). The energy pooling machine has a certain energy loss during the collection process. If a plant is connected to an energy collection machine with a K plant on its line, the loss of energy is 2k + 1. For example,

Analysis: http://www.cnblogs.com/huhuuu/archive/2011/11/25/2263803.htmlNote: Two points have been harvested from this question1, the number of the hour on the first quadrant (x, y) to (0,0) line is gcd (x, y)2, the newly learned gcd (x, y) =k How many pairs of posture, known 0Make X=min (n,m), so f[i] represents the logarithm of gcd (x, y) =i,Then the complexity of O (XLOGX) can be obtained f[1] to F[n] (reverse loop)Ordinary Möbius (ie, inversion) is actually O (XLOGX), just need to sift throug

intMU[MAXN];Ten intn,m; One A intGETMU () { -mu[1]=1; - for(intI=2; i){ the if(!Flag[i]) { -mu[i]=-1; -prime[++cnt]=i; - } + for(intj=1; i*prime[j]){ -flag[i*prime[j]]=1; + if(i%prime[j]==0){ Amu[i*prime[j]]=0; at Break; - } -mu[i*prime[j]]=-Mu[i]; - } - } - } in -ll work (intx) { toll ret=0; + intn=n/x,m=m/x; - for(intI=1; i) theret+=1ll*mu[i]* (n/i) * (m/i); * returnret; $ }Panax Notoginseng - i

) for (int i=j;i - #defineD (i,j,n) for (int i=j;i>=n;--i) + using namespacestd; -typedefLong LongLL; +InlineintGetint () { A intR=1, v=0;CharCh=GetChar (); at for(;! IsDigit (CH); Ch=getchar ())if(ch=='-') r=-1; - for(; isdigit (ch); Ch=getchar ()) v=v*Ten+ch-'0'; - returnr*v; - } - Const intn=1e5+Ten, inf=~0u>>2; - /*****************template**********************/ in intphi[n],prime[n],tot,n,m; - BOOLCheck[n]; to voidGetphi (intN) { +memset (check,0,sizeofcheck); -phi[1]=1; t