2019年华南理工大学程序设计竞赛 H Parco_Love_GCD(区间gcd_裸)

最新推荐文章于 2021-08-13 13:21:41 发布

Jayson_Jun

最新推荐文章于 2021-08-13 13:21:41 发布

阅读量193

点赞数

CC 4.0 BY-SA版权

分类专栏： acm 数论

本文链接：https://blog.youkuaiyun.com/weixin_38772011/article/details/89366514

acm 同时被 2 个专栏收录

35 篇文章

订阅专栏

数论

3 篇文章

订阅专栏

本文主要解析2019年华南理工大学程序设计竞赛中H题——Parco_Love_GCD问题，涉及区间最大公约数（GCD）的裸题解法。通过链接可以访问到官方题目详情。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

url:https://ac.nowcoder.com/acm/contest/625/H

1.1

2.2 1.2

3.3 2.3 1.3

4.4 3.4 2.3 1.3

.......

	#include<bits/stdc++.h>
	using namespace std;
	typedef long long ll;
	#define maxn 500005
	const int mod = 1000000007;
	typedef pair<int,int>P;
	vector<P>vv[maxn];
	ll n,nn,a[maxn];
	ll l,r;
	ll res;
	
	ll gcd(ll a,ll b)
	{
	    return b?gcd(b,a%b):a;
	}
	void update(ll v)
	{
	    ll last=r;
	    ll tmp=0;
	    for(ll i=0;i<vv[r].size();i++)
	    {
	        ll g=vv[r][i].first,p=vv[r][i].second;
	        tmp+=1ll*(last-p+1)*g,last=p-1;
	//        cout<<tmp<<endl;
	    }
	//    cout<<endl;
	    res = (res%mod + 1ll*v*tmp%mod)%mod;
	}
	int main()
	{
	    scanf("%lld",&n);
	    nn=(int)(sqrt(1.0*n)+0.5);
	    for(ll i=1;i<=n;i++)
	    {
	        scanf("%lld",&a[i]);
	        vv[i].clear();
	    }
	    vv[0].clear();
	    for(ll i=1;i<=n;i++)
	    {
	        ll x=a[i],y=i;
	        for(ll j=0;j<vv[i-1].size();j++)
	        {
	            ll g=gcd(vv[i-1][j].first,a[i]);
	            if(g!=x)
	            {
	                vv[i].push_back(P(x,y));
	//                cout<<x<<" "<<y<<endl;
	                x=g;
	            }
	            y=vv[i-1][j].second;
	        }
	        vv[i].push_back(P(x,y));
	//        cout<<x<<" "<<y<<endl;
	//        cout<<endl;
	    }
	    res=0;
	    while(r<=n)
	    {
	        r++,update(1);
	    }
	    printf("%lld\n",res);
	    return 0;
	}