hdu-5768 Lucky7 容斥

http://acm.hdu.edu.cn/showproblem.php?pid=5768

题意是要求给定范围内,满足模7为零,模pi不为ai的数的个数(pi,ai已知条件)

可以先求出范围内所有模7为零的元素,在求出所有模7为0并且模pi为ai的数的个数,最后减一下

题目相当于转化为解多个同余方程组加容斥的问题

数的个数比较少,可以直接dfs暴力容斥

#include<bits/stdc++.h>
using namespace std;
long long ans=0;
int n;
long long X,Y;
long long f(long long mo,long long yu)
{
    return (Y-yu+mo)/mo-(X-yu-1+mo)/mo;
}
int p[30],a[30];
int Extended_Euclid(long long a,long long b,long long &x,long long &y)
{
    long long d;
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    d=Extended_Euclid(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
long long quick_add(long long k,long long x,long long mm)
{
    if(x<0)
        k=-k,x=-x;
    long long ans=0;
    while(x)
    {
        if(x&1)
            ans=(ans+k)%mm;
        k=(k+k)%mm;
        x>>=1;
    }
    return ans;
}
void dfs(int num,int step,long long mo,long long yu)
{
    if(step==n)
    {
     //   cout<<num<<" "<<mo<<" "<<yu<<endl;
        if(num)
        ans+=((num&1)?-1:1)*(f(mo,yu));
        return ;
    }
    long long xx,yy,ret,nn;
    Extended_Euclid(mo,p[step],xx,yy);
  //  cout<<xx<<" "<<yy<<endl;
    nn=mo*p[step];
    ret=quick_add(xx,quick_add(a[step]-yu,mo,nn),nn)+yu;

    dfs(num+1,step+1,nn,(nn+ret%nn)%nn);
    dfs(num,step+1,mo,yu);
}
int main()
{
    int tca=1,T,i;
    long long xx,yy;
    Extended_Euclid(7ll,3ll,xx,yy);
   // cout<<xx<<" "<<yy<<endl;
    cin>>T;
    while(T--)
    {
        cin>>n>>X>>Y;
        for(i=0;i<n;i++)
        {
            scanf("%d %d",&p[i],&a[i]);
        }
        ans=f(7,0);
        dfs(0,0,7,0);
        printf("Case #%d: %lld\n",tca++,ans);
    }
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值