习题日常第二十练

1恋与程序员（题目链接）

这个题一开始我以为是个最短路，但是老感觉怪怪的，后来发现它并不满足用最短路的条件，因为卡片重复用过后只算一次的价值。直接用dfs深搜，加一点点剪枝。代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define  ll long long int
/*const int mod=100003;*/
using namespace std;
int n,m,k,p,vis[110],a[110][110],z[110],ans;
struct node
{
    int f,w,h;
}t[110][110];
void dfs(int x,int tep)
{
    int i;
    if(tep>=ans)
        return;
    if(x==p)
    {
        ans=min(ans,tep);
        return;
    }
    for(i=1;i<=n;i++)
    {
        if(t[x][i].f)
        {
            if(!vis[t[x][i].h])
            {
                vis[t[x][i].h]=1;
                dfs(i,tep+t[x][i].w);
                vis[t[x][i].h]=0;
            }
            else
            {
                 dfs(i,tep);
            }
        }
    }
}
int main()
{
    ios::sync_with_stdio(false);
    int i,j,x,y,l,r,lr,rr,mod;
    while(cin>>n>>m>>k>>p)
    {
        ans=1e8;
        memset(a,0,sizeof(z));
        memset(vis,0,sizeof(vis));
        memset(z,0,sizeof(z));
        for(i=0;i<=n;i++)
        {
            for(j=0;j<=n;j++)
            {
                t[i][j].f=0;
            }
        }
        for(i=0;i<m;i++)
        {
           cin>>x>>y>>l;
           t[x][y].f=1;
           t[x][y].h=l;
        }
        for(i=0;i<k;i++)
        {
            cin>>x>>y;
            z[x]=y;
        }
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=n;j++)
            {
                if(t[i][j].f)
                {
                    t[i][j].w=z[t[i][j].h];
                }
            }
        }
        dfs(1,0);
        cout<<ans<<endl;
    }
  return 0;
}

2LP 钱不够（题目链接）

这个题就是个简单的动态规划，每一个格子（除了边界）的值都加上它行上的前一个和列上的前一个中的比较小的一个即可。代码如下。

 memset(a,0,sizeof(a));
        cin>>n;
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=n;j++)
                cin>>a[i][j];
        }
        for(i=1;i<=n;i++)
        {
            a[i][1]+=a[i-1][1];
            a[1][i]+=a[1][i-1];
        }
        for(i=2;i<=n;i++)
        {
            for(j=2;j<=n;j++)
            {
                a[i][j]+=min(a[i-1][j],a[i][j-1]);
            }
        }
        cout<<a[n][n]<<endl;

3牛牛的幂运算（题目链接）

数论题，自己不会做，不会的主要原因是想不到如何去重。看过别人的题解以后才发现可以通过拆底数的方法。具体代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define  ll long long int
const int mod=1e9+7;
using namespace std;
ll mypow(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1)
            ans*=a;
        a*=a;
        b>>=1;
    }
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    ll n,j,k,l,i;
    cin>>n;
    ll ans=n*n;
    for(i=1;i<33;i++)
        for(j=1;j<33;j++)
    {
        if(__gcd(i,j)==1)
        {
            k=pow(n,1.0/i);
            l=pow(n,1.0/j);
            for(;mypow(k,i)<=n;k++);
            for(;mypow(l,j)<=n;l++);
            k--;
            l--;
            ans=(ans+(min(k-1,l-1)*min(n/i,n/j)%mod))%mod;
        }
    }
    cout<<ans;
  return 0;
}

4旅行（题目链接）

这个题算是个最短路的板题吧。我们可以通过枚举中转点的方式，找出离中转点最远的两个点作为起始和终止点。最后得出最大值即可。代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define  ll long long int
const int mod=1e9+7;
using namespace std;
ll n,m,k,tot,a[2010],dis[2010],b[2010],sum;
struct node
{
    ll to,next,w;
}t[2010];
void add(ll x,ll y,ll z)
{
    t[++tot].to=y;
    t[tot].w=z;
    t[tot].next=a[x];
    a[x]=tot;
    return;
}
priority_queue<pair<ll,ll> >p;
ll dist(ll x)
{
    ll i;
    for(i=1;i<=n;i++)
        dis[i]=1e8;
        memset(b,0,sizeof(b));
    dis[x]=0;
    p.push(make_pair(0,x));
    while(!p.empty())
    {
        ll xx=p.top().second,i;
        p.pop();
        if(b[xx])
            continue;
        b[xx]=1;
        for(i=a[xx];i!=-1;i=t[i].next)
        {
            if(dis[t[i].to]>dis[xx]+t[i].w)
            {
                dis[t[i].to]=dis[xx]+t[i].w;
                p.push(make_pair(-dis[t[i].to],t[i].to));
            }
        }
    }
    ll ans1=0,ans2=0,flat=0;
    for(i=1;i<=n;i++)
    {
        if(dis[i]<1e8&&dis[i]>ans1)
            ans1=dis[i];
    }
    for(i=1;i<=n;i++)
    {
        if(flat)
            break;
        if(dis[i]==ans1)
        {
            dis[i]=0;
            flat=1;
        }
    }
    for(i=1;i<=n;i++)
    {
         if(dis[i]<1e8&&dis[i]>ans2)
            ans2=dis[i];
    }
    if(ans1!=0&&ans2!=0)
        return ans1+ans2;
    else
        return -1;
}
int main()
{
    ios::sync_with_stdio(false);
    ll j,k,l,i,tz;
    cin>>tz;
    while(tz--)
    {
        tot=0;
        sum=-1;
       cin>>n>>m;
       for(i=0;i<=n;i++)
        a[i]=-1;
       for(i=0;i<m;i++)
       {
           cin>>j>>k>>l;
           add(j,k,l);
           add(k,j,l);
       }
       for(i=1;i<=n;i++)
       {
           sum=max(dist(i),sum);
       }
       cout<<sum<<endl;
    }
  return 0;
}