首页 > 代码库 > HDU 5726 - GCD
HDU 5726 - GCD
RMQ应用题。
引用别人的解题思路:
1 #include<cstdio> 2 #include<cmath> 3 #include<algorithm> 4 #include<map> 5 #define MAXN 100000+5 6 using namespace std; 7 map<int,long long> m; 8 int n,a[MAXN],q; 9 int dp[MAXN][17]; 10 int gcd(int a,int b){return b==0?a:gcd(b,a%b);} 11 void RMQ_init() 12 { 13 for(int i=1;i<=n;i++) dp[i][0]=a[i]; 14 double j_max=(log(n)/log(2.0)); 15 for(int j=1;j<=j_max;j++) 16 { 17 for(int i=1;i<=n;i++) 18 if(i+(1 << j)-1 <= n) dp[i][j]=gcd(dp[i][j-1],dp[i+(1<<(j-1))][j-1]); 19 } 20 } 21 int query(int l,int r) 22 { 23 int k=log(r-l+1)/log(2.0); 24 return gcd(dp[l][k],dp[r-(1<<k)+1][k]); 25 } 26 void pretreat()//这个预处理是需要好好理解的地方 27 { 28 m.clear(); 29 for(int l=1;l<=n;l++) 30 { 31 int now=l; 32 while(now <= n) 33 { 34 int d=query(l,now); 35 int st=now,ed=n,mid; 36 while(ed-st>1) 37 { 38 mid=(st+ed)>>1; 39 if(query(l,mid) == d) st=mid; 40 else ed=mid; 41 } 42 m[d]+=st-now+1; 43 now=st+1; 44 } 45 } 46 } 47 int main() 48 { 49 int t; 50 scanf("%d",&t); 51 for(int kase=1;kase<=t;kase++) 52 { 53 scanf("%d",&n); 54 for(int i=1;i<=n;i++) scanf("%d",&a[i]); 55 RMQ_init(); 56 pretreat(); 57 scanf("%d",&q); 58 int l,r; 59 printf("Case #%d:\n",kase); 60 for(int i=1;i<=q;i++) 61 { 62 scanf("%d%d",&l,&r); 63 int ans=query(l,r); 64 printf("%d %I64d\n",ans,m[ans]); 65 } 66 } 67 }
HDU 5726 - GCD
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。