BestCoder Round #78 (div.1) CA Loves GCD

本文介绍了一种通过因子分解来求解最大公约数的算法,并提供了完整的C++实现代码。该算法首先计算出每个数的所有约数,然后利用这些约数递归地找出所有可能的最大公约数组合。

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

因子分解,然后把约数为1到1000的方案都求出来,比如约数为2的方案为2^(n个数中能被2整除的个数)-1,然后从1000到1逆推,推的过程要更更新小约数的值比如推完am[8]时,am[2],am[4],am[1]都要减去am[8](因为m个数的公约数为8时那公约数一定也为2和4,)逆推可以保证最大公约数。
这么水的题都要想很长时间我真是太渣了。。。

#include<cstdio>
#include<cstring>
#define mod 100000007
int am[1005],as[1005];
long long ans;
int main()
{
   int n,t;
   scanf("%d",&t);
   while(t--)
   {
        ans = 0;
        memset(as,0,sizeof(as));
        scanf("%d",&n);
         while(n--)
          {
               int a;
             scanf("%d",&a);
          for(int j=1;j<=a;j++)
             if(!(a%j))as[j]++;   
          }
        for(int i=1;i<=1000;i++)
         {
             long long w = 1;
          for(int j=0;j<as[i];j++)
            {
               w*=2;
               w%=mod;
            }
            w-=1;
            w=(w%mod+mod)%mod;
          am[i] = w;
      }
      for(int i=1000;i>=1;i--)
      {
          ans=ans+((long long)am[i]*i)%mod;
          ans%=mod;
          for(int j=1;j<=i;j++)
           if(!(i%j))
           {
             am[j]-=am[i];
             am[j]=(am[j]%mod+mod)%mod;
           }
      }
      printf("%I64d\n",ans);
   }
   return 0;    
} 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值