[数论 DAG最长路] Codeforces 772C VK Cup 2017 - Round 2 C. Vulnerable Kerbals

最新推荐文章于 2019-10-02 10:47:15 发布

里阿奴摩西

最新推荐文章于 2019-10-02 10:47:15 发布

阅读量517

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/u014609452/article/details/73196068

版权

数论专栏收录该内容

52 篇文章 0 订阅

订阅专栏

$ax\equiv b \pmod m$ 有解就是 $(a,m) | b$ 也就是 $(a,m)|(b,m)$
那么把数按照和 m <script type="math/tex" id="MathJax-Element-110">m</script>的gcd的整除关系建成DAG 跑一个最长链

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#define pb push_back 
using namespace std;
typedef long long ll;

const int N=200005;

int n,m;
int a[N],vst[N];
vector<int> b[N];
int f[N],g[N],pre[N];

#define read(x) scanf("%d",&(x))
int pnt,ans[N];

typedef pair<int,int> abcd;

pair<int,int> Gcd(int a,int b){
  if (!b) return abcd(1,0);
  abcd t=Gcd(b,a%b);
  return abcd(t.second,t.first-a/b*t.second);
}

inline int solve(ll a,ll b){
  abcd t=Gcd(a,m); int d=t.first*a+t.second*m;
  return ((b/d)*(t.first%m+m)%m)%m;
}

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(n); read(m);
  for (int i=1;i<=n;i++) read(a[i]),vst[a[i]]=1;
  for (int i=0,d;i<m;i++)
    if (!vst[i])
      d=__gcd(i,m),g[d%m]++,f[d%m]++,b[d%m].pb(i);
  for (int j=0;j<m;j++) pre[j]=-1;
  for (int i=1;i<m;i++)
    for (int j=i+i;j<=m;j+=i)
      if (f[j%m]<g[j%m]+f[i])
    f[j%m]=max(f[j%m],g[j%m]+f[i]),pre[j%m]=i;
  int k=0;
  for (int j=0;j<m;j++) if (!vst[j] && f[j]>f[k]) k=j;
  while (k!=-1){
    for (int x:b[k]) ans[++pnt]=x;
    k=pre[k];
  }
  reverse(ans+1,ans+pnt+1);
  for (int i=pnt;i>1;i--)
    ans[i]=solve(ans[i-1],ans[i]);
  printf("%d\n",pnt);
  for (int i=1;i<=pnt;i++)
    printf("%d ",ans[i]);
  return 0;
}