【模拟】NEERC15 G Generators(2015-2016 ACM-ICPC)(Codeforces GYM 100851)

题目链接：

　　http://codeforces.com/gym/100851

题目大意：

　　n个序列。每个序列有4个值x,a,b,c,之后按照x=(a*x+b)%c扩展无穷项。

　　求每个序列各取一个数之后求和不是K的倍数的最大值。

　　（x,a,b,c<=1000,n<=10000,K<=10⁹）

题目思路：

　　【模拟】

　　先暴力把每个序列能够获得的值都求出来。存下最大的两个%K不相等的值。

　　接下来先取每个序列最大的值，如果%K不为0则为答案。

　　否则把其中一个换成次优值。因为前面满足%K不相等所以只用换一个。

 //
 //by coolxxx
 //#include<bits/stdc++.h>
 #include<iostream>
 #include<algorithm>
 #include<string>
 #include<iomanip>
 #include<map>
 #include<stack>
 #include<queue>
 #include<set>
 #include<bitset>
 #include<memory.h>
 #include<time.h>
 #include<stdio.h>
 #include<stdlib.h>
 #include<string.h>
 //#include<stdbool.h>
 #include<math.h>
 #define min(a,b) ((a)<(b)?(a):(b))
 #define max(a,b) ((a)>(b)?(a):(b))
 #define abs(a) ((a)>0?(a):(-(a)))
 #define lowbit(a) (a&(-a))
 #define sqr(a) ((a)*(a))
 #define swap(a,b) ((a)^=(b),(b)^=(a),(a)^=(b))
 #define mem(a,b) memset(a,b,sizeof(a))
 #define eps (1e-8)
 #define J 10
 #define mod 1000000007
 #define MAX 0x7f7f7f7f
 #define PI 3.14159265358979323
 #define N 10004
 #define M 1004
 using namespace std;
 typedef long long LL;
 int cas,cass;
 int n,m,lll,ans;
 struct xxx
 {
     int num,c;
 }q[N][];
 int s[N];
 bool cmp1(xxx aa,xxx bb)
 {
     return aa.c>bb.c;
 }
 void print()
 {
     int i,j;
     for(i=;i<=n;i++)
     {
         for(j=;j<=s[i];j++)
         {
             printf("%d:%d  ",q[i][j].c,q[i][j].num);
         }
         puts("");
     }
 }
 int main()
 {
     #ifndef ONLINE_JUDGE
 //    freopen("1.txt","r",stdin);
 //    freopen("2.txt","w",stdout);
     #endif
     int i,j,k;
     int a,b,c,x;
 //    for(scanf("%d",&cass);cass;cass--)
 //    for(scanf("%d",&cas),cass=1;cass<=cas;cass++)
 //    while(~scanf("%s",s+1))
     while(~scanf("%d",&n))
     {
         ans=;
         scanf("%d",&m);
         for(i=;i<=n;i++)
         {
             scanf("%d%d%d%d",&x,&a,&b,&c);
             mem(s,);q[i][].num=q[i][].num=q[i][].c=q[i][].c=-;
             s[x]=;
             for(j=,x=(a*x+b)%c;!s[x];x=(a*x+b)%c,j++)
                 s[x]=j;
             for(x=c-,j=;x>= && j<;x--)
             {
                 if(!s[x])continue;
                 if(x>=q[i][].c)
                 {
                     if(x%m!=q[i][].c%m)q[i][]=q[i][];
                     q[i][].c=x,q[i][].num=s[x];
                 }
                 else if(x>q[i][].c && x%m!=q[i][].c%m)
                     q[i][].c=x,q[i][].num=s[x];
                 if(q[i][].c!=- && q[i][].c!=-)break;
             }
         }
         for(i=;i<=n;i++)ans+=q[i][].c;
         if(ans%m)
         {
             printf("%d\n",ans);
             for(i=;i<=n;i++)
                 printf("%d ",q[i][].num-);
             puts("");
         }
         else
         {
             cas=;q[][].c=MAX;q[][].c=;
             for(i=;i<=n;i++)
             {
                 if(q[i][].c==-)continue;
                 if(q[i][].c-q[i][].c<q[cas][].c-q[cas][].c)
                     cas=i;
             }
             if(cas==){puts("-1");continue;}
             ans=ans-q[cas][].c+q[cas][].c;
             printf("%d\n",ans);
             q[cas][]=q[cas][];
             for(i=;i<=n;i++)
                 printf("%d ",q[i][].num-);
             puts("");
         }
     }
     return ;
 }
 /*
 //

 //
 */