day 15

我变强了!,今天模拟赛我第11!!!,得分200/300,感谢良心数据!;

以下是总结:

T1贪心,排序后第K大乘第K小即可,long long可过;

T2我不会做,得了三十分暴力分,听说有原题,题面是  1/x+1/y=1/(n!)  ,求有多少种x,y符合题意;

T3其实是二分图匹配,但我不会,写了个贪心+搜索,还好数据没卡我,得了70分;

于是得分100+30+70=200!;

开心开心。

T1

库特的向量(code) 

#include<iostream>
#include<cstdio>
#include<cctype>
#include<algorithm>
using namespace std;
long long a[1010],b[1010],n;
long long ans;
inline int read()
{
    int x=0,f=1;char c=getchar();
    while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
    while(isdigit(c)){x=(x<<3)+(x<<1)+(c^48);c=getchar();}
    return x*f;
}
inline int cmp(int x,int y)
{
    return x<y;
}
inline int cmp1(int x,int y)
{
    return x>y;
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)a[i]=read();
    for(int i=1;i<=n;i++)b[i]=read();
    sort(a+1,a+n+1,cmp);sort(b+1,b+n+1,cmp1);
    for(int i=1;i<=n;i++)
        ans+=a[i]*b[i];
    printf("%lld",ans);
    return 0;
}
View Code

 T2

恭介的法则(rule)

经过我苦思冥想一顿调试,找出了一个方法:

1.分解n!的质因数,由于要分解的是(n!)2

但只需要将分解出的质数指数乘2即可,由约数个数定理可得出答案;

2.但这和原题(bzoj2721)(洛谷1445)不太一样,它不取模!;

于是,我们还需要高精度乘法;而且要压位!不然存不下;

于是,一段好看丑陋的代码产生了

#pragma GCC optimize(2)
#include<iostream>
#include<cstdio>
using namespace std;
const int maxn=3000000;
long long n,num[maxn],ans[3000000]={0,1},size=1;
void dec(int x)
{
    while(x%2==0){
        num[2]++;x/=2;
    }
    for(int i=3;i*i<=x;i+=2)
    while(x%i==0)
    num[i]++,x=x/i;
    if(x!=1) num[x]++;
}
void mul(long long x)
{
    for(int i=1;i<=size;i++)
    ans[i]*=x;
    for(int i=1;i<size;i++)
    {
        if(ans[i]>=10000)
        {
            int temp=ans[i]/10000;
            ans[i]%=10000;
            ans[i+1]+=temp;
        }
    }
    while(ans[size]>=10000)
    {
        int temp=ans[size]/10000;
        ans[size]%=10000;
        ans[++size]+=temp;
    }
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    dec(i);
    for(int i=1;i<=n;i++)
    if(num[i])
    {
        mul(num[i]*2+1);
    }
    printf("%lld",ans[size]);
    for(int i=size-1;i>=1;i--)
    {
        printf("%lld",ans[i]/1000);
        printf("%lld",ans[i]/100%10);
        printf("%lld",ans[i]/10%10);
        printf("%lld",ans[i]%10);
    }
    printf("\n");
    return 0;
}
View Code

但是,它tle了!于是,我继续修改;

经过同宿舍dalao——lx指点;

我AC了!

想法:由于每个num中的值可能很小,所以我们只要累积到一定大小再乘入答案;

这样就可以跑的飞快!

AC代码

#pragma GCC optimize(2)
#include<iostream>
#include<cstdio>
using namespace std;
const int maxn=3000000;
unsigned long long n,num[maxn],ans[3000000]={0,1},size=1;
void dec(int x)
{
    while(x%2==0){
        num[2]++;x/=2;
    }
    for(int i=3;i*i<=x;i+=2)
    while(x%i==0)
    num[i]++,x=x/i;
    if(x!=1) num[x]++;
}
void mul(long long x)
{
    for(int i=1;i<=size;i++)
    ans[i]*=x;
    for(int i=1;i<size;i++)
    {
        if(ans[i]>=10000)
        {
            int temp=ans[i]/10000;
            ans[i]%=10000;
            ans[i+1]+=temp;
        }
    }
    while(ans[size]>=10000)
    {
        int temp=ans[size]/10000;
        ans[size]%=10000;
        ans[++size]+=temp;
    }
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    dec(i);
    unsigned long long p=1;
    for(int i=1;i<=n;i++)
    if(num[i])
    {
        p*=(num[i]*2+1);
        if(p>50000)
        {
            mul(p);
            p=1;
        }
    }
    mul(p);
    printf("%lld",ans[size]);
    for(int i=size-1;i>=1;i--)
    {
        printf("%lld",ans[i]/1000);
        printf("%lld",ans[i]/100%10);
        printf("%lld",ans[i]/10%10);
        printf("%lld",ans[i]%10);
    }
    printf("\n");
    return 0;
}
View Code

T3

沙耶的玩偶(doll) 

#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;
int n,m,r,c,cnt;
vector<int>q[2505];
bool v[2505];
int num[60][60],l[2505];
char a[60][60];
int ans;
void build(int x,int y){q[x].push_back(y);}
bool xyl(int x)
{
    for(int i=0;i<q[x].size();i++)
    {
        if(!v[q[x][i]])
        {
            int j=q[x][i];
            v[j]=1;
            int w=l[j];
            l[j]=x;
            if(!w||xyl(w))return 1;
            l[j]=w;
        }
    }
    return 0;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&r,&c);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
     {
       cin>>a[i][j];
       if(a[i][j]=='.')
        num[i][j]=++cnt; 
     }
    int x,y;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      if(a[i][j]=='.')
      {
        x=i+r;y=j+c;
        if(x>=1&&x<=n&&y>=1&&y<=m&&a[x][y]=='.') build(num[i][j],num[x][y]);
        x=i+r;y=j-c;
        if(x>=1&&x<=n&&y>=1&&y<=m&&a[x][y]=='.') build(num[i][j],num[x][y]);
        x=i+c;y=j+r;
        if(x>=1&&x<=n&&y>=1&&y<=m&&a[x][y]=='.') build(num[i][j],num[x][y]);
        x=i+c;y=j-r;
        if(x>=1&&x<=n&&y>=1&&y<=m&&a[x][y]=='.') build(num[i][j],num[x][y]);
      }
    for(int i=1;i<=cnt;i++)
    {
      memset(v,0,sizeof(v));
      ans+=xyl(i);
    }
    cout<<cnt-ans<<endl;
    return 0;
}
View Code

完。

posted @ 2019-08-15 14:51 Frost_Delay 阅读(...) 评论(...) 编辑 收藏