每日一练 19.2.16

NOIP2012 提高组

同余方程

题目描述
求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解。

输入
输入只有一行，包含两个正整数 a, b，用一个空格隔开。

输出
输出只有一行，包含一个正整数 x0，即最小正整数解。输入数据保证一定有解。

样例输入
3 10

样例输出
7

提示
对于 100%的数据，2 ≤a, b≤ 2,000,000,000。

思路
通过扩展欧几里得算法实现

扩展欧几里得算法及证明

代码实现

#include <iostream>
#include <cstdio>
#include <string>
#include <queue>
#include <cstring>
#include <cmath>
#include <cctype>
#include <vector>
#include <algorithm>
#include <map>
#include <stack>
 
using namespace std;
typedef long long ll;
typedef pair<int,int> pa;
 
const int N=105;
int a,b,exx,exy;
int ex_gcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    int ans=ex_gcd(b,a%b,x,y);
    int temp=x;
    x=y;
    y=temp-(a/b)*x;
    return ans;
}
int main()
{
    scanf("%d%d",&a,&b);
    ex_gcd(a,b,exx,exy);
    while(exx<0) exx+=b;
    printf("%d\n",exx);
    return 0;
}
 

借教室

题目描述
在大学期间，经常需要租借教室。大到院系举办活动，小到学习小组自习讨论，都需要向学校申请借教室。教室的大小功能不同，借教室人的身份不同，借教室的手续也不一样。
面对海量租借教室的信息，我们自然希望编程解决这个问题。
我们需要处理接下来n天的借教室信息，其中第i天学校有ri个教室可供租借。共有m份订单，每份订单用三个正整数描述，分别为dj, sj, tj ，表示某租借者需要从第sj 天到第tj 天租借教室（包括第sj天和第tj天），每天需要租借dj个教室。
我们假定，租借者对教室的大小、地点没有要求。即对于每份订单，我们只需要每天提供dj个教室，而它们具体是哪些教室，每天是否是相同的教室则不用考虑。
借教室的原则是先到先得，也就是说我们要按照订单的先后顺序依次为每份订单分配教室。如果在分配的过程中遇到一份订单无法完全满足，则需要停止教室的分配，通知当前申请人修改订单。这里的无法满足指从第sj天到第tj天中有至少一天剩余的教室数量不足dj个。
现在我们需要知道，是否会有订单无法完全满足。如果有，需要通知哪一个申请人修改订单。

输入
第一行包含两个正整数n, m，表示天数和订单的数量。
第二行包含n个正整数，其中第i个数为ri，表示第i天可用于租借的教室数量。
接下来有m行，每行包含三个正整数dj, sj, tj ，表示租借的数量，租借开始、结束分别在第几天。
每行相邻的两个数之间均用一个空格隔开。天数与订单均用从1开始的整数编号。

输出
如果所有订单均可满足，则输出只有一行，包含一个整数 0。否则（订单无法完全满足）输出两行，第一行输出一个负整数-1，第二行输出需要修改订单的申请人编号。

样例输入
4 3
2 5 4 3
2 1 3
3 2 4
4 2 4

样例输出
-1
2

提示
第1份订单满足后，4天剩余的教室数分别为0，3，2，3。第2份订单要求第2天到第4天每天提供3个教室，而第3天剩余的教室数为2，因此无法满足。分配停止，通知第2个申请人修改订单。

对于100%的数据，有1≤n,m≤10⁶,0≤ri,dj≤10⁹,1≤sj≤tj≤n。

思路
可以用线段树或二分的算法解决这道问题

代码实现
线段树

#include <iostream>
#include <cstdio>
#include <string>
#include <queue>
#include <cstring>
#include <cmath>
#include <cctype>
#include <vector>
#include <algorithm>
#include <map>
#include <stack>

using namespace std;
typedef long long ll;
typedef pair<int,int> pa;

const int N=1000005;
struct rec
{
    int sum;
    int inc;
} tree[2*N];
bool flag;
int i,n,m;
void build(int l, int r, int rt)
{
    if(l+1 == r)
    {
        int temp;
        scanf("%d",&temp);
        tree[rt].sum=temp;
        return;
    }
    int mid = (l + r)>>1;
    build(l, mid, rt <<1);
    build(mid, r, (rt <<1) + 1);
    tree[rt].sum=min(tree[rt<<1].sum,tree[(rt <<1) + 1].sum);
}
void update(int rt)
{
    if(!tree[rt].inc) return;
    tree[rt<<1].inc+=tree[rt].inc;
    tree[(rt<<1)+1].inc+=tree[rt].inc;
    tree[rt<<1].sum+=tree[rt].inc;
    tree[(rt<<1)+1].sum+=tree[rt].inc;
    tree[rt].inc=0;
}
void modify(int rt,int l,int r,int ll,int rr,int delta)
{
    if(flag) return;
    if(l>=ll && r<=rr)
    {
        tree[rt].sum+=delta;
        if(tree[rt].sum<0)
        {
            printf("-1\n%d\n",i);
            flag=true;
            return;
        }
        tree[rt].inc+=delta;
        return;
    }
    update(rt);
    int mid=(l+r)>>1;
    if(ll<mid) modify(rt<<1,l,mid,ll,rr,delta);
    if(rr>mid) modify((rt<<1)+1,mid,r,ll,rr,delta);
    tree[rt].sum=min(tree[rt<<1].sum,tree[(rt<<1)+1].sum);
}
int main()
{
    scanf("%d%d",&n,&m);
    build(1,n+1,1);
    for(i=1; i<=m; i++)
    {
        int d,s,t;
        scanf("%d%d%d",&d,&s,&t);
        modify(1,1,n+1,s,t+1,-d);
        if(flag) return 0;
    }
    printf("0\n");
    return 0;
}