T1 首先筛法暴力出1到30w中具有公因数最多的数量(算出只有180),然后我们惊喜地发现我们可以直接做,最小表示法枚举第k位,然后向前判断是否合法,处理处最小的能放进去的就可以了,同时判断不可能是看这一位能否放进去,然而数据中并没有impossible
T2 这题暴力的时间复杂度为O(n^2*2^k),理论上极限数据即使给6s都跑不过,但是我们发现这是随机数据,于是我们可以水过而且很快【思考熊,具体就是状压二进制表示选择了哪些然后暴力验证
T3 原型POJ2323,递推,但是我们要写高精度,公式见代码
http://pan.baidu.com/s/1eQm8gIQ 内有题目和题解以及代码,以下是个人代码
T1
#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#define N 300005
#define M 180
using namespace std;
int p[M+5],n,f=0;
char a[N];
bool v[27];
int main()
{
freopen("arilang.in","r",stdin);
freopen("arilang.out","w",stdout);
scanf("%d\n",&n);memset(p,0,sizeof(p));
for (int i=1;i<=n/2+1;i++)
if (n%i==0) p[++p[0]]=i;
for (int i=1;i<=n;i++)
{
memset(v,1,sizeof(v));
for (int j=1;j<=p[0];j++)
{
if (p[j]<i)
{
//printf("%d ",a[i-p[j]]-'a'+1);
if (i+p[j]>n)
{
v[a[i+p[j]-n]-'a'+1]=0;
}
v[a[i-p[j]]-'a'+1]=0;
}
else break;
}
char l='a';int j;
//for (int j=1;j<=26;j++) cout<<v[j];cout<<' '<<j<<endl;
for (j=1;j<=26;j++)
if (v[j]!=0) break;
else ++l;
a[i]=l;
if (j==27) {f=1;break;}
}
if (f==1) printf("Impossible");
else
for (int i=1;i<=n;i++)
printf("%c",a[i]);
return 0;
}
T2
#include <iostream>
#include <cstdio>
#include <cstring>
#define N 35
#define M 25
#define inf 1000000007
using namespace std;
int g[N][N];
int v[N],q[40];
int ans,anss,K,n,m;
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
bool bfs(int x)
{
memset(v,0,sizeof(v));
int h=0,t=1;q[t]=0;
while (h<t)
{
int p=q[++h];
for (int i=1;i<n;i++)
if (!v[i]&&(g[p][i]&x))
{
v[i]=1;
q[++t]=i;
}
}
return v[1];
}
int main()
{
freopen("license.in","r",stdin);
freopen("license.out","w",stdout);
K=read(),n=read(),m=read();ans=inf;
for (int i=1;i<=m;i++)
{
int x=read(),y=read(),w=read();
g[x][y]=g[y][x]|=1<<w;
}
int x;
for (int i=1;i<1<<K;i++)
{
x=0;
for (int j=i;j;j>>=1)
{
x+=j&1;
}
if (x<ans&&bfs(i))
{
ans=x;
anss=i;
}
}
printf("%d\n",ans);
int k=0;
if (ans!=inf)
{
for (;anss;anss>>=1,k++)
if (anss&1) printf("%d ",k);
}
else printf("Impossible");
return 0;
}
T3
#include<iostream>
#include<cstdio>
#include<string.h>
using namespace std;
const int MAXT=100+9,MAXN=50+9,MAXM=1500+9,E=100000000;
int n,kk;
struct bign
{
int a[9],l;
void operator = (int r)
{
for(l=0;r;r/=E)
a[l++]=r%E;
}
void operator += (const bign &r)
{
int i;
for(i=0;i<r.l;++i)
a[i]+=r.a[i];
if(r.l>l)
l=r.l;
for(i=0;i<l;++i)
if(a[i]>=E)
{
a[i]-=E;
++a[i+1];
}
if(a[l])
++l;
}
}dp[MAXN][MAXM];
void print(bign x)
{
if (x.l)
{
printf("%d",x.a[x.l-1]);
for(int i=x.l-2;i>=0;--i)
printf("%08d",x.a[i]);
putchar(10);
}
else puts("0");
}
int main()
{
freopen("inversion.in","r",stdin);
freopen("inversion.out","w",stdout);
/*dp[1][0]=1;
for(int i=1;i<51;i++)
for(int j=0;dp[i][j];j++)
for(int k=0;k<=i;k++)
dp[i+1][j+k]+=dp[i][j];*/
dp[0][0]=1;dp[1][0]=1;dp[2][0]=1;dp[2][1]=1;
for (int i=3;i<=50;i++)
{
dp[i][0]=1;
for (int j=1;j<=i*(i-1)/2;j++)
{
dp[i][j]=0;
for (int k=1;k<=i;k++)
dp[i][j]+=dp[i-1][j-k+1];
}
}
for(int t=1;scanf("%d%d",&n,&kk),n;++t)
{
print(dp[n][kk]);
}
}