【JZOJ3387】【NOIP2013模拟】【NEERC2009(POJ3872)】终极武器（区间问题+二分+完全图）

Problem

  给定N（≤10000）段闭区间[ai，bi]，满足1<=b[i-1]＜ai<=bi<=10^18。对于两个数字p和q，如果对于这N个区间内的任意一个整数num，把它在十进制表示下的后k（≤19）位中某一位上的p换成q（或者q换成p），都满足得到的整数仍然在这N个波段内，数字p和q是k等价的。我们称两两之间k等价的数字组成一个k等价类。
  求1~9这9个数字分成了哪些k等价类。

Solution

• 首先要喷一喷OJ，一开始的题竟然是残缺的。。。
• 然后，改回来之后，这道题给我的感觉依然是虚无缥缈。感觉无从下手。
• 不过，这依然不失为一道好题。它并不困难，但是考验逻辑思维的清晰性；并且，它还让我懂得了——做题时一定要保证思维清晰。

---

• 下面切入正题。
• 首先，我们考虑枚举一个区间i。枚举后k位中的第j位，并将ai、bi分别分成三部分。比如对于ai=2345，j=3（倒数第3位，也即正数第2位），我们会如下划分：
  
• 如图，x1（第j位以前）=2，x2（第j位上的数）=3，x3（第j位以后）=45。bi的划分同理。
• 这样，我们就可以进行分类讨论。

---

• 情况1：当ai和bi的x1、x2均相同时，就要看x3的取值。譬如当ai=2345，bi=2389时，x3的取值范围为[45,89]。考虑暴力枚举第j位修改成什么，设修改成l，则[2345,2389]中的数均可被修改成[2l45,2l89]。
• 我们遇到了一个非常棘手的问题——判断[2l45,2l89]这段区间内的每一个数是否都在n个区间中出现过。

---

• 首先，题目给定了b[i-1]＜ai<=bi；也就是说，当前的区间肯定在前一个区间右边，并且所有区间两两之间的交集为空集。
• 考虑将所有首尾相连的区间合并。这样，任意两个相邻的区间之间都会有至少一个空格。
• 那么，我们如要判断区间[l,r]中的每一个数是否都在n个区间中出现过，只需判断l和r是否都在这n个区间中，并且它们是否同处一个区间。因为如果它们在不同区间，它们之间肯定会有空格，那么它们之前的某些数就没有在这n个中出现过。
• 那么，如何求解一个数出现在哪个区间呢？这个其实很简单，二分一下就好了。
• 至此，情况1被完美解决。

---

• 情况2：当ai和bi的x1相同，但x2不同时，x2的取值也会影响到x3的取值范围。譬如当ai=2345，bi=2689时，x2可以=3,4,5,6。当x2=3时，x3=[45,99]；当x2=4~5时，x3=[0,99]；当x3=6时，x3=[0,89]。
• 考虑将x2的取值范围劈成三部分：第一部分是ai的x2，此时x3为[ai的x3,99]；第二部分是ai的x2+1~bi的x2-1，此时x3为[0,99]；第三部分是bi的x2，此时x3为[0,bi的x3]。
• 于是可以分三部分处理。

---

• 情况3：当ai和bi的x1不同时，x1的取值会影响到x2、x3的取值。譬如当ai=2345，bi=6789时，x1可以为2,3,4,5,6。这时问题就看似变得有些复杂了。
• 设x1和x2并起来为x4，则x4的取值范围为[23,67]。这看似不太好搞。
• 但当x4=[30,59]时，x3的取值范围都是[0,99]；不论你将x4变成什么，比如你将34修改成39，它对应的x3的取值范围始终是[0,99]。所以你可以不用管x4=[30,59]这种情况。
• 当x4=[23,29]时，类似情况2，不过我们可将x2的取值范围只劈成2部分：第一部分，x2=3，此时x4=[45,99]；第二部分，x2=4~9，此时x4=[0,99]。
• 当x4=[60,67]时，同上述方法，也是将x2的取值范围劈成2部分：第一部分，x2=0~6，此时x4=[0,99]；第二部分，x2=7，此时x4=[0,89]。

