2019/10/26 TZOJ

1001 Flooded Island 

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=4521

把陆地四面其中三面被海洋环绕的‘@’标记。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
char ma[15][15];
int dir[4][2]={-1,0,1,0,0,-1,0,1};
int main()
{
    memset(ma,'.',sizeof(ma));
    int n,m;scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        getchar();
        for(int j=1;j<=m;j++)
            scanf("%c",&ma[i][j]);
    }
    int mini=30,maxi=-1,minj=30,maxj=-1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            int cnt=0;
            if(ma[i][j]=='X')
            {
                for(int k=0;k<4;k++)
                {
                    int nowi=i+dir[k][0],nowj=j+dir[k][1];
                    if(ma[nowi][nowj]=='.') cnt++;
                }
                if(cnt>=3) ma[i][j]='@';
            }
            if(ma[i][j]=='X')
            {
                    mini=min(mini,i);
                    maxi=max(maxi,i);
                    minj=min(minj,j);
                    maxj=max(maxj,j);
            }
        }
    //printf("%d %d %d %d
",mini,maxi,minj,maxj);
    for(int i=mini;i<=maxi;i++)
    {
        for(int j=minj;j<=maxj;j++)
            printf("%c",ma[i][j]=='X'?'X':'.');
        printf("
");
    }
}

1002 Buy Tickets

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6018

线段树，我先跑了

1003 Strange fuction

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6052

F(x) = 6 * x^7+8*x^6+7*x^3+5*x^2-y*x (0 <= x <=100)，给一个y解出F(x)min(x∈[0,100])

求导，二分。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
double f(double x,double y)
{
    return 42*pow(x,6.0)+48*pow(x,5.0)+21*pow(x,2.0)+10*x-y;
}
double ans(double x,double y)
{
    return 6*pow(x,7.0)+8*pow(x,6.0)+7*pow(x,3.0)+5*x*x-y*x;
}
int main()
{
    int t;scanf("%d",&t);
    while(t--)
        {
            double y;scanf("%lf",&y);
            double mid,left=0.0,right=100.0;
            while(right-left>1e-8)
            {
                mid=(left+right)/2.0;
                if(f(mid,y)<1e-8) left=mid;
                else right=mid;
            }
            printf("%.4f
",ans(mid,y));
        }
}

太难了。自动忽略1004。

1005 A sequence of numbers

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6063

看a、b、c是等差数列还是等比数列，然后算第k个。等差数列没话说，等比数列快速幂，注意mod

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=200907;
ll f(ll a,ll b)
{
    ll res=1;
    while(b)
    {
        if(b&1) res=(res*a)%mod;
        b>>=1;
        a=(a*a)%mod;
    }
    return res;
}
int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
        ll a,b,c,k,ans;scanf("%lld%lld%lld%lld",&a,&b,&c,&k);
        if((a+c)==2*b)
        {
            ll d=(b-a)%mod;
            ans=(d*(k-1)+a)%mod;
        }
        else{
            ll q=(b/a)%mod;
            ans=(f(q,k-1)*a%mod)%mod;
        }
        printf("%lld
",ans);
    }
}

1006 GCD

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6066

求满足__gcd(i,n)>=k的i的个数，将满足的因子欧拉函数累加

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int Eular(int n)
{
    int res=n;
    for(int i=2;i*i<=n;i++)
        if(n%i==0)
        {
            res-=res/i;
            while(n%i==0) n/=i;
        }
    if(n>1) res-=res/n;
    return res;
}
int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
        ll ans=0;
        int i,n,k;scanf("%d%d",&n,&k);
        for(i=1;i*i<n;i++)
        {
            if(n%i==0)
            {
                if(i>=k) ans+=Eular(n/i);
                if(n/i>=k) ans+=Eular(i);
            }
        }
        if(i*i==n&&i>=k) ans+=Eular(i);
        printf("%lld
",ans);
    }
}

1007 Following Orders

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6080

我又双叒叕跑了

1008 ST Task

角速度v来转这个杆子，第i秒的影子Lcos(v*PI/180.0*i)，注意向下取整。

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=2632

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define PI acos(-1)
int main()
{
    int l,v;scanf("%d%d",&l,&v);
    for(int i=0;i<=180/v;i++)
    {
        double l1=abs(l*cos(v*PI/180.0*i));
        int ll=(int)(l1+0.5);
        if(v*i<=90)
        {
            for(int j=l-ll;j>=1;j--) printf("S");
            for(int j=ll;j>=1;j--) printf("T");
            for(int j=0;j<l;j++) printf("S");
            printf("  %.5f
",l1*1.0/l*1.0/2.0);
        }
        else{
            for(int j=0;j<l;j++) printf("S");
            for(int j=ll;j>=1;j--) printf("T");
            for(int j=l-ll;j>=1;j--) printf("S");
            printf("  %.5f
",l1*1.0/l*1.0/2.0);
        }
    }
}