• <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
            

            


            		
	
            

            

            
	
            
		
			
            
				數據加載中……
			
            
		            		
		
			
            
				
					
	
            
		
            
			pku 2478 Farey Sequence 歐拉函數
		
            
		
            
		
            #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 閱讀(206) 評論(0)  編輯 收藏 引用  
		
            
	
            
	
	


	


            
	    
    




            






            

				
			
            
		            		
	
            
	
            
		
			

		
	
            

            

    
    







            
	   
	



99久久er这里只有精品18|
久久精品夜色噜噜亚洲A∨|
久久午夜无码鲁丝片午夜精品|
国产色综合久久无码有码|
久久精品国产2020|
99久久国产综合精品麻豆|
99久久国产综合精品网成人影院|
久久久久久久精品妇女99
|
国产精品久久久久蜜芽|
狠狠综合久久综合88亚洲|
成人精品一区二区久久|
精品国产乱码久久久久久呢|
国产精品一久久香蕉国产线看观看
|
久久97久久97精品免视看|
亚洲色欲久久久久综合网|
久久婷婷五月综合色高清
|
噜噜噜色噜噜噜久久|
久久久91精品国产一区二区三区
|
久久这里只精品99re66|
夜夜亚洲天天久久|
国产99久久久久久免费看|
91精品国产91久久综合|
亚洲精品无码久久千人斩|
亚洲第一极品精品无码久久|
久久99国产精品久久99小说|
亚洲欧美成人久久综合中文网
|
久久这里只有精品久久|
亚洲αv久久久噜噜噜噜噜|
亚洲国产精品综合久久网络|
91麻豆国产精品91久久久|
人妻无码αv中文字幕久久琪琪布
人妻无码精品久久亚瑟影视
|
久久亚洲国产精品成人AV秋霞|
久久午夜福利无码1000合集|
亚洲午夜久久久久久久久电影网|
亚洲精品美女久久777777|
潮喷大喷水系列无码久久精品
|
亚洲国产天堂久久久久久|
一本色道久久88综合日韩精品
|
久久久一本精品99久久精品88|
久久精品国产亚洲AV蜜臀色欲|
大美女久久久久久j久久|