hdoj-1796 容斥水题

原创于 2016-03-11 22:11:50 发布 · 288 阅读

0 ·

CC 4.0 BY-SA版权

数学专栏收录该内容

17 篇文章

订阅专栏

本文介绍了一个使用C++实现的算法示例，通过深度优先搜索（DFS）结合最小公倍数（LCM）计算的方法来解决特定问题。代码展示了如何通过递归调用进行状态遍历，并在过程中计算不同整数的最小公倍数。

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

第一发略坑

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
__int64 n,m,ans,sum,a[20],b;
__int64 lcm(__int64 a,__int64 b)
{
    int a1 = a,b1 = b;
    while(b){
        long long c = b;
        b = a%b;
        a = c;
    }
    return a1/a*b1/a*a;
}
void dfs(long long pre,long long am,long all,long pre2)
{
    if(am==all)
    {
        long long g = 1,re=1;
        for(int i=0;i<m;i++)
        {
           if(((1<<i)&pre))
           g = lcm(g,a[i]);
           if(g>=n)return;  
        }
        if(all%2){ans+=(n-1)/g;}
        else {ans-=(n-1)/g;}
    }
    else
    {
        for(int i=pre2;i<m;i++)
          if(!((1<<i)&pre)){
             dfs(pre|(1<<i),am+1,all,i+1);
          }
    }
}
int main()
{
    while(scanf("%I64d%I64d",&n,&m)!=EOF)
    {
        int i1 = 0;
       for(int i=0;i<m;i++)
     {
          scanf("%I64d",&b);
          if(b>0&&b<n){
           a[i1] = b;
           i1++;
         }
     }
       m = i1;
         sort(a,a+m);
       ans = 0;
       for(int i=1;i<=m;i++)
            dfs(0,0,i,0);  
        printf("%I64d\n",ans);
    }
    return 0;
}