• 【Vijos】【1164】曹冲养猪


    中国剩余定理

      没啥重要的……模板题,中国剩余定理就是解出模线性方程组的一个可行解(好像也是唯一解?)

      这是一种神奇的构造方法……明白了为什么这样构造是对的就行了=。=至于怎么想到这种构造方法的……去问孙子去→_→

     1 //Vijos 1164
     2 #include<cstdio>
     3 #include<cstdlib>
     4 #include<cstring>
     5 #include<iostream>
     6 #include<algorithm>
     7 #define rep(i,n) for(int i=0;i<n;++i)
     8 #define F(i,j,n) for(int i=j;i<=n;++i)
     9 #define D(i,j,n) for(int i=j;i>=n;--i)
    10 using namespace std;
    11 
    12 int getint(){
    13     int v=0,sign=1; char ch=getchar();
    14     while(ch<'0'||ch>'9') {if (ch=='-') sign=-1; ch=getchar();}
    15     while(ch>='0'&&ch<='9') {v=v*10+ch-'0'; ch=getchar();}
    16     return v*=sign;
    17 }
    18 typedef long long LL;
    19 /*******************tamplate********************/
    20 
    21 void exgcd(LL a,LL b,LL &d,LL &x,LL &y){
    22     if(!b){ d=a; x=1; y=0; }
    23     else { exgcd(b,a%b,d,y,x); y-=x*(a/b); }
    24 }
    25 LL china(int n,int *a,int *m){
    26     LL M=1,d,y,x=0;
    27     F(i,1,n) M*=m[i];
    28     F(i,1,n){
    29         LL w=M/m[i];
    30         exgcd(m[i],w,d,d,y);
    31         x=(x+y*w*a[i])%M;
    32     }
    33     return (x+M)%M;
    34 }
    35 int a[11],b[11];
    36 int main(){
    37     int n=getint();
    38     F(i,1,n) a[i]=getint(),b[i]=getint();
    39     printf("%lld
    ",china(n,b,a));
    40     return 0;
    41 }
    View Code
  • 相关阅读:
    以正确的方式开源 Python 项目
    一个备胎的自我修养
    关于我们 | 读书马上
    基于libevent, libuv和android Looper不断演进socket编程
    libuv 与 libev 的对比
    OCaml Language Sucks
    Practical Common Lisp
    learning
    WebApi系列~QQ互联的引入(QConnectSDK)
    知方可补不足~用xsl来修饰xml
  • 原文地址:https://www.cnblogs.com/Tunix/p/4269016.html
Copyright © 2020-2023  润新知