---

• 这样，我们就可以处理出哪些数字不属于同一个等价类。
• 考虑维护一个邻接矩阵e，e[i][j]表示数字i与数字j是否等价。初始化为1。
• 每次我们做的时候，枚举一个区间i，并枚举一个位置j，设位置j上的数为x2，若发现将它变成某个y会出锅，就令e[x2][y]=e[y][x2]=0。
• 然后，我们发现，对于每一个等价类，个中的任意两个数字x,y，都必须满足e[x][y]=1。也就是说，它们两两之间必须要右边。那么可以转化为求解所有完全图。
• 一个简洁的dfs即可解决问题。

---

• 时间复杂度：$O(n*k*9^2)$。

Code

#include <cstdio>
#include <cstring>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;

const int N=10001;
int i,j,k,l,n,K,x,y;
ll a[N],b[N],pow,L,R,prx,pry,px,py,ax,ay,bx,by;
bool e[10][10],bz[10],now[10];

int getpos(ll x)
{
    int l=1,r=n;
    while(l<r)
    {
        int mid=l+r>>1;
        if(x<=b[mid])
                r=mid;
        else    l=mid+1;
    }
    return a[l]<=x&&x<=b[l] ? l : 0;
}

void calc(int k,int l)
{
    x=getpos(px); y=getpos(py);
    e[k][l]=e[l][k]= x && y && x==y;
    if(!e[1][9])
    {
        k=-1; l=3;
    }
}

void dfs(int x)
{
    bz[x]=1; printf("%d",x);
    int y;
    fo(y,1,9) now[y]&=e[x][y];
    fo(y,1,9) if( !bz[y] && now[y] ) dfs(y);
}

int main()
{
    scanf("%d%d",&n,&K);
    fo(i,1,n) 
    {
        scanf("%lld%lld",&a[i],&b[i]);
        if( i>1 && a[i]==b[i-1]+1 ) b[--i]=b[i+1], n--;
    }

    memset(e,1,sizeof e);
    fo(i,1,n)
    {
        pow=1;
        fo(j,1,K)
        {
            L=a[i]/pow%10;//ai的x2
            R=b[i]/pow%10;//bi的x2
            prx=a[i]/(pow*10)*(pow*10);//ai的x1
            pry=b[i]/(pow*10)*(pow*10);//bi的x1

            if(prx==pry)
            {
                fo(k,L+1,R-1)
                    fo(l,1,9) 
                        if( k!=l && e[k][l] )
                        {
                            px=prx+1ll*l*pow;
                            py=pry+1ll*(l+1)*pow-1;
                            calc(k,l);
                        }
            }
            else
            {
                fo(k,L+1,9)
                    fo(l,1,9)
                        if( k!=l && e[k][l] )
                        {
                            px=prx+1ll*l*pow;
                            py=prx+1ll*(l+1)*pow-1;
                            calc(k,l);
                        }
                fo(k,1,R-1)
                    fo(l,1,9)
                        if( k!=l && e[k][l] )
                        {
                            px=pry+1ll*l*pow;
                            py=pry+1ll*(l+1)*pow-1;
                            calc(k,l);
                        }
            }

            ax=a[i]%pow; ay=pow-1;
            by=b[i]%pow; bx=0;
            if( prx==pry && L==R ) ay=by, bx=ax;
            fo(k,1,9)
            {
                if( k!=L && e[L][k])
                {
                    px=prx+k*pow+ax; 
                    py=prx+k*pow+ay;
                    calc(L,k);
                }
                if( k!=R && e[R][k])
                {
                    px=pry+k*pow+bx;
                    py=pry+k*pow+by;
                    calc(R,k);
                } 
            }

            pow*=10;
        }
    }

    fo(i,1,9) 
        if(!bz[i]) 
        {
            memset(now,1,sizeof now);
            dfs(i);
            putchar('\n');
        }
}