Codeforces Round #448 (Div. 2)

https://codeforces.com/contest/895/problem/C

题意：
一个集合选出一个子集 他们的乘积是一个平方数
问又多少个这样的集合

题解：
数的范围只有70，可以将其素因子状压一下，平方数只与质因子的奇偶性有关，dp[i][j]表示前i个数j状态的方案数。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pb push_back
#define mem(a,b) memset(a,b,sizeof(a))

const int N = 1e5 + 5;
const int MOD = 1e9 + 7;
int a[N];
int cnt[75];
int dp[75][ ( 1 << 19 ) + 5];
int s[75];
int prime[19] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67};
int _2p[N];
int main()
{
    ios::sync_with_stdio ( false );
    cin.tie ( 0 );
    for ( int i = 1; i <= 70; i++ )
    {
        int t = i;
        for ( int j = 0; j < 19; j++ )
        {
            while ( t % prime[j] == 0 )
                t /= prime[j], s[i] ^= ( 1 << j );
        }
    }
    _2p[0] = 1;
    for ( int i = 1; i < N; i++ )
        _2p[i] = ( _2p[i - 1] * 2 ) % MOD;
    int n;
    cin >> n;
    for ( int i = 1; i <= n; i++ )
        cin >> a[i], cnt[a[i]]++;
    dp[0][0] = 1;
    for ( int i = 1; i <= 70; i++ )
    {
        if ( !cnt[i] )
        {
            for ( int j = 0; j < ( 1 << 19 ); j++ )
                dp[i][j] = dp[i - 1][j];
        }
        else
        {
            for ( int j = 0; j < ( 1 << 19 ); j++ )
            {
                dp[i][j ^ s[i]] = ( ( ll ) dp[i][j ^ s[i]] + ( ll ) _2p[cnt[i] - 1] * dp[i - 1][j] ) % MOD; //从cnt[i]个数个选奇数个，C(n，1)+C(n,3)+...=2^(n-1)
                dp[i][j] = ( ( ll ) dp[i][j] + ( ll ) _2p[cnt[i] - 1] * dp[i - 1][j] ) % MOD; //从cnt[i]个数个选偶数个,C(n,0)+C(n,2)*C(n,4)+...=2^(n-1)
            }
        }
    }
    cout << ( dp[70][0] - 1 ) % MOD << endl; //减去0的情况
    return 0;
}

http://codeforces.com/contest/895/problem/D

题意：给你二个等长 ( len <= 1e6 ) 字符串 a, b, 并且 a 的 字典序比 b 小， 让你构造字符串 c , 问 c 的情况有多少种 (mod 1e9+7)

要求：1. c 的字母组成 和 a 一样

2. a 的字典序比 c 小

3. c 的字典序比 b 小

题解：由于题目要求的是 c 的组合有多少种，那可以先求出 c 的排列，再除以每个字母个数的阶乘积即可

求比字符串 b 字典序小的排列个数,

dfs( id, b) 表示 c[ 1~(id-1) ] = b [ 1~(id-1) ] 情况下，比 b 字典序小的有多少

那么 若 c[id] = b[id] 取和字符 b[id] 相等的一个，问题就转化成 dfs(id+1, b)

若 c[id] < b[id] 取比字符 b[id] 小的一个，剩下的字符全排列

若 c[id] > b[id] 则 c 的字典序比 b 的大，舍去

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#define mod 1000000007
using namespace std;
const int MAXN = 1000000+10;
long long A[MAXN];
int num[26];
void init()
{
    A[0] = 1;
    for(int i = 1; i < MAXN; ++i)
        A[i] = A[i-1]*i%mod;
    memset(num, 0x00, sizeof(num));
}
int sum(int k)
{
    int ans = 0;
    for(int i = 0; i < k; ++i)
        ans += num[i];
    return ans;
}
long long dfs(int id, string &s) //使用引用避免爆栈
{
    if( id >= s.size() ) return 0;
    long long ans = 0;
    if(num[s[id]-'a'])
    {
        --num[s[id]-'a'];
        ans = dfs(id+1, s) * (++num[s[id]-'a']) % mod;
    }
    long long k = sum(s[id]-'a');
    return (ans + k * A[ s.size()-id-1 ] % mod) % mod;
}
long long sta[MAXN];
long long silimar()
{
    long long ans = 1;
    for(int i = 0; i < 26; ++i)
        ans = ans * A[num[i]] % mod;
    return ans;
}
long long ex(long long x, long long n) //return  x^n
{
    long long sum = 1;
    x %= mod;
    while(n)
    {
        if( n&1 ) sum *= x, sum %= mod;
        n /= 2;
        x *= x;
        x %= mod;
    }
    return sum;
}
int main (void)
{
    init();
    string a, b, c;
    cin >> a >> b;
    for(int i = 0; i < a.size(); ++i)
        ++num[ a[i]-'a' ];
   long long ans = dfs(0, b) - dfs(0, a); //递归爆栈 MLE
    long long k = silimar();
    ans = ans * ex(k, mod-2) % mod;
    cout << ( ans - 1 + mod ) % mod;
    return 0;
}

