• <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)  編輯 收藏 引用  
		
            
	
            
	
	


	


            
	    
    




            






            

				
			
            
		            		
	
            
	
            
		
			

		
	
            

            

    
    







            
	   
	



久久AV高清无码|
久久人人爽人人澡人人高潮AV|
亚洲人成伊人成综合网久久久|
国产成人精品综合久久久|
国产精品视频久久久|
久久精品夜色噜噜亚洲A∨|
国内精品伊人久久久久妇|
精品熟女少妇a∨免费久久|
久久久久人妻一区精品|
亚洲国产精品久久电影欧美|
久久99国产一区二区三区|
性色欲网站人妻丰满中文久久不卡|
久久久久一区二区三区|
午夜精品久久久久久中宇|
久久人人爽人人爽AV片|
中文字幕亚洲综合久久|
色综合久久久久综合体桃花网
|
久久99国产精品久久99|
久久久久亚洲精品日久生情|
精品欧美一区二区三区久久久|
欧洲精品久久久av无码电影|
欧美久久一区二区三区|
91亚洲国产成人久久精品|
奇米综合四色77777久久|
久久久久久免费视频|
人妻无码久久精品|
国产综合成人久久大片91|
国产麻豆精品久久一二三|
日韩人妻无码一区二区三区久久|
久久99这里只有精品国产|
久久精品人妻一区二区三区|
国产成人无码精品久久久免费|
老司机国内精品久久久久|
国产亚洲综合久久系列|
久久久噜噜噜www成人网|
久久亚洲精品人成综合网|
人妻少妇久久中文字幕一区二区|
久久综合狠狠综合久久综合88
|
要久久爱在线免费观看|
亚洲精品高清一二区久久|
久久久久久国产精品无码下载|