[NOI2014]随机数生成器

最新推荐文章于 2018-06-04 14:21:40 发布

原创最新推荐文章于 2018-06-04 14:21:40 发布 · 314 阅读

0 ·

CC 4.0 BY-SA版权

NOI 专栏收录该内容

1 篇文章

订阅专栏

本文通过具体实例详细解析了一种矩阵路径构造算法的过程。通过定义矩阵A和路径集合B、未经过路径集合C，逐步推导出必须经过的路径及不可经过的路径集合，最终确定了完整的路径集合。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

直接照着题目构造 $X_i,i\in [1,nm]$


12	9	1	7
5	11	6	2
4	10	3	8

我们把矩阵数组记为 $A$
路径集合 $B$ ，路径不经过的集合 $C$
首先：
$A_{1,3}=1$ 必须经过
那么 $A_{2..3,1..2}(=\{4,5,11,10\})$ 和 $A_{1..0,4..4}(\not O)$ 不能经过
然后 $A_{2,4}=2$ 也必须经过。
$A_{3..3,1..3}(=\{4,10,3\})$ 和 $A_{1..1,5..4}(\not O)$ 不能经过
现在为止：
$B=\{1,2\}\ C=\{3,4,5,10,11\}$
显然， $6$ 也必须经过（ $3,4,5$ 不能经过）
那么 $A_{3..3,1..2}(=\{4,10\})$ 和 $A{1..1,4..4}(=\{7\})$ 不能经过
$B=\{1,2,6\}\ C=\{3,4,5,7,10,11\}$
后面同理，最后得：
$B=\{1,2,6,8,9,12\}\ C=\{3,4,5,7,10,11\}$
所以不就是贪心？？？

ll x,a,b,c,d;
#define __ 5000
#define _ __*__
int n,m,g[_+10],k[_+10],l[__+10],r[__+10];
void swap(int x,int y)
{
    int ___=g[x];
    g[x]=g[y];
    g[y]=___;
}
int main(){
    x=read();
    a=read();
    b=read();
    c=read();
    d=read();
    n=read();
    m=read();
    fr(i,1,n*m)
        g[i]=i;
    fr(i,1,n*m)
    {
        x=(a*x*x+b*x+c)%d;
        int k=x%i+1;
        swap(i,k);
    }
    fr(i,1,read())
    {
        int u=read(),v=read();
        swap(u,v);
    }
    fr(i,1,n*m)
        k[g[i]]=i-1;
    fr(i,1,n)
    {
        l[i]=1;
        r[i]=m;
    }
    fr(i,1,n*m)
        if(k[i]%m+1>=l[k[i]/m+1]&&k[i]%m+1<=r[k[i]/m+1])
        {
            printf("%d ",i);
            fr(j,1,k[i]/m)
                r[j]=min(r[j],k[i]%m+1);
            fr(j,k[i]/m+2,n)
                l[j]=max(l[j],k[i]%m+1);
        }
    rt 0;
}