题目链接 Eyes Closed

题意 两个人玩一个游戏，现在有两种操作：

1、两个人格子挑选一个区间，保证两个的区间不相交。在这两个区间里面各选出一个数，交换这两个数。

2、挑选一个区间，求这个区间的和的期望。

对于第一种操作，先求出两个区间的长度len1len1和len2len2，再求出两个区间的期望和s1s1和s2s2。

对于第一个区间，我们先把这个区间里的所有数（期望值）乘上(len1−1)/(len1)(len1−1)/(len1)，再加上s2/len1/len2s2/len1/len2

对于第二个区间，我们先把这个区间里的所有数（期望值）乘上(len2−1)/(len2)(len2−1)/(len2)，再加上s1/len1/len2s1/len1/len2

线段树维护这两个操作即可。

#include <bits/stdc++.h>
 
using namespace std;
 
#define rep(i, a, b)    for (int i(a); i <= (b); ++i)
#define dec(i, a, b)    for (int i(a); i >= (b); --i)
#define ls      (i << 1)
#define rs      (i << 1 | 1)
#define mid     ((L + R) >> 1)
#define lson        i << 1, L, mid
#define rson        i << 1 | 1, mid + 1, R
 
typedef long long LL;
 
const int N = 4e5 + 10;
 
double add[N], mul[N], s[N];
int n, q;
 
void pushup(int i){ s[i] = s[ls] + s[rs]; }
 
void pushdown(int i, int L, int R){
    s[ls] = mul[i] * s[ls] + add[i] * (mid - L + 1);
    mul[ls] *= mul[i];
    add[ls] = mul[i] * add[ls] + add[i];
        s[rs] = mul[i] * s[rs] + add[i] * (R - mid);
    mul[rs] *= mul[i];
    add[rs] = mul[i] * add[rs] + add[i];
    add[i] = 0;
    mul[i] = 1.0;
}  
 
void build(int i, int L, int R){
    add[i] = 0, mul[i] = 1;
    s[i] = 0;
    if (L == R){ scanf("%lf", s + i); return; }
    build(lson);
    build(rson);
    pushup(i);
}
 
void add_update(int i, int L, int R, int l, int r, double val){
    if (l <= L && R <= r){
        s[i] += (R - L + 1) * val;
        add[i] += val;
        return;
    }
 
    pushdown(i, L, R);
    if (l <= mid) add_update(lson, l, r, val);
    if (r  > mid) add_update(rson, l, r, val);
    pushup(i);
}
 
void mul_update(int i, int L, int R, int l, int r, double val){
    if (l <= L && R <= r){
        s[i] *= val;
        mul[i] *= val;
        add[i] *= val;
        return;
    }
 
    pushdown(i, L, R);
    if (l <= mid) mul_update(lson, l, r, val);
    if (r  > mid) mul_update(rson, l, r, val);
    pushup(i);
}
 
double query(int i, int L, int R, int l, int r){
    if (l <= L && R <= r) return s[i];
    double ret = 0;
    pushdown(i, L, R);
        if (l <= mid) ret += query(lson, l, r);
    if (r  > mid) ret += query(rson, l, r);
    pushup(i);
    return ret;
}  
 
int main(){
 
    scanf("%d%d", &n, &q);
    build(1, 1, n);
 
    while (q--){
        int op;
        scanf("%d", &op);
        if (op == 1){
            int l1, r1, l2, r2;
            scanf("%d%d%d%d", &l1, &r1, &l2, &r2);
            double s1 = query(1, 1, n, l1, r1);
            double s2 = query(1, 1, n, l2, r2);
            double len1 = r1 - l1 + 1;
            double len2 = r2 - l2 + 1;
            mul_update(1, 1, n, l1, r1, 1.0 * (len1 - 1) / len1);
            mul_update(1, 1, n, l2, r2, 1.0 * (len2 - 1) / len2);
            add_update(1, 1, n, l1, r1, 1.0 / len1 / len2 * s2);
            add_update(1, 1, n, l2, r2, 1.0 / len1 / len2 * s1);
        }
 
        else{
            int l, r;
            scanf("%d%d", &l, &r);
            printf("%.12f\n", query(1, 1, n, l, r));
        }
    }
 
    return 0;
}