大步小步算法（网格涂色，uva 11916）

算法详解

http://blog.miskcoo.com/2015/05/discrete-logarithm-problem

参考代码

http://blog.csdn.net/u011345136/article/details/38796399

M至少应当等于不能涂色的格子的行编号的最大值。。。好咯。。。就算这样也还是不会做，还有太多知识点不知道了。

主要问题就是求解一个模方程

a^x=b(mod m)

的最小解。

要用到大步小步算法。

就是令          x=A*sqrt(m)+B。

原方程就变成了          a^(A*sqrt(m)+B)=b(mod m)          或者          a^(A*sqrt(m))*a^(B)=b(mod m)

然后把a^(A*sqrt(m))对m求一下逆，得到v。

方程两边同时乘以v，得到          a^B=b*v(mod m)

枚举所有可能的a^B，把对应的最小的B存到map里。

然后从小到大枚举A，计算出b*v然后在表里找一下有没有a^B，如果有那么A*sqrt(m)+B就是最小解，直接返回。否则就继续找，直到枚举完，返回-1。

关于枚举范围：

首先题目给出条件1<=M,N<=1e8。那么也就是说1<=x<=1e8。而x=x=A*sqrt(m)+B。容易推出1<=A,B<=sqrt(1e8)。这就是A,B枚举的范围。

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=100000007;
typedef pair<ll,ll> pii;

//
//hash
const ll hashsize=sqrt(mod+0.5)+10;
const ll maxstate=sqrt(mod+0.5)+10;
struct Node
{
    ll x,y,next;
}node[maxstate];
ll head[hashsize],tot;
void init()
{
    tot=0;
    memset(head,-1,sizeof(head));
}
void add(ll x,ll y)
{
    ll h=x%hashsize;
    for(ll i=head[h];i!=-1;i=node[i].next)
        if(node[i].x==x)
        {
            node[i].y=min(node[i].y,y);
            return;
        }
    node[tot].x=x;
    node[tot].y=y;
    node[tot].next=head[h];
    head[h]=tot++;
}
ll fd(ll x)
{
    ll h=x%hashsize;
    for(ll i=head[h];i!=-1;i=node[i].next)
        if(node[i].x==x)
            return node[i].y;
    return -1;
}
/
/
//运算专题
ll mul(ll a,ll b)
{
    return (a%mod)*(b%mod)%mod;
}
ll sum(ll a,ll b)
{
    return ((a%mod)+(b%mod))%mod;
}
ll mp(ll x,ll n)
{
    ll ret=1;
    while(n)
    {
        if(n&1) ret=(ret*x)%mod;
        x=(x*x)%mod;
        n>>=1;
    }
    return ret;
}
void ex_gcd(ll a,ll b,ll& d,ll& x,ll& y)
{
    if(!b) {d=a;x=1;y=0;}
    else {ex_gcd(b,a%b,d,y,x);y-=x*(a/b);}
}
ll inv(ll a)
{
    ll d,x,y;
    ex_gcd(a,mod,d,x,y);
    return d==1?(x+mod)%mod:-1;
}
ll china(ll n,ll* a,ll* m)
{
    ll M=1,d,y,x=0;
    for(ll i=0;i<n;i++)
        M*=m[i];
    for(ll i=0;i<n;i++)
    {
        ll w=M/m[i];
        ex_gcd(m[i],w,d,d,y);
        x=(x+y*w*a[i])%M;
    }
    return (x+M)%M;
}
/
//大步小步算法
ll BSGS(ll a,ll b)
{
    ll m=sqrt(mod+0.5);
    init();
    ll e=1;
    ll v=inv(mp(a,m));
    for(ll i=0;i<m;i++)
    {
        add(e,i);
        e=(e*a)%mod;
    }
    for(ll i=0;i<m;i++)
    {
        ll ans=fd(b);
        if(ans!=-1) return i*m+ans;
        b=(b*v)%mod;
    }
    return -1;
}
/
ll N,M,K,B,R,kase;
set<pii>S;
ll solve()
{
    ll cnt=N;
    ll c=0;
    for(set<pii>::iterator it=S.begin();it!=S.end();++it)
    {
        if(it->first!=M&&!S.count(make_pair(it->first+1,it->second)))
            cnt++;
        if(it->first==1)
            cnt--;
        if(it->first==M) c++;
    }
    ll ret=(mul(mp(K,cnt),mp(K-1,M*N-cnt-B)))%mod;
    if(ret==R) return M;
    M++;
    ret=mul(ret,mul(mp(K,c),mp(K-1,N-c)));
    if(ret==R) return M;
    ll v=inv(ret);
    ll cheng=mp(K-1,N);
    return BSGS(cheng,(R*v)%mod)+M;
}

int main()
{
    ll T,x,y;
    scanf("%lld",&T);
    while(T--)
    {
        M=1;
        S.clear();
        scanf("%lld %lld %lld %lld",&N,&K,&B,&R);
        for(ll i=0;i<B;i++)
        {
            scanf("%lld %lld",&x,&y);
            S.insert(make_pair(x,y));
            M=max(M,x);
        }
        printf("Case %lld: %lld\n",++kase,solve());
    }
    return 0;
}