【数论】 HDOJ 1573 X问题
来源:互联网 发布:时时彩预测软件 编辑:程序博客网 时间:2024/05/17 23:49
非互质的中国剩余定理。。。用那个叫合并上升法的东西。。。orzorzorz。。。
#include <iostream> #include <sstream> #include <algorithm> #include <vector> #include <queue> #include <stack> #include <map> #include <set> #include <bitset> #include <cstdio> #include <cstring> #include <cstdlib> #include <cmath> #include <climits> #define maxn 800005#define eps 1e-6 #define mod 10007 #define INF 99999999 #define lowbit(x) (x&(-x)) //#define lson o<<1, L, mid //#define rson o<<1 | 1, mid+1, R typedef long long LL;using namespace std;void extend_gcd(int a, int b, int &d, int &x, int &y){if(b == 0) { d = a, x = 1, y = 0; }else { extend_gcd(b, a%b, d, y, x), y -= x*(a/b); }}int a[maxn], b[maxn];int main(void){int _, a1, a2, b1, b2, tmp, i, n, m, ok, g, x, y;while(scanf("%d", &_)!=EOF) {while(_--) {scanf("%d%d", &m, &n);for(i = 0; i < n; i++) scanf("%d", &a[i]);for(i = 0; i < n; i++) scanf("%d", &b[i]);a1 = a[0], b1 = b[0], ok = 0;for(i = 1; i < n; i++) {a2 = a[i], b2 = b[i];extend_gcd(a1, a2, g, x, y);if((b2-b1)%g) { //无解的情况ok = 1;break;}tmp = a2/g;x = x*(b2-b1)/g; //x 的所有解x = (x%tmp+tmp)%tmp; //找到最小正解xb1 = a1*x+b1; //更新b1a1 = (a1*a2)/g; //更新a1b1 = (b1%a1+a1)%a1; //找到最小正解b1}//解出的答案是b1,循环是a1.if(ok || m<b1) printf("0\n");else printf("%d\n", (m-b1)/a1+1-(b1==0 ? 1 : 0));}}return 0;}
0 0
- 【数论】 HDOJ 1573 X问题
- HDU 1573 X问题 数论
- HDU 1573 X问题(数论)
- HDU-1573 X问题(数论:最小公倍数)
- HDOJ-----1573---X问题(扩展欧几里得)
- HDU 1573 X问题 数论-(中国剩余定理)
- HDU 1573 X问题 [一元线性同余方程组]【数论】
- hdoj-1573-X问题【中国剩余定理 & 除数不互质】
- HDOJ 1573 X问题(中国剩余定理非互质版本)
- HDOJ 1573 X问题 (余数不互质的中国剩余定理)
- HDOJ 1573 X问题 (中国剩余定理CRT)
- HDOJ-1211 RS(数论)
- HDOJ-1576(数论)
- hdoj-5120-Intersection【数论】
- HDOJ-2669 Romantic(数论)
- HDU/HDOJ 1573 X问题 非互质情况下的中国剩余定理
- HDOJ 题目1573 X问题(中国剩余定理,拓展欧几里得)
- hdoj 1573 X问题 【CRT】【求解方程组在 N内正整数解的个数】
- [杂谈]余刚21年前"内心独白":中南海不是世外桃源(全文)
- python tornado window 平台安装
- One-Based Arithmetic - CF 440C 递推
- ios 循环结构
- 盘点近几年的有名黑客行为
- 【数论】 HDOJ 1573 X问题
- 高性能动画!HTML5 Canvas JavaScript框架KineticJS
- Android Wear OS
- mysql表结构复制
- jmeter使用
- ios for循环实现九九乘法方程式
- linux最近的学习
- ios 数据类型
- C# Media Player 控件问题处理