ACM study day3

今天练了二分和快速幂，题目挺难的，挑几个我做上的说一下吧。

先给出几个二分和快速幂的模板函数；

二分

void BS(int m)
{
    int x=,y=a[m-]-a[];
    while(y-x>)
    {
        int mid=(x+y)/;
        if(judge(mid))
            x=mid;
        else
            y=mid;
    }
    cout<<x<<endl;
}

快速幂

LL PowerMod(LL a, LL b)//a^b
{
    LL ans = ;
    while (b > )
    {
        if (b & )
        {
            ans = (ans*a) % MOD;
        }
        b >>= ;
        a = (a*a) % MOD;
    }
    return ans;
}  

矩阵乘法加矩阵快速幂取模

struct mat
{
    ll a[][];
};
mat mat_mul(mat x,mat y)
{
    mat res;
    memset(res.a,,sizeof(res.a));
    for(int i=;i<;i++)
        for(int j=;j<;j++)
        for(int k=;k<;k++)
        res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;
    return res;
}  

void mat_pow(int n)
{
mat c,res;
c.a[][]=c.a[][]=c.a[][]=;
c.a[][]=;
memset(res.a,,sizeof(res.a));
for(int i=;i<;i++) res.a[i][i]=;
while(n)
{
if(n&) res=mat_mul(res,c);
c=mat_mul(c,c);
n=n>>;
}
printf("%I64d\n",res.a[][]);
}

A  POJ 2456

这道题就是二分，但是二分检验的地方比较有趣，全部贴上来，那个judge函数我想了一阵子。。。（比较笨）。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <string>
#include <map>
#include <vector>
#include <set>
#include <queue>
using namespace std;
int a[];
int m,n;
bool judge(int d)
{
    int fir,last=;
    for(int i=;i<n-;i++)
    {
    fir=last+;
    while(fir<m&&a[fir]-a[last]<d)
        fir++;
    if(fir==m)
        return false;
    last=fir;
    }
    return true;
}
void BS(int m)
{
    int x=,y=a[m-]-a[];
    while(y-x>)
    {
        int mid=(x+y)/;
        if(judge(mid))
            x=mid;
        else
            y=mid;
    }
    cout<<x<<endl;
}

int main()
{

cin>>m>>n;
for(int i=;i<m;i++)
{
    scanf("%d",&a[i]);
}
sort(a,a+m);
BS(m);
return ;
}