• <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 閱讀(198) 評論(0)  編輯 收藏 引用  
		
            
	
            
	
	


	


            
	    
    




            






            

				
			
            
		            		
	
            
	
            
		
			

		
	
            

            

    
    







            
	   
	



欧美日韩精品久久久久|
久久精品国产第一区二区三区|
一级a性色生活片久久无少妇一级婬片免费放
|
国产成人精品久久亚洲高清不卡 |
精品国产青草久久久久福利|
久久人人爽人人爽人人爽|
成人久久精品一区二区三区|
色天使久久综合网天天|
99久久777色|
亚洲国产精品一区二区三区久久|
亚洲色婷婷综合久久|
精品久久久无码21p发布|
亚洲精品无码久久久久sm|
国产精品久久久福利|
伊人伊成久久人综合网777|
午夜精品久久久久久99热|
国内精品久久久久久久coent|
97精品依人久久久大香线蕉97|
久久99精品九九九久久婷婷|
69久久精品无码一区二区|
国产69精品久久久久APP下载|
亚洲午夜久久久精品影院|
97精品伊人久久大香线蕉app|
国产精品久久久久久五月尺|
精品无码人妻久久久久久|
久久精品九九亚洲精品天堂|
精品久久久久香蕉网|
久久综合久久综合久久综合|
国产成人精品久久免费动漫|
久久综合亚洲色一区二区三区|
久久久受www免费人成|
久久中文字幕人妻熟av女|
日韩欧美亚洲综合久久影院Ds|
久久se精品一区二区影院|
99精品久久久久中文字幕|
久久国产精品无码HDAV
|
久久国产精品成人影院|
亚洲AV无码1区2区久久|
伊人久久综合无码成人网|
中文字幕久久久久人妻|
久久精品中文騷妇女内射|