1213 解的个数[一中数论随堂练]

1213 解的个数

时间限制: 1 s

空间限制: 128000 KB

题目等级 : 黄金 Gold

题解

查看运行结果

题目描述 Description

已知整数x,y满足如下面的条件:

ax+by+c = 0

p<=x<=q

r<=y<=s

求满足这些条件的x,y的个数。

输入描述 Input Description

第一行有一个整数nn<=10),表示有n个任务。n<=10

以下有n行,每行有7个整数,分别为:a,b,c,p,q,r,s。均不超过108。

输出描述 Output Description

n行,第i行是第i个任务的解的个数。

样例输入 Sample Input

2

2 3 -7 0 10 0 10

1 1 1 -10 10 -9 9

样例输出 Sample Output

1

19

数据范围及提示 Data Size & Hint

分类标签 Tags 点此展开

欧几里德定理 数论

两种AC代码:

1、扩展欧几里得算法,需要考虑特殊数据

#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
typedef long long ll;
double z[10];
void exgcd(ll a,ll b,ll &d,ll &x,ll &y){
    if(!b){d=a;x=1;y=0;return;}
    exgcd(b,a%b,d,y,x);y-=x*(a/b);
}
int main(){
    ll a,b,c,p,q,r,s,x,y,d,k;
    int n;cin>>n;
    while(n--){
        cin>>a>>b>>c>>p>>q>>r>>s;c=-c;
        if(p>q||r>s||(!a&&!b&&c)){cout<<0<<endl;continue;}
        if(!a||!b){
            if(!a) a=q-p+1;
            else if(c%a==0&&(c/a)<=q&&(c/a)>=p) a=1;
            else a=0;
            if(!b) b=s-r+1;
            else if(c%b==0&&(c/b)<=s&&(c/b)>=r) b=1;
            else b=0;
            cout<<a*b<<endl;
            continue;
        }
        exgcd(a,b,d,x,y);
        if(c%d!=0){cout<<0<<endl;continue;}
        k=c/d;x*=k;y*=k;
        a=a/d;b=b/d;
        z[0]=(q-x)*1.0/b;
        z[1]=(p-x)*1.0/b;
        z[2]=(y-s)*1.0/a;
        z[3]=(y-r)*1.0/a;
        sort(z,z+4);
        a=floor(z[2]);b=ceil(z[1]);
        cout<<a-b+1<<endl;
    }
    return 0;
}

2、直接枚举

#include<cstdio>
#include<iostream>
#define ll long long
using namespace std;
ll n,a,b,c,p,q,r,l,ans;
int main(){
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a>>b>>c>>p>>q>>l>>r;
        ans=0;
        if(b==0){
            for(int x=p;x<=q;x++){
                if(x*a==-c) ans++;
            }
            cout<<ans*(r-l+1)<<endl;
        }
        else{
            for(int x=p;x<=q;x++){
                if((-c-a*x)%b==0&&(-c-a*x)*1.0/b>=l&&(-c-a*x)*1.0/b<=r) ans++;
            }
            cout<<ans<<endl;
        }
    }
    return 0;
}
时间: 2024-10-04 14:14:12

1213 解的个数[一中数论随堂练]的相关文章

1200 同余方程[一中数论随堂练]

1200 同余方程 2012年NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题解 查看运行结果 题目描述 Description 求关于 x 同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入描述 Input Description 输入只有一行,包含两个正整数 a, b,用 一个 空格隔开. 输出描述 Output Description 输出只有一行包含一个正整数x0,即最小正整数解,输入数据保证一定有解. 样例输入

3040 中国余数定理 1[一中数论随堂练]

3040 中国余数定理 1 时间限制: 1 s 空间限制: 32000 KB 题目等级 : 青铜 Bronze 题解 题目描述 Description 摘自算法导论...... 找出第k个被3,5,7除的时候,余数为2,3,2的数: 输入描述 Input Description 一个数k. 输出描述 Output Description 求出第k个符合条件的数. 样例输入 Sample Input 1 样例输出 Sample Output 23 数据范围及提示 Data Size & Hint

1702 素数判定 2[[一中数论随堂练]

1702 素数判定 2 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题解 题目描述 Description 一个数,他是素数么? 设他为P满足(P<=263-1) 输入描述 Input Description P 输出描述 Output Description Yes|No 样例输入 Sample Input 2 样例输出 Sample Output Yes 数据范围及提示 Data Size & Hint 算法导论——数论那一节注意Carmich

4939 欧拉函数[一中数论随堂练]

4939 欧拉函数 时间限制: 1 s 空间限制: 1000 KB 题目等级 : 钻石 Diamond 题解 题目描述 Description 输入一个数n,输出小于n且与n互素的整数个数 输入描述 Input Description 包含多组数据,n=0时结束 测试数据组数不会很多,不必先打表后输出 输出描述 Output Description 一组数据一行 样例输入 Sample Input 364684 346 5432 11 24 0 2333333 233333333 0 23333

3990 中国余数定理 2[一中数论随堂练]

3990 中国余数定理 2 时间限制: 1 s 空间限制: 1000 KB 题目等级 : 白银 Silver 题解 查看运行结果 题目描述 Description Skytree神犇最近在研究中国博大精深的数学. 这时,Sci蒟蒻前来拜访,于是Skytree给Sci蒟蒻出了一道数学题: 给定n个质数,以及k模这些质数的余数.问:在闭区间[a,b]中,有多少个k?最小的k是多少? Sci蒟蒻数学能力差了Skytree三条街,所以他只好寻求计算机的帮助.他发邮件给同为oier的你,你能帮他解决这个问

扩展gcd codevs 1213 解的个数

codevs 1213 解的个数 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 黄金 Gold 题目描述 Description 已知整数x,y满足如下面的条件: ax+by+c = 0 p<=x<=q r<=y<=s 求满足这些条件的x,y的个数. 输入描述 Input Description 第一行有一个整数n(n<=10),表示有n个任务.n<=10 以下有n行,每行有7个整数,分别为:a,b,c,p,q,r,s.均不超过108. 输出描述 O

codevs 1213 解的个数(我去年打了个表 - -)

#include<iostream> #include<cstdio> #include<cstring> using namespace std; int T,x,y,ans,g,l1,r1,l2,r2; int init() { int x=0;char s=getchar();bool f=0; while(s<'0'||s>'9'){if(s=='-')f=1;s=getchar();} while(s>='0'&&s<=

扩展欧几里德--解的个数

解的个数 已知整数x,y满足如下面的条件: ax+by+c = 0 p<=x<=q r<=y<=s 求满足这些条件的x,y的个数. 输入描述 Input Description 第一行有一个整数n(n<=10),表示有n个任务.n<=10 以下有n行,每行有7个整数,分别为:a,b,c,p,q,r,s.均不超过108. 输出描述 共n行,第i行是第i个任务的解的个数. 样例输入 2 2 3 -7 0 10 0 10 1 1 1 -10 10 -9 9 样例输出 1 19

HDU1299&amp;&amp;SPOJ - KPEQU(找因子来找解的个数)

题目链接:HDU 1299 传送门 题意: 求方程1/x+1/y=1/n的解的个数 1/3+1/2 与1/2+1/3看作是一组解. 分析: 1/x+1/y = 1/n 设y = n + k; ==>1/x + 1/(n+k)=1/n; ==>x = n^2/k + n; 因为x为整数,k就是n^2的约数. 代码如下: #include <iostream> #include <cstring> #include <cstdio> #include <a