2021年广东工业大学第十五届文远知行杯程序设计竞赛（同步赛 DHG

**

D动态序列

**
传送门
题目大意：
给定一个序列，有n个数，进行如下操作

1 b：序列中所有数乘以整数b

2 b：序列中所有数增加整数b

3 b：在序列头部添加一个整数b

4 b：在序列尾部添加一个整数b

5 p：输出序列的第p个数，并将结果对1000000007取模
1 ≤ b ≤ 1e9， 1 ≤ n ≤ 1e5
思路： 线段树思想
对于乘法和加法的修改只要用懒标记即可，最后输出把结果带上懒标记。
对于添加的数字，我们把它去除前面的懒标记，设加法标记plu , 乘法标记mul ，再用一下乘法逆元，记录的数字即为 x = (x - plu + mod ) %mod * ksm(mul , mod-2 ) %mod。
注意数组开3e5以上，然后输入的时候从下标1e5+1开始存就行。
（一开始想复杂了，稀里哗啦打了一长串线段树，后来一点用没有，还是没掌握主要思想QAQ

#include <iostream>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <vector>
using namespace std;
#define IOS ios::sync_with_stdio(false);cin.tie(0)
#define _for(i,a,b) for(int i=(a) ;i<=(b) ;i++)
#define _rep(i,a,b) for(int i=(a) ;i>=(b) ;i--)
#define mst(x,y) memset(x,y,sizeof(x))
#define all(v) v.begin() , v.end()
#define pb(v) push_back(v)
#define INF 0x3f3f3f3f
#define int long long
#define lson 2*p,l,mid
#define rson 2*p+1,mid+1,r
typedef long long ll;
const int N = 3e5+10;
const int mod =1e9+7;
int n,q;
int a[N];
//a*(x+b) = a*x +a*b
int mul=1;
int plu;
int top=1e5+1,tail;//数组首尾指针
ll qsm( ll a, ll b )
{
    ll ans = 1, temp = a;
    while( b )
    {
        if( b&1 )   ans = (ans * temp ) %mod;
        temp = ( temp * temp ) %mod;
        b>>=1;
    }
    return ans;
}
void solve(int op, ll b)
{
    if( op==1 )
    {
        mul = ( mul * b )%mod;
        plu = ( plu * b) %mod;
    }
    else if( op==2 )
    {
        plu = ( plu + b) %mod;
    }
    else if( op==3 )
    {
        b = (b - plu + mod ) %mod * qsm(mul , mod-2 ) %mod;
        a[--top] = b;
    }
    else if( op==4 )
    {
        b = (b - plu + mod ) %mod * qsm(mul , mod-2 ) %mod;
        a[++tail] = b;
    }
    else
    {
        ll ans = (a[top+b-1] * mul %mod + plu ) %mod;
        cout<<ans<<endl;
    }
}
signed main()
{
//    freopen("data.txt",1"r",stdin);
    IOS;

    cin>>n>>q;
    tail = 1e5+n;
    _for(i,1+1e5,n+1e5)  cin>>a[i];

    while( q-- )
    {
        int x,y;
        cin>>x>>y;
        solve(x,y);
    }
}

**

H H-有多短

传送门

**
思路： 首先明确，最远两个点指经过路径的权值和最大的两个点，
要想使最大值最小值，一般想着贪心二分，不过这里好像没什么好贪好二分的。
因为图是一棵树，最远的两点显然是两个叶子结点，度数为1，将所有k都平均分配到所有叶子结点上，就会让最大值最小，有n个叶子,每个叶子连着的边权值为 1.0 * k / n， 最大值是 2 * 1.0 * k / n ，只会经过两次（起点和终点），平均分配会使得所有的叶子路径都是相同权值和（假设其他非叶子结点的边权为0）。

#include <iostream>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <vector>
using namespace std;
#define IOS ios::sync_with_stdio(false);cin.tie(0)
#define _for(i,a,b) for(int i=(a) ;i<=(b) ;i++)
#define _rep(i,a,b) for(int i=(a) ;i>=(b) ;i--)
#define mst(x,y) memset(x,y,sizeof(x))
#define all(v) v.begin() , v.end()
#define pb(v) push_back(v)
#define INF 0x3f3f3f3f
#define int long long
#define lson 2*p,l,mid
#define rson 2*p+1,mid+1,r
typedef long long ll;

const int N =1e5+10;
const int mod =1e9+7;
int n,k;
int a[N];
signed main()
{
//    freopen("data.txt","r",stdin);
//    IOS;
    int n;
    cin>>n>>k;
    _for(i,1,n-1)
    {
        int x,y;
        cin>>x>>y;
        a[x]++;
        a[y]++;
    }
    int cnt=0;
    _for(i,1,1e5)
    {
        if( a[i] == 1)//统计叶子结点个数
        {
            cnt++;
        }
    }
    printf("%.6f\n",2 * 1.0 * k / cnt);
}

G 分割

传送门

思路：
总面积 = 所有的矩形长的和 * 所有矩形宽的和
对于所有x的和，有如下计算公式

考虑每个xx[k]的贡献，当它作为xx[j]时，它被加了i-1次，当它作为xx[i]时，它被减了n-i次，所以每个xx[i]的贡献为xx[i] * ( i-1 -( n-i) ) ， y同理

#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <vector>
using namespace std;
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
#define _for(i,a,b) for(int i=(a) ;i<=(b) ;i++)
#define _rep(i,a,b) for(int i=(a) ;i>=(b) ;i--)
#define mst(x,y) memset(x,y,sizeof(x))
#define all(v) v.begin() , v.end()
#define pb(v) push_back(v)
#define INF 0x3f3f3f3f
typedef long long ll;

const int N = 1e5+10;
const int mod = 1e9+7;
int n,m;
ll xx[N];
ll yy[N];

signed  main()
{
    IOS;
    cin>>n>>m;
    _for(i,1,n) cin>>xx[i];
    _for(i,1,m) cin>>yy[i];

    sort(xx+1 , xx+1+n);
    sort(yy+1 , yy+1+m);
    ll ans=0;
    ll cnt=0;
    _for(i,1,n)
    {
        ans += xx[i] * ( i-1 -( n-i) )%mod;
        ans %=mod;
    }
    _for(i,1,m)
    {
        cnt += yy[i] * ( i-1 - (m-i) )%mod;
        cnt %=mod;
    }
    cout<<ans*cnt%mod<<endl;
}