[SGU200 Cracking RSA]

最新推荐文章于 2019-01-08 08:27:47 发布

weixin_30326745

最新推荐文章于 2019-01-08 08:27:47 发布

阅读量68

点赞数

原文链接：http://www.cnblogs.com/procedure2012/archive/2012/03/22/2412106.html

版权

[关键字]：数学高斯消元异或方程组

[题目大意]：给出m个数，保证每个数的质因数都不超过第t大的质数。从中选出若干数使其乘积为完全平方数，问有多少选择方法。

//=================================================================================

[分析]：首先可以想到如果一个数是完全平方数那它的所有质因数的指数都应该是偶数。根据这个可以列出方程组：

a11*x1+a12*x2+……+a1n*xn=d1%2=0

a21*x1+a22*x2+……+a2n*xn=d2%2=0

……

at1*x1+at2*x2+……+atn*xn=d%t=0

其中xi={0,1}代表选不选，aij={0,1}就是奇偶性，然后只要求出这个异或（mod 2=xor 1）方程组有多少组解就行了。又因为每个未知数只有{0,1}两种答案，所以2的所有不定元（就是可以任意取值）的个数次方-1（非空集）。

[代码]：

View Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

int sum,t,m,ans,line;
int p[110],g[210],d[110][110];
bool v[600];

void Prepare()
{
     memset(v,0,sizeof(v));
for (int i=2;i<=1000;++i)
     {
if (!v[i]) p[++sum]=i;
if (sum==100) break;
for (int j=1;(j<=sum && i*p[j]<=1000);++j)
         {
             v[i*p[j]]=1;
if (i%p[j]==0) break;
         }
     }
//printf("%d\n",p[sum]);
}

void Init()
{
     scanf("%d%d",&t,&m);
int x;
for (int i=1;i<=m;++i)
     {
         scanf("%d",&x);
for (int j=1;j<=t;++j)
while (x%p[j]==0)
             {
                   d[j][i]^=1;
                   x/=p[j];
             }
     }
}

void Add(int *a,int *b,int *c)
{
int temp=0,l=max(a[0],b[0]);
for (int i=1;i<=l;++i)
     {
         c[i]=a[i]+b[i]+temp;
         temp=c[i]/10;
         c[i]=c[i]%10;
     }
if (temp>0) c[++l]=temp;
     c[0]=l;
}

void Solve()
{
if (t<m) t=m;
for (int i=1;i<=m;++i)
     {
         ++line;
for (int j=line;j<=t;++j)
if (d[j][i]==1)
             {
for (int k=1;k<=m;++k) swap(d[line][k],d[j][k]);
break;
             }
if (d[line][i]==0) ++ans,--line; else
         {
for (int j=line+1;j<=t;++j)
if (d[j][i]==1)
for (int k=1;k<=m;++k)
                                       d[j][k]^=d[line][k];
         }
     }
     g[0]=1,g[1]=1;
for (int i=1;i<=ans;++i) Add(g,g,g);
     --g[1];
for (int i=g[0];i>=1;--i) printf("%d",g[i]);
     printf("\n");
}

int main()
{
    Prepare();
    Init();
    Solve();
    system("pause");
return 0;
}