• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>

    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>

            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            

    
            

            




            







												


        

            
	
            
		
			

            

            
	
            
		
			
            Drolca
            
			Apologize To Drolca
		            		
	
            

            

            

            

            隨筆 - 28, 文章 - 1, 評論 - 6, 引用 - 0
            

            


            		
	
            

            

            
	
            
		
			
            
				數(shù)據(jù)加載中……
			
            
		            		
		
			
            
				
					
	
            
		
            
			pku 2478 Farey Sequence 歐拉函數(shù)
		
            
		
            
		
            #include <iostream>
            using namespace std;

            const int size=1000005;
            bool ok[size];
            int ph[size];
            int prime[size];
            __int64 sum[size];
            int pn;
            int n;
            void getph()
            {
                            int i,j;
                pn            =0;
                            for(i=2;i<size;i++)
                            {
                                if(!ok[i])
                                {
                        prime[pn            ++]=i;
                        ph[i]            =i-1;
                    }
                                for(j=0;(j<pn)&&(i*prime[j]<size);j++)
                                {
                        ok[i            *prime[j]]=1;
                                    if((i%prime[j])==0)
                                    {
                            ph[i            *prime[j]]=ph[i]*prime[j];
                                        break;
                        }
                                    else
                            ph[i            *prime[j]]=ph[i]*(prime[j]-1);

                    }
                }
            }
            int main()
            {
                ph[            1]=1;
                getph();
                            for(int i=2;i<size;i++)
                    sum[i]            =sum[i-1]+ph[i];
                            while(scanf("%d",&n)!=EOF&&(n!=0))
                    printf(            "%I64d\n",sum[n]);
                            return 0;
            }
            

		
            
		
            
			posted on 2009-09-18 13:20 Drolca 閱讀(198) 評論(0)  編輯 收藏 引用  
		
            
	
            
	
	


	


            
	    
    




            






            

				
			
            
		            		
	
            
	
            
		
			

		
	
            

            

    
    







            
	   
	



亚洲午夜久久久久久噜噜噜|
狠狠色婷婷久久一区二区三区|
精品无码久久久久久久动漫|
久久久国产精品|
一本一本久久aa综合精品|
久久综合九色综合网站|
国产精品亚洲综合专区片高清久久久|
久久久久无码专区亚洲av|
久久无码专区国产精品发布|
精品久久久久久国产|
亚洲午夜福利精品久久|
2021久久国自产拍精品|
亚洲а∨天堂久久精品9966|
久久99亚洲网美利坚合众国|
久久综合伊人77777|
国产精品久久国产精品99盘|
一本色道久久综合|
亚洲国产精品人久久|
亚洲欧美成人综合久久久|
狠狠人妻久久久久久综合|
91精品国产高清久久久久久io|
亚洲午夜无码AV毛片久久|
婷婷综合久久狠狠色99h|
人妻无码αv中文字幕久久|
亚洲精品综合久久|
久久精品国产精品亚洲人人|
99久久综合国产精品二区|
久久精品中文騷妇女内射|
色诱久久久久综合网ywww|
久久久久久免费视频|
欧美性大战久久久久久|
久久国产精品无码网站|
久久精品二区|
久久久久99精品成人片牛牛影视|
色成年激情久久综合|
久久综合九色综合97_久久久|
99久久无色码中文字幕|
国产精品久久自在自线观看|
狠狠色丁香婷婷综合久久来|
青青热久久综合网伊人|
精品久久久久久无码免费|