Yuan
|
首頁
|
發(fā)新隨筆
|
發(fā)新文章
|
聯(lián)系
|
聚合
|
管理
noi 2010 能量采集 八中2005 ★★★ 公約數(shù)確定的對數(shù)
/**/
/*
題意:原題可轉(zhuǎn)化為求在n*m范圍內(nèi)
n m
∑∑ 2*(gcd(x,y)-1)+1
x=1 y=1
感覺挺像visible trees的用容斥做,但不會容斥做
看了
http://hi.baidu.com/570193465/blog/item/d4219303d7b43e1c738b6547.html
O(nsqrt(n))
對于上式,重點是求出 t=gcd(x,y) 時的(x,y)對數(shù)
可以枚舉gcd
記錄cnt[i] = (n/i)*(m/i) 即公約數(shù)是i的倍數(shù)(k*i)的對數(shù)
然后再調(diào)整cnt[i],使其表示公約數(shù)是i的對數(shù)
cnt[i]-=cnt[k*i]即可!
*/
#include
<
cstdio
>
#include
<
cstring
>
inline
int
min(
int
a,
int
b)
{
return
a
<
b
?
a:b;}
const
int
MAXN
=
100010
;
long
long
cnt[MAXN];
int
main()
{
int
n,m;
while
(
~
scanf(
"
%d%d
"
,
&
n,
&
m))
{
int
t
=
min(n,m);
for
(
int
i
=
2
;i
<=
t;i
++
)
//
gcd = i
cnt[i]
=
(
long
long
)(n
/
i)
*
(m
/
i);
//
cnt[i] : the number of whose gcd is k*i
for
(
int
i
=
t;i
>=
1
;i
--
)
for
(
int
k
=
2
;k
*
i
<=
t;k
++
)
cnt[i]
-=
cnt[k
*
i];
//
to get the real number of whose gc is i
long
long
ans
=
0
;
for
(
int
i
=
1
;i
<=
t;i
++
)
ans
+=
2
*
(i
-
1
)
*
cnt[i];
printf(
"
%I64d\n
"
,ans
+
(
long
long
)n
*
m);
}
return
0
;
}
發(fā)表于 2010-09-07 23:59
_Yuan
閱讀(566)
評論(2)
編輯
收藏
引用
所屬分類:
OJ解題報告
評論
#
re: noi 2010 能量采集 八中2005 ★★★ 公約數(shù)確定的對數(shù)
贊一個~代碼好短
LitIce
評論于 2010-11-15 14:57
回復
更多評論
#
re: noi 2010 能量采集 八中2005 ★★★ 公約數(shù)確定的對數(shù)
@LitIce
那個是看別人寫的 #_#
_Yuan
評論于 2010-11-15 15:49
回復
更多評論
刷新評論列表
只有注冊用戶
登錄
后才能發(fā)表評論。
【推薦】100%開源!大型工業(yè)跨平臺軟件C++源碼提供,建模,組態(tài)!
相關(guān)文章:
SRM 239 HiddenTriangles ★★★★
CodeForces 59E 以邊為狀態(tài)bfs ★★★★
TCO'10 Wildcard Round 500pt CalculationCards
zoj 3462 bitset
SRM 496 PalindromfulString 容斥寫法 ★★★★
CodeForces 57D
CodeForces 55D 數(shù)位統(tǒng)計 記憶化搜索 跟pre有關(guān) ★★★★
CodeForces 55E Very simple problem
zoj 3455 統(tǒng)計出現(xiàn)次數(shù) 判斷相等 用l[i]記錄字母出現(xiàn)i次的個數(shù) ★★★★
zoj 3354 映射 環(huán) 計數(shù) ★★★
網(wǎng)站導航:
博客園
IT新聞
BlogJava
博問
Chat2DB
管理
常用鏈接
我的隨筆
我的評論
我參與的隨筆
隨筆分類
Dp(27)
(rss)
OJ解題報告(153)
(rss)
OThers(17)
(rss)
TopCoder
(rss)
計算幾何(2)
(rss)
枚舉(4)
(rss)
數(shù)據(jù)結(jié)構(gòu)(6)
(rss)
數(shù)論(5)
(rss)
搜索(2)
(rss)
貪心(4)
(rss)
圖論(10)
(rss)
學習筆記(6)
(rss)
學習總結(jié)(19)
(rss)
組合數(shù)學(3)
(rss)
Links
Lord Li
Lord zeus
搜索
最新評論
1.?re: 雙向BFS[未登錄]
博主,只用一個隊列不就可以解決你第一個問題了嗎
--jason
2.?re:nvgagkguaioguaiiananfajfofajiosfgoasoajgia[未登錄]
cscdcuis
--1
3.?re: zoj 3436 逆推 搜
評論內(nèi)容較長,點擊標題查看
--ZH
4.?re: zoj 2318 計算幾何 spfa判負環(huán)
寫得好!
--ipqhjjybj
5.?re: Poj 1066
@楊書鑒
你寫的排序好像不對啊。。。
--小猊
Powered by:
博客園
模板提供:
滬江博客
Copyright ©2025 _Yuan
亚洲国产精品综合久久一线
|
国产成人精品白浆久久69
|
久久综合伊人77777麻豆
|
看久久久久久a级毛片
|
国产精品成人99久久久久
|
久久er99热精品一区二区
|
日本精品久久久久中文字幕8
|
亚洲?V乱码久久精品蜜桃
|
亚洲AV无码久久精品成人
|
精品亚洲综合久久中文字幕
|
少妇内射兰兰久久
|
一本久久精品一区二区
|
久久精品夜色噜噜亚洲A∨
|
色综合久久天天综合
|
久久电影网2021
|
国产成人无码精品久久久久免费
|
久久精品国产99久久久古代
|
久久亚洲精品国产亚洲老地址
|
激情久久久久久久久久
|
国产高潮国产高潮久久久91
|
久久香蕉国产线看观看99
|
国产精品伦理久久久久久
|
久久精品中文字幕第23页
|
久久精品综合网
|
国产精品久久午夜夜伦鲁鲁
|
A狠狠久久蜜臀婷色中文网
|
亚洲国产成人久久综合碰碰动漫3d
|
久久亚洲精品无码VA大香大香
|
精品久久8x国产免费观看
|
国产精品久久久久一区二区三区
|
久久人做人爽一区二区三区
|
性欧美大战久久久久久久久
|
AV狠狠色丁香婷婷综合久久
|
国产2021久久精品
|
久久亚洲私人国产精品vA
|
观看 国产综合久久久久鬼色 欧美 亚洲 一区二区
|
国产国产成人精品久久
|
精品久久久久久久久免费影院
|
97精品久久天干天天天按摩
|
深夜久久AAAAA级毛片免费看
|
久久精品成人国产午夜